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SUMMARY 


A computer code has been developed for analyzing two-dimensional flow fields in 
supersonic combustion ramjet^(scramjet) inlets* The code is written in the CYBER 200 
FORTRAN language for the CDC CYBER 203 vector-processing computer system. It solves 
the two-dimensional Euler or Navier-Stokes equations in conservation form by an 
unsplit, explicit, two-step finite -difference method. A more recent explicit- 
implicit, two-step scheme has also been incorporated in the code for the viscous flow 
analysis. An algebraic, two-layer eddy-viscosity model is used for the turbulent 
flow calculations. Ihe code can analyze both inviscid and viscous flows with no 
strut, one strut, or multiple struts embedded in the flow field. Although the code 
is primarily written for supersonic internal flow, it can be used for a variety of 
other flow problems provided that suitable modifications are made to the boundary 
conditions. This report contains a listing of the program, descriptions of the input 
and output parameters, and a sample flow-field solution. 


INTRODUCTION 

A two-dimensional computer code NASCRIN (^merical ^alysis of SCRa mjet INl et) 
has been developed to analyze scramjet-inlet flow fields. The code is based on the 
analysis developed in references 1 and 2* The analysis as such uses the two- 
dimensional Euler or Navier-Stokes equations in conservation form to describe the 
inlet flow. A two-layer eddy-viscosity model due to Baldwin and Lomax (ref. 3) is 
used for the turbulent flow calculation. In order to facilitate the use of a general 
inlet geometry with embedded bodies, an algebraic coordinate transformation is used 
which generates a set of boundary-fitted curvilinear coordinates (ref. 4). It trans- 
forms the physical domain into a rectangular domain with uniform mesh spacing, and 
embedded bodies in the flow field are transformed into slits (fig. 1). The transfor- 
mation also allows for concentrating mesh lines in regions of high gradients. The 
transformed governing equations are solved by an unsplit, explicit, two-step finite- 
difference method due to MacCormack (ref. 5). This method is highly efficient on 
vector-processing computers, as it allows a high degree of vectorization. A more 
recent explicit-implicit scheme (refs. 6 and 7) has also been incorporated in the 
code for the viscous flow analysis and can be used upon the option of the user. As 
discussed in references 1 and 2, the code can also be used in a quasi-three- 
dimensional sense for a class of scramjet inlets under certain simplifying 
assumptions. 

The code is written in the CYBER 200 FORTRAN language for the CDC® CYBER 203 
vector-processing computer system at the Langley Research Center. In its present 
form, it can analyze two-dimensional inviscid and viscous (laminar and turbulent) 
flows with no strut, one strut, or multiple struts in the flow field. Although the 
code is primarily written for supersonic internal flow, it can be used for a variety 
of other flow problems, including subsonic inflow with supersonic outflow, supersonic 
inflow with subsonic outflow, or fully subsonic flow, by suitably modifying the 
boundary condition subroutines. 

This report contains a brief description of the code structure, various input 
parameters, and output flow quantities. A listing of the program can be found in 
appendix A. A sample case is also presented to illustrate the use of the code. 



CODE DESCRIPTION 


The code listing given in appendix A is primarily set up for supersonic flow, 
although by making suitable changes in the boundary condition subroutines, it can be 
used for supersonic-subsonic or fully subsonic flow. For the supersonic flow case, 
the conditions at the inflow boundary are held fixed at some known values, and 
extrapolation is used at the outflow boundary. The code assumes a solid top boundary 
from the second grid point to the last grid point in the axial direction. In the 
case of a symmetric problem, the code solves the upper half of the flow field. ^ (See 
examples in appendix B.) 

The code consists of a main program, NASCRIN, and 20 subroutines. It has been 
structured in such a way that the user has to make changes only in the first few 
subroutines for a particular problem. A brief description of program NASCRIN and the 
various subroutines is given here in the order they appear in the code. 


Program NASCRIN 

Most input for parameters such as type of geometry, flow, and grid is done in 
NASCRIN. Each of these input parameters is discussed below. SI units are used 
throughout the code. 


Line no. 


Parameter Description 


NASCRIN 52 N1 

NASCRIN 53 Ml 

NASCRIN 57 SWEEP 

NASCRIN 59 NSTRUT 

NASCRIN 62 NS1 

NASCRIN 63 NS2 

NASCRIN 64 MS 

NASCRIN 69 LSYM 

NASCRIN 73 NCWL 


Number of grid points in x-coordinate direction. It is 
also the second dimension in various arrays. 

Number of grid points in y-coordinate direction. It is 
also the first dimension in various arrays. 

Sidewall sweep angle, A (fig. 1). If the code is to be 
used in a quasi-three-dimensional sense, define SWEEP in 
degrees; otherwise, set it to be zero. 

NSTRUT = 0 for no strut embedded in the flow field; 

NSTRUT = 1 for a strut embedded in the flow field. 

If NSTRUT = 0, NS1 , NS2, and MS are not required. For 
NSTRUT = 1 , NS1 and NS2 are grid points in the x-direction 
between which the strut lies. MS is the grid point in the 
y-direction through which an axial mesh line passes that is 
coincident with the top surface of the strut. (See exam- 
ples in appendix B. ) 

A parameter for specifying the type of lower boundary. If 
the flow is symmetric about the lower boundary and only the 
upper half of the flow is being calculated, set LSYM = 1; 
otherwise, LSVM = 0. 

If the lower boundary is fully or partially solid, set 
NCWL =1. If it is a free boundary, set NCWL = 0. 
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Line no, 


Parameter Description 


NAS GRIN 74 NCWL1 

NASCRIN 75 NCWL2 

NASCRIN 1 1 7 NFLOW 

NASCRIN 123 LFI 

NASCRIN 1 24 LFO 

NASCRIN 127 ILT 

NASCRIN 132 LMAX 

NASCRIN 133 LW 

NASCRIN 136 LEXIT 

NASCRIN 1 39 CRIT 


If NCWL = 0, NCWL1 and NCWL2 are not required. If 
NCWL = 1 , NCWL1 and NCWL2 are the grid points in the 
x-direction between which there is a solid lower boundary 
(cowl, strut). 

NFLOW = 0 for inviscid flow; NFLOW = 1 for laminar flow; 
and NFLOW = 2 for turbulent flow. If NFLOW = 2 , it is 
necessary to modify calls to the subroutine EDDY, which 
calculates eddy viscosity. These calls are made in 
subroutine TURB. 

LFI is the unit number from which the solution is read to 
restart the program. It is used only for ILT = 2, and its 
name must be defined on the program card. 

LFO is the unit number on which the output of the current 
run is saved. It also must be defined on the program card. 

ILT = 1 for the first run of the computer program; 

ILT = 2 for second and subsequent restarts. 

LMAX is the maximum number of time-steps the current pro- 
gram will run. The program will stop after LMAX time- 
steps, print the complete solution, and save it on LFO. 

The following message is printed at the end of the 
solution: "PROGRAM TERMINATED ON LMAX." 

LW is the number of time-steps after which an intermediate 
solution is printed. LMAX should be equal to or an integer 
multiple of LW. 

A parameter used to apply boundary conditions at the out- 
flow boundary. LEXIT = 0 for zeroth-order extrapolation; 
LEXIT = 1 for first-order extrapolation. 

Local convergence criterion. Convergence is defined as 

- New density - old density 1 1 

Old density FDT LSTEP 

where LSTEP = Number of time-steps after which the conver- 
gence is checked. In the present code, LSTEP = 2. That 
is, after every two time-steps, quantity ERROR is calcu- 
lated at each grid point, and if it is less than CRIT at 
all the grid points, the solution is assumed converged 
based on the local convergence criterion. The program then 
stops, prints the complete solution, and saves it on LFO. 
The following message is printed at the end of the solu- 
tion: "PROGRAM TERMINATED ON LOCAL ERROR CRITERION." FDT 

is defined on line NASCRIN 219. 


3 



Line no» Parameter Description 

NASCRIN 140 CRITAVG Average or global convergence criterion* Average conver- 
gence is defined as 

ERRORAVG = r- ; 

Number of grid points 

If ERRORAVG becomes less than CRITAVG, the program again 
stops based on the average error criterion. It prints the 
solution and saves it on LFO. The following message is 
printed at the end of the solution: "PROGRAM TERMINATED ON 

MEAN SORT ERROR CRITERION." CRITAVG should be one or two 
orders of magnitude smaller than CRIT. 

NASCRIN 144 FM These are the Mach number, pressure, and temperature, 

NASCRIN 145 PF respectively, at the inflow boundary. SI units are used. 

NASCRIN 1 46 TF 

NASCRIN 170 IMET IMET = 1 if transformation metric coefficients are to be 

printed at each grid point; otherwise, IMET =0. 

NASCRIN 171 BETA BETA controls grid clustering near upper and lower bound- 

aries. Ihe program uses a stretching function that clus- 
ters the mesh lines equally near the boundaries. The 
stretching function can be changed by making changes in 
subroutine XY. BETA should be greater than 1. Ihe closer 
it gets to 1, the higher the stretching (i.e., more points 
are clustered near the boundaries). 

NASCRIN 217 CCIMPY Ihis parameter controls the artificial damping in the 

implicit stage of the ejq>licit-implicit method (refs. 6 
and 7). CCIMPY = 0 turns off the damping, and CCIMPY = 1 
turns it on. 

NASCRIN 218 FDTL A parameter that allows the use of local time-step marching 

to speed up the convergence to steady state. If the local 
time-step is not used, set FDTL = 1 . If the local time- 
step is to be used, set FDTL greater than 1 . For example, 
if FDTL = 2, the time-step at a grid point cannot exceed 
twice the global minimum time-step. If at a grid point, 
the local time-step is less than two times the global mini- 
mum, the local time-step is used; otherwise, the time-step 
at the grid point is set equal to two times the global 
minimum. 


NASCRIN 219 FDT 


FDT is a constant multiplier used in the time-step calcula- 
tions. In the present code. 



over the entire grid. Here Ay is the grid spacing in 
the y-coordinate direction, v is the velocity in the 
y-direction, and a is the speed of sound at a grid point. 
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Line no. 


Parame ter Description 


NASCRIN 219 FDT 
(cont . ) 


NASCRIN 220 CCP 

NASCRIN 223 CCT 


NASCRIN 224 LKK 


NASCRIN 237 TIMED 


Set FDT < 1 for the fully explicit scheme. For the case 
of viscous flow calculations, if FDT > 1, the explicit- 
implicit scheme of MacCormack is turned on automatically. 

Coefficients in the artificial damping e>^ression in the 
explicit stage. The damping expression has two parts, one 
of which depends on the pressure gradient and the other on 
the temperature gradient. CCP multiplies the first part, 
and CCT multiplies the second part. Typical values of CCP 
and CCT are in the neighborhood of 0.5 and should be kept 
as small as possible to keep the artificial damping contri- 
bution to a minimum. 

This parameter allows the user to try various values of FDT 
or CCP and CCT. For LKK = 2, only one set of values can 
be used. But, for LKK = 3, the user can specify two 
values of FDT or CCP and CCT. If the program detects a 
negative temperature in the flow field for the first value, 
it will restart the calculations with the second value. 

But, if the program proceeds normally with the first value, 
it will ignore the second value. (See appendix C for 
examples in using LKK. ) 

TIMED is the time required by the free stream to travel 
3 times the length of the flow domain. For supersonic 
flows, it has been observed that the flow reaches steady 
state if the calculations have advanced in time to TIMED. 

If the program does not stop because of local or average 
convergence criterion and LMAX is sufficiently large, then 
it will stop if the physical time TIME > TIMED. Physical 
time TIME is obtained by summing up the global minimum 
time-step size from each time-step. The program again 
prints the solution and saves it also on LFO. The 
following message is printed at the end of the solution: 
"PROGRAM TERMINATED ON PHYSICAL TIME CONVERGENCE 
CRITERION." 


Subroutine START 

The X, y coordinates of the upper and lower boundaries are prescribed in subrou- 
tine START between lines START 27 and START 78. For example, if the strut EFGH is 
removed from figure 1(c), define the x,y coordinates of lines APQRD and BC between 
the above-mentioned lines in START. However, if there is a strut in the flow field 
as shown in figure 1(c), define the x, y coordinates of lines APQRD and UEHGV between 
lines START 27 and START 61, and define the x, y coordinates of lines UEFGV and BC 
between lines START 64 and START 78. Grid points on lines UE and GV for lower and 
upper portions should have the same x, y coordinates. All the distances should be 
defined in meters. 

For a quasi-three-dimensional application of the code in scramjet-inlet calcula- 
tions, geometry should be prescribed in a plane normal to the sidewall sweep. 
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Subroutine TURB 


Subroutine TURB is called in subroutine VISCOUS at line VISCOUS 222 for turbu- 
lent flow calculations (i.e., when NFLOW = 2). Subroutine TURB, in turn, calls sub- 
routines VORT and EDDY to calculate eddy viscosity. It may be necessary to make 
several calls to subroutine EDDY in order to cover the entire flow field. (See exam- 
ples in appendix B. ) 

There are six parameters in the call to subroutine EDDY. These parameters are 
described below; 


NNX1,NNX2 
NY1, NY2 


NWALL 


NSTR 


Grid numbers in the x, y coordinate directions, respec- 
tively, for the flow domain being considered in a particu- 
lar call to subroutine EDDY. 

NWALL = 1 when there is only one solid boundary (i.e., a 
wall) at NY2; NWALL = 2 when there are solid boundaries 
at both NY1 and NY2. 

NSTR = 1 for the region in which the lower strut surface 
forms the upper solid boundary of the flow domain. In this 
case, set NY2 = MSI ; otherwise, NSTR = 0. 


Subroutine XY 

Subroutine XY performs the algebraic coordinate transformation and stretches the 
grid toward the upper and lower boundaries. Transformation metric coefficients are 
also calculated in this subroutine. 


Subroutines VISCOUS and INVICID 

Subroutine VISCOUS is called at line NASCRIN 241 for viscous (laminar and turbu- 
lent) flow calculations, and subroutine INVICID is called at line NASCRIN 240 for 
inviscid flow calculations. Both subroutines generate a starting solution for the 
first run of the program (ILT = 1) and read the solution from file LFI for second and 
subsequent runs (ILT = 2), Subroutines for flux vectors, damping, and boundary con- 
ditions are called in these subroutines, and the finite differencing for predictor 
and corrector steps is also done here. After the program is terminated either on a 
convergence criterion or on LMAX, subroutine PRINT is called to print the solution, 
and the solution is saved on file LFO, 

In the case of subroutine VISCOUS, if FDT > 1, subroutines IMPY and VECIM2 are 
called to apply the implicit stage of the explicit-implicit method. This method is 
applied only in the y-coordinate direction because for a viscous flow calculation, 
the grid is highly compressed in the y-coordinate direction near the boundaries, 
which makes the minimum grid spacing in the y-direction much smaller than the minimum 
grid spacing in the x-direction. Because of this large difference in spacings, even 
if an implicit stage was incorporated in the x-direction, it would not be used. 

Sometimes it is necessary to use a small value of FDT or a large value of damp- 
ing coefficients CCP and CCT to keep the calculations stable in the initial stages, 
probably because the starting solution is a poor guess. But as the calculations 
progress, it is possible to gradually increase the time-step by increasing FDT or to 
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decrease damping by decreasing CCP and CCT to a desired value over a number of time- 
steps* In the present program, these changes in FDT or CCP and CCT can be incorpo- 
rated by inserting a few cards after line VISCOUS 51 and after line INVICID 60* 
Several examples of implementing these changes are given in appendix D. 


Subroutines DVISCP and DVISCC 

Subroutines DVISCP and DVISCC are called only for viscous flow calculations. 
Both subroutines calculate the viscous dissipation terms in the Navier-Stokes equa- 
tions. Subroutine DVISCP uses backward differencing, and subroutine DIVSCC uses 
forward differencing in the calculation of these terms. 


Subroutines VEC1 and IVEC1 

Subroutines VEC1 and IVEC1 calculate the flux vectors for viscous and inviscid 
flow, respectively. 


Subroutines DAMPP and DAMPC 

Subroutines DAMPP and DAMPC are used to add artificial damping in the flux vec- 
tors. Subroutine DAMPP uses backward differencing, and subroutine DAMPC uses forward 
differencing in the calculation of damping terms. 


Subroutine VEC2 

Subroutine VEC2 is called after the predictor or corrector step. It calculates 
local flow quantities from the conserved variables. 


Subroutine BOUND 

Subroutine BOUND applies the boundary conditions for the viscous flow. Ihe 
subroutine is set up for supersonic inflow and outflow boundaries. It applies a 
no-slip boundary condition on the upper boundary from the second grid point to the 
last grid point. The lower boundary can be either a solid wall, a free boundary, or 
a symmetry line. Ihe conditions are held fixed at some prescribed values at the 
inflow boundary, and zeroth- or first-order extrapolation is used at the outflow 
boundary. 

To use the program for supersonic inflow and subsonic outflow, the user needs to 
prescribe a back pressure at the outflow boundary. This can be done by changing line 
BOUND 76 for LEXIT = 1 or line BOUND 82 for LEXIT = 0 to the line 
P(1,N1;M1) = PB where PB is the back pressure in pascals. 

It is also possible to modify the code for a subsonic inflow boundary. The user 
should insert the proper boundary conditions for subsonic inflow after line BOUND 33. 
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Subroutine IBOUND 


Subroutine IBOUND applies the boundary conditions for the inviscid flow. This 
subroutine is also set up for supersonic inflow and outflow boundaries/ although 
changes similar to those suggested in subroutine BOUND can be made to modify it to 
handle a subsonic inflow or subsonic outflow boundary. 

Both subroutines BOUND and IBOUND check for negative temperature in the field 
during each predictor and corrector step. Anytime a negative temperature occurs in 
the field, the code prints the following message; "NEGATIVE TEMPERATURE IN THE 
FIELD." It then prints the number of time-steps at which the negative temperature 
occurred and the values of all the flow variables at each grid point. This allows 
the user to examine the flow variables in the region where the calculations became 
unstable. After this, the code returns to statement 1000 at line NASCRIN 221. It 
either restarts with a new value of time-step or damping parameters, or stops. 


Subroutines VORT and EDDY 

Subroutines VORT and EDDY calculate the eddy viscosity from the algebraic turbu- 
lence model of reference 3. Subroutine VORT calculates the vorticity at all the grid 
points. The vorticity is then used in subroutine EDDY to calculate the eddy viscos- 
ity. The parameters in subroutine EDDY were discussed earlier in subroutine TURB. 


Subroutines IMPY and VECIM2 

Subroutines IMPY and VECIM2 are used only for the viscous flow calculations to 
apply the implicit stage in the y-direction of the explicit-implicit method. Subrou- 
tine IMPY actually applies the implicit stage, and then subroutine VECIM2 calculates 
the flow variables at each grid point from the conserved variables, calls subroutine 
BOUND, and saves the flux quantities on the appropriate upper or lower boundary for 
the implicit stage of the next time-step. The implicit stage is used only when 
FDT > 1 • 


Subroutine SPILL 

Subroutine SPILL is called after every LW time-steps at line VISCOUS 244 for the 
viscous flow and at line INVICID 172 for the inviscid flow. The purpose of this 
subroutine is to calculate the spillage in a quasi-three-dimensional application of 
the code. It takes the solution in the plane normal to the sidewall sweep and pro- 
jects it back to the cowl plane. The constant tangential velocity component is then 
properly superimposed on the projected solution to provide the solution in the cowl 
plane. Knowing the component of velocity normal to the plane of the cowl and the 
density field, flow spillage (i.e., the amount of flow leaving the inlet ahead of the 
cowl) can be obtained. For viscous flow calculations, an arbitrary boundary-layer 
profile is prescribed on the constant tangential velocity component before superim- 
posing it on the projected solution. In subroutine SPILL, it is assumed that the 
tangential velocity changes from zero to the constant inviscid value over 10 to 
12 grid points from the solid boundary. For the case of zero sweep, obviously the 
spillage in the quasi-three-dimensional application of the code will be zero. For 
this case, the subroutine will skip the spillage calculations. 
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Subroutine SPILL also calculates and prints the nondimensional pressure on the 
upper and lower boundaries and on the upper and lower strut surfaces if there is a 
strut in the flow field. Ihe pressure is nondimensionalized with respect to the 
pressure at the inflow boundary. 


Subroutine PRINT 

Subroutine PRINT is called after every LW time-steps and after the normal termi- 
nation of the program. It prints the coitplete flow-field solution at all the grid 
points. It also prints certain input parameters and the mass flux across a given 
axial location. 


SAMPLE PROBLEM AND OUTPUT QUANTITIES 


To illustrate the use of the code, the listing given in appendix A is set up for 
a quasi-three-dimensional application of the code for the inlet geometry shown in 
figure 1 • The geometry input in subroutine START is in a plane normal to the sweep 
line (fig. 1(c)). Only the upper half of the problem is solved. The problem is set 
up to calculate turbulent flow at Mach 5. The following values of input parameters 
have been used: 

N1 = 55 

Ml = 61 

SWEEP =30 

NSTRUT = 1 

NSl =16 

NS2 = 46 

MS = 29 

LSYM = 1 


NCWL = 1 
NCWLl = 24 
NCWL2 = 44 
NFLOW = 2 
LFI = 

LFO = sJ 
ILT = 1 
LMAX = 25000 


(The program card equates UNIT9 to CASE20 at line 
NASCRIN 2. CASE20 is the name of the file on which 
the solution is saved.) 
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LW = LMAX 


LEXIT = 1 
GRIT = 0,0001 


CRITAVG =0.1 * GRIT 

FM = 5. (FM is the Mach number at the face of the inlet in 

the plane normal to the sweep line. ) 

PF = 10000. 

TF = 300. 

IMET = 1 
BETA = 1.005 


FDTL = 1. 

FDT = 1 

GGP = 0.5 

GGT = GGP 

TIMEL = 3. * x(1,N1)/UF (UF is the velocity at the face of the inlet 

corresponding to Mach 5. ) 

To calculate eddy viscosity, five calls to subroutine EDDY are made in subroutine 
TURB as shown in the listing in appendix A. As mentioned earlier, for the quasi- 
three-dimensional application of the code, the flow is calculated in a plane normal 
to the sidewall sweep. The code prints the detailed solution in this plane. It then 
calculates the spillage from the inlet by projecting the solution to a plane parallel 
to the cowl plane and properly superimposing on it the component of velocity parallel 
to sweep line. The sample outputs are presented in appendix E, and a brief descrip- 
tion of output quantities is given below. 

Since parameter IMET is set equal to 1, the code prints the metric coefficients 
at all the grid points. In this portion of the output, 

grid number in x-direction 

grid number in y-direction 

distance in x-direction, m 

distance in y-direction, m 

_ 

■ ac 
_ ^ 


N 

M 

X 

Y 

XXI 
YXI 
1 0 



XETA 


5x 

5r| 


YETA 


AJ 


_ ^ 

5ti 

” 5ti “ 5t] 


where ^,r] are the coordinate directions in the transformed plane. 

The code then prints the maximum error and mean square root error after every 
50 time-steps till it stops either on LMAX or on one of the convergence criteria. It 
then prints the total number of time-steps (or iterations), the physical time 
advanced in seconds, and the time-limit in seconds. (Time-limit is set by the param- 
eter TIMEL. ) It also prints the values of CRIT (local error), CRITAVG (average 
error), and LEXIT (order of extrapolation at outflow)* The complete solution is then 
printed at each body station in the x-direction. It includes 

X,Y coordinates of the grid point, m 

U,V axial and normal velocities, m/sec 

P pressure. Pa 

T temperature, K 

RO density, kg/m^ 

SH static enthalpy, m^/sec^ 

VIST turbulent viscosity, kg/m-sec 

ERRO local error at a grid point 

For NSTRUT = 1 , flow variables at M = MS between N = NS1 to NS2 correspond to 
the top surface of the strut. Flow variables at the lower surface are printed 
separately. At the lower surface of the strut, the code prints total enthalpy HS and 
laminar viscosity VISS in place of VIST and ERRO. 

With each call to subroutine PRINT, after printing the detailed solution, the 
code prints the mass flow rate at each axial station. This allows the user to check 
for mass conservation inside the inlet. 

After the code stops on LTiAX or any of the convergence criteria, it calls sub- 
routine SPILL, which calculates flow spillage. If the sweep angle is zero, spillage 
calculations are not done. A sample output for spillage calculations is shown in 
appendix E. In this output, XP is the projected axial distance in the cowl place 
(i.e., XP = X/cos A). Local mass spilled (in kilograms per meter-second) is the 
integrated value of the product of density and downward velocity component at a given 
XP location. Total mass spilled (in kilograms per second) is obtained by integrating 
the local values up to a given XP location. If, for example, the cowl closure starts 
at XP = 0.132791 m, corresponding to body station 40, the total mass spilled ahead 
of the cowl is equal to 0.01193 kg/sec. The user can then calculate the spillage in 



terms of the percentage of the incoming mass flow rate at the face of the three- 
dimensional inlet. 

Subroutine SPILL also calculates and prints the surface pressure distribution on 
the top and bottom (or centerline in the case of a symmetric problem) boundaries and 
on the upper and lower strut surfaces for NSTRUT = 1. Ihese pressures are nondimen- 
sionalized with respect to the pressure at the face of the inlet. 

Finally, the code prints a message that tells on what basis the calculations 
were terminated. 

From the solution of the sample problem presented here, pressure contours and 
the velocity vector field were plotted. These plots are shown in figures 2 and 3, 
respectively. Although only the upper half of the inlet was solved, the results were 
projected to the lower half of the inlet for plotting purposes. 


COMPUTATIONAL RATE AND TIME REQUIREMENTS 

The code is operational on the CDC CYBER 203 computer system and is highly vec- 
torized to take full advantage of the vector-processing capability of the system. 
With the explicit method, it has a computational rate of 2 x 10“^ second per grid 
point per time-step for viscous flow and 1.3 x 10 ^ second per grid point per time- 
step for inviscid flow when the 64-bit word arithmetic is used. The code can easily 
be modified to use the 32— bit word arithmetic, in which case, the computational rate 
is improved by a factor of about 2.5. For example, the code can perform approxi- 
mately 50 time— steps per second for viscous flow and 75 time— steps per second for 
inviscid flow on a grid of 51 x 5i when the explicit method is used with the 32-bit 
word arithmetic. An inviscid problem using a 51 x 5i grid and requiring 1000 time- 
steps for convergence can be solved in about 15 seconds. 

Finally, in the case of viscous flow calculations with very fine mesh spacing 
near the solid boundaries, a significant improvement in computational rate can be 
obtained by using the explicit-implicit method. 


CONCLUDING REMARKS 

A vectorized computer program NASCRIN has been developed that solves the two- 
dimensional Euler or Navier-Stokes equations in conservation form. It primarily uses 
MacCormack's explicit method, but for viscous flow calculations, MacCormack's 
explicit-implicit method can also be used at the user's discretion. An algebraic, 
two-layer eddy- viscosity model is used for the turbulent flow calculations. The code 
can analyze problems with or without embedded bodies in the flow field. It is highly 
user oriented and is structured in such a way that for most supersonic flow problems, 
the user has to make changes only in the main program and first two subroutines. 
Although the code is primarily written for supersonic internal flow, it can be used 
with some modifications for a variety of other flow problems. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 1 , 1 983 
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PROGRAM LISTING 


PRQGRAM MASCRIN( INPUT> OUT PUT ^ TAPE 5«I NPUT^ TAP E6-0UTPUT^ 
1UNIT9-CASE20) 

NASCRIN stands for NUMERICAL ANALYSIS OF SUPERSO 
RAMJET INLETS. 


THIS PROGRAM IS BASED ON AIaA PAPER NO. 61-105^ 
ANU NASA TP 1934# ALL BY AJAY KuHAR OF NASA LANG 

THE PkOGRaM solves THE TWU-D I MENS I ONAL NAvIER-ST 
LAMINAR AND TURBULENT FLOW AND T wO-DI MENS lONAL E 
INVISCIO FLOW IN FULLY CONSERVATIVE FORM. 


NASA TP 1940# 

LEY RESEARCH CENTER. 


OKES 

ULER 


EQUATIONS 

EQUATIONS 


IT USES MaCCORMaCK'S JNSPLIT# EXPLICIT FINITE-DIFFERENCE METHOD. 

FOR VISCOUS FLOW# MACCORMACK'S EX PL IC IT- 1 MPL I C IT SCHEME CAN ALSO BE 
USED at USER'S OPTION. 

TURBULENCE IS MODELED BY B ALD WlN-L DMA X MODEL. 

AN algebraic COORDINATE TRANSFORMATION IS USED TO GENERATE A SET OF 
BOUNDARY FITTED CURVILINEAR COORDINATES. 

OPERATIONAL ON C DC-C Y B E K-203 VECTOR PROCESSING COMPUTER. 

BIT 81 

COMMON/ FI/ NX M#NXM1#NXM2#NXM3#NXM4#N1#M1#N11# 
1#NCWL#NCWL1#NCWL2#NCWLM#NC»LP#NSTRUT 


A1 

/M 

12#M 

13 

# N12 # N1 3 

# 

NXM5,N 

XM6 

#MS2 

#MS3 

#M 

S4,iSPl 

#MS 

P2 

IP 

1# 

NS2M 

1# 

NS1M1#NS 

2P1# 

NNS 

#NS 

1M2# 

NXM7 






2/ 

RO 

(61# 

55 

)#U(61#5 

a 

> 

) #v 

(61 

#55 

)#F( 

61#5 

5) 

#T(6 

1#55)# 


5) 

#H 

(61# 

55 

)#RJS(55 

) 

#US 

(55 

)# V 

S(55 

)# 






TS 

(5 

5)#S 

HS 

(55)#HS( 

5 

5) 










3/ 

VI 

SL(6 

1# 

55)# VIST 

( 

61# 

55) 

#VI 

S(61 

#55) 

#V 

ISS( 

55) 



4/ 

SI 

GX(6 

1# 

55)#SIGY 

( 

61# 

55) 

#TA 

UXY( 

61#5 

5) 

#QX( 

61# 

55) 

# 

5) 

#s 

IGXS 

(5 

5 )#olGYS 

( 

55) 

>TAUXY 

S(55 

)#QX 

S( 

55), 

QYS 

(55 

) 

5/AU 

(61# 

55 

#4 )# AM(6 

1 

#55 

#4) 

#AN 

(61# 

55#4 

)# 

AUK 

61# 

55# 

4) 


C0MM0N/F6/AUSC5!?,4)#AMS(55#4)# ANS(55#4)# AUSl (55# 4) 

C0rtr1UN/F7/X( 61# 55)#Y(61#55), XXI(61#5 5)#YXI (bl#55)# 

1 XETA(61#55)#YETA(61#55)#AJ( 61#55 ) 

C0MM0N/F8/XS (55)#YS(55)#)(SX1 ( 5 5 ) # YSX I ( 55 ) # XS ET A ( 55 ) # 

1 YSETA(55)#AJS(55) 

CuMM0N/F9/INT(55 )#TE1(55)#1 E2(55)#TE3 (55)#TE4(55) 

C0MM0N/F10/DT(61#55)#0TM 

CaMM0N/Fll/ERR0(61#55)#Bl(61#55) 

C0MM0N/F12/LSYM#NFL0W# LF I # LF 0# ILT # L MAX# L W# LE XIT# CR IT#FDT #FDTL# 
1 TIME#CCP#CCT#L#LK#T1MEL#CR1TAVG 
C0MMaN/F13/GAMA#RGAS#CP#PR#PKT#CVIS# FM# PF# TF # OF# ROF# SHF # HF 
C0MM0N/F14/DAU(61# 55#4)#B0AU(55#4)#BDAuS(55# 4)#T1( 55#5)# 
iCCIiPY#6Ml#GM2#GM3#GM4,GM5 
N1*55 
Ml- 61 

DEFINE SIDEwALL SWEEP. IF THE FLOW IS BEING CALCULATED IN A 
NORMAL TO SIDEWALL SWEEP# DEFU4E SIDEWALL SWEEP IN DEGREES. 
OTHERWISE EQUATE IT TD ZERO. 

SWEEP-30. 

NSTAUT-0 FOP. NO STRUT CASE. NSTRuT-1 FOR A STRUT 
NSTRUT-1 

STRUT IS LOCATED BETWEEN NSl AND NS2 MESH POINTS 
AND XS COINCIDENT WITH MESH MS IN ETA DIRECTION. 

NSl-16 
NS2-46 
MS-29 

LSYM IS A PARAMETER FOR THE LOWER BOUNDARY. 

IF THE LOWER BOUNDARY IS A SYMMETRY LINE AND 
wE ARE SOLVING ONLY UPPER HALF OF THE FLOW DOMAIN# LSYM-1 


NASCRIN 
NASCRIN 

NIC COMBUSTION NASCRIN 

NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
FOR NASCRIN 

FUR NASCRIN 

NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 

NASCRIN 

hll#N51#NS2#NS# Mb#MSl NASCRIN 

NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
PLANENASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 
NASCRIN 


IN THE FLOW. 


IN XI DIRECTION 


1 

2 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 
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C 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


IF THE LOWER BOUNDARY IS NOT A SYflHETRY LINE^ LSYM ■ 0 
LSYrt-1 

NCWLl and NCWL2 ARE THE HESH POINTS BETWEEN WHICH COWL OK THE 
CENTER STRUT IS LOCATED. IF THERE IS NO COWL OR STRUT ON THE 
CENTER LINE> SET NCWL « 0 OTHERWISE SET « 1 . 

NCWL-1 

NCWLl-2^ 

NCWL2-AA 

Nll-Nl-1 

N12-N1-2 

N13-N1-3 

M12-M1-2 

M13-M1-3 

NXri-Ni*Ml 

NXNi-NXM-1 

NAN2-NXN-2 

NXf13-NXM-2*Ml 

NXN^»NXM3-2 

NX.i5-NXM-4^P,l-2 

NXrt6*NXM-2*Hl-4 

nXH7»NXM-A^N1-^ 

NSlHi-NSl-1 

NSiH2«NSl-2 

NSlPl-NSl+1 

NS2H1-NS2-1 

NS2P1-NS2+1 

NNS-Nl-NSl^a 

N:>-NS2-NS1 + 1 

NSill-NS-l 

Mb2*M:>-2 
M33-NS-3 
MS^*MS-4 
HSPi-(lS + l 
MSP2-NS+2 
MSP3T1S + 3 
NCWLM-NCWLl-1 
NCwlP«NCWL2+1 

NFlJw IS A PARAMETER FOR THE TYPE OF FLOW. NFLOW- 0^ 1 OR 2 
OEPlEnDING on WHETHER THE FLuw IS INVISCID> LAMINAR JR TURBULENT. 

IF NFLQW«2> IT IS NECESSARY TO MODIFY THE CALLS TO SUBROUTINE E DuY 
WHICH ARE FOUND IN SUBROUTINE TUKB . 

THE ENTIRE FLOW FIELD MAY HAv'E TU BE DIVIDED INTO SEVERAL SUB- 
DOMAINS AND MORE THAN ONE CALL TO SUBROUTINE EDDY MAY BE REQUIRED. 

SEE EXAMPLES IN USER'S MANUAL FUR FURTHER CLARIFICATION. 


NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

naScrin 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCKIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCKiN 

NASCRIN 

NASCRIN 

NASCRIN 




NFLOW-2 

LFi IS THE PERMANENT FILE iNUMbER FROM nHICH THE SOLUTION IS 
READ TO RESTART THE PROGRAM. LFI IS REQUIRED ONLY WHEN ILT-2. 
LFO IS THE PERMANENT FILE NUMBER ON WHICH THE OUTPUT OF THE 

Current computer run is stored, both lfi and lfo should be 

DEFINED ON THE PROGRAM CARD. 

LFl-9 

LF0«9 

SET ILT-1 FOR THE FIRST RUN OF THE COMPUTER PROGRAM. 

TO RESTART THE PROGRAM, SET ILT-2 
ILT-1 

LMAX IS THF MAXIMUM NUMBER JF TImE-STEPS. THE PROGRAM WILL 
STOP AFTER LMAX TIME-STEPS AND WILL STOkE THE SOLUTION ON LFO. 

AN INTERMEDIATE SOLuTiON IS PRINTED AFTER EVERY Lw TIME-STEPS. 
lMAX should be equal to or AN INTEGER MULTIPLE OF LW. 

LMAX-25000 

LW-LMAX 

LEXlT-0 FOR ZERO OkOER EXTRAPOLATION AnD LEXIT-1 FOR FIRST ORDER 
extrapolation at THE OUTFLOW BOUNDARY. 

LEXIT-1 


NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 


68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
ICO 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 
127 
126 

129 

130 

131 

132 

133 

134 

135 

136 
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c 

ERRJR CRITERION 

NASCRIN 

137 

c 

critavg is the average error over the entire elow field 

NASCKIN 

138 


CRIT-.OOOl 

NASCRIN 

139 


CRITAvG»0.1*CR1T 

NASCRIN 

140 

c 


NASCRIN 

141 

c 

CONDITIONS AT THE FACE OF THE INLET. UiE SI UNITS. 

NASCRIN 

142 

c 

(P IN N/M2 AND T IN K) 

NASCRIN 

143 


Fh«5. 

NASCRIN 

144 


PF-10000. 

NASCRIN 

145 


TF-300. 

NASCRIN 

146 


G Art A ■ 1 . 4 

NASCRIN 

147 


RGAS«287. 

NASCRIN 

146 


C P»GAMA*RGAS/(6ArtA-l. ) 

NASCRIN 

149 


Grtl*G Art A“1 . 

NASCRIN 

150 


Grt2»Grtl/GAMA 

NASCRIN 

151 


Gh3«CP/GAMA 

NASCRIN 

152 


PR-.72 

NASCRIN 

153 


PRT-.9 

NASCRIN 

154 


Gh4«2.*GAMA/PR 

NASCRIN 

155 


GM5«GAMA*RGAS 

NASCRIN 

156 

c 

CVIS IS THE COEFFICIENT USED IN SUTHERLAND'S VISCOSITY FORMULA 

NASCRIN 

157 

c 

FOR AIR. 

NASCRIN 

158 


CVIS-I.A58E-6 

NASCRIN 

159 


UF«FM*{GAMA*RGAS*TF)**.5 

NASCRIN 

160 


ROF-PF/RGAS/TF 

NASCRIN 

161 


SHF-GAMA*RGAS*TF/(GAMA-1. ) 

NASCRIN 

162 


HF«SHF+JF*UF/2. 

NASCRIN 

163 

c 


NASCRIN 

164 

c 

GEOMETRY INPUT AND METRIC COEFFICIENTS CALCULATIONS. 

NASCRIN 

165 

c 


NASCRIN 

166 

c 

DEFINE IMET-0 IF METRIC COEFFICIENTS NEED NOT BE PRINTED. 

NASCRIN 

167 

c 

IMET-1 IF METRIC COEFFICIENTS ARE TO BE PRINTED, 

NASCRIN 

168 

c 

BETA-STRETCHING FACTOR. IT SHOULD BE GREATER THAN 1. 

NASCRIN 

169 


IMET-1 

NASCRIN 

170 


BETA-1.005 

NASCKIN 

171 


MM-0 

NASCRIN 

172 


1F(NSTRUT.EQ.0)MS-M1 

NASCRIN 

173 


CALL START(MM>Si<EEP) 

NASCRIN 

174 


CALL XY(l»rt.S,BETA) 

NASCRIN 

175 


1F(NSTRUT.EQ.0)G0 TO 500 

NASCRIN 

176 


MM-1 

NASCRIN 

177 


DO 10 NN-NS1>NS2 

NASCKIN 

178 


N-NN-NSl+1 

NASCRIN 

179 


XS(N)-X(MS»NN) 

NASCRIN 

180 


YS(N)-Y(MS>NN) 

NASCRIN 

181 


XSETA(N)-XETA(MS>NN) 

NASCRIN 

182 


YSErA(N)-YETA(MS»NN) 

NASCRIN 

183 


XSXI (n)-XXI (MS>NN) 

NASCRIN 

184 

10 

YSXI (N)-YXI (MS»NN) 

NASCRIN 

185 


CALL STARTIMM, SWEEP) 

NASCRIN 

186 


CALL XY(MS»M1»BETA) 

NASCRIN 

187 


DO 20 N-1»NS1M1 

NASCRIN 

188 


YETA(MS»N)- ( YIMS+lfNl-YIMS-lfN) )/2. 

NASCRIN 

189 

20 

XETA(MS*N)-(X(rtS + l»N)-X(MS-l»Nn/2. 

NASCRIN 

190 


DO 30 N-NS2P1>N1 

NASCRIN 

191 


YETA(MS»N)-(Y(MS'-1>N)-Y(MS-1/N) )/2, 

NASCRIN 

192 

30 

XETA(MS» N)-( X(MS + 1>N)-X(MS-1^N) )/2, 

NASCRIN 

193 


AJS(1;NS)-XSXI(1)NS)''YSETA(1}NS)-XSETA(1;NS)*YSXI(1}.NS) 

NASCRIN 

194 

500 

CONTINUE 

NASCRIN 

195 


AJ(1»1JNXH)-XXI( 1»1JNXM)*YETA(1»1JNXM)-XETA( 1»UNXM)*YXI{1» IJNXM) 

NASCRIN 

196 


IF(1MET.EQ.O)GO TO A5 

NASCRIN 

197 


DO AO N-1>N1 

NASCRIN 

196 


WRITE(6»1100) 

NASCRIN 

199 


Ou AO MM-1»M1 

NASCRIN 

200 


M-MI-MM+I 

NASCRIN 

201 

40 

WRirE(6,1200)N»H.X(rt>N)>Y(M»N) ,XX 1 ( M< N) > XETA (M» N )» YXI (M»N ) , YET A ( M» NASCRIN 

202 


1N)>AJ(M»N) 

NASCRIN 

203 


IF(NSTRUT.EQ.O)SO TO A5 

NASCRIN 

204 


M-MS 

NASCRIN 

205 
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WRITt<6>1100) 

i<RITfc(6,1200 ) (M,M,AS(N)>irS(N)»XSXI(N)>X 5 ETA(H)>YSXI(N),YjETA(K)> 
IS (N)/M«1»NS) 

1100 format (/, 3X,'M*>3X» ' rt • / 1 OX/ • X • > 1 0X< ' Y • « 8 X> ' XX I S 7 X/ • XE TA • > 8 X» 
l‘YXI'/7X/ 'YETA'»9X/ 'AJ'»/ ) 

1200 FQRMAT(2X>2I3>7F11.8) 

A5 CaNTINJE 
C 

C INTERVAL VECTOR. 

INT( IjND-OeVISTLd/Mlj INTUiNl )) 

lKK-0 

CCMPY-1. 

EDTL-1. 

FDT-1. 

CCP-.5 

1000 CONTINUE 
LKK-LKK+1 
CCT-CCP 

IF (lKK.EQ.2)ST0P 
LK«1 

00 50 I«1^4 
dDAJi>(l,I;Nl )«D. 

6DAO(l,im)»0. 

OAUdd# i;nxm)-o. 

bO CONdNUE 

VdT dd5NXM)-0. 

ViSLdd JNXM)-0. 

visdd;NXM)»o. 

viss(i;ni)*o. 

C riNEL IS THE TIME kEQJIRED 3Y THE FREE STREAM TO TRAVEL 3 TIMES 

C OVER THE LENGTH OF THE FLOW DOMAIN. 

TlMEL-3.*Xd/Nl)/UF 

TIME-0. 


L-0 

IF(NFLOW. EO.O) 
IF(NFlOW.NF.O) 

I F ( L K . E Q . 2 ) GO 
STOP 
ENu 

:>UaRUUTINE START(MM^SWEEP ) 


CALL INVi:iD(SWEEP ) 
CALL VISCOUS(SWEEP) 
TO 1000 


NASCRIN 

AJNASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

NASCRIN 

START 


CUMMQn/F 1/NXM/NXM1^ NXM2/NXMi/NXMA#NiiMif Nll> M 1> NS 1# NS2 > NS# MS# MSI START 


1>NC WL >NCWLdNC4L2>NC^ LM^NChLP^NSTRUT START 

CUMM0N/FA1/M12# M13#N12^M3/NXM5#NXM6#MS2#MS3#MSA#MSP1#M5P2# MSP3 START 

l#NSMi#NSlPI>N$2Ml.NSlMl#NS2Pl/NNS#NSlM2#NXM7 START 

C0MM0N/F7/X(61,5 5),Y(61,5 5)^ XXI(6lo5),YXI (6 1# 55)# START 

1 XETA(61#55)#YETA(61#55)#AJ(61#55) START 

COMMON/FB/XS ( 55 ) # YS ( 55 ) # X S X I ( 5 5 ) # YSX I ( 55 ) # XS ET A { 5 5 ) # START 

i YSETa(55),AJS(55) START 

FOR MM«0# THIS SUBROUTINE DESCRIBES THE X#Y COORDINATES OF THE START 

LOWER BOUNDARY AND THE MESH LINES PASSING THROUGH THE LOwER SURFACSTART 

OF THE STRUT. FOR MM-1, IT DESCRIBES THE COORDINATES OF THE START 

MESH LINE PASSING THROUGH THE UPPER SUkFaCE OF THE STRUT AND THE START 
UPPER boundary. if THERE d NO STRUT IN THE FLOWFIELD# THEN IT START 
U£:>CRI*BFS THE COORDINATES OF UPPER AND LOWER BOUNDARIES. START 

Pl-4 . *ATAN ( 1 . ) START 

SW«Sw£EPtPI/180. START 

^t*^*^^***r^*****f START 

START 

DESCRIBE X(dN), X(M1,N), Y ( 1 # N )# Y ( Ml # N) # FOR N* 1 TO N1 START 

AND# IF MM-1# DESCRIBE START 

X(MS#N)/ Y(MS#N) START 

START 

FOR CJASi-THREF-DIMENSTIUNAL C A LC ULA TI UNS# PROVIDE THE GEOMETRY IN START 
A PLANE NORMAL TO THE SIDEWALL SWEEP. START 

IF IN INCHES# MULTIPLY BY ♦OZSA TO GET METRIC UNITS START 

START 

THE GEOMETRY GIVEN BELOW IS IN A PLANE NORMAL TO SWEEP LINE. START 

SWEEP ANGLE IS 33 DEGREES. START 


206 

207 

208 

209 

210 
211 
212 
213 
21A 

215 

216 

217 

218 

219 

220 
221 
222 
223 

225 

226 

227 

228 

229 

230 

231 

232 

233 
23A 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


16 



APPENDIX A 




TANC1«TAN(6. ♦PI/180. ) 




START 

30 



TANCZ-TANIT.S+PI/ISO. ) 




START 

31 



TANC3-TA'<(6.271178*PI/160.) 




START 

32 



X(l,l)-0. 




START 

33 



DO 1 N-2»6 




START 

34 


1 

X(l»N)»X(l^N-l)+.005 




start 

35 



X(1,7J»X(1,6)+.0045 




START 

36 



xa»8)-X(l»7) + .004 




START 

37 



X(i»9)«X(l,8 )t-.0035 




START 

38 



X(l»10)»X(l,9)+.003 




START 

39 



00 2 N«11»N1 




START 

40 


2 

X<l,N)-X(l,S-l)+.0025 




START 

41 



DX-.O025 




START 

42 

c 






START 

43 



IFIrtM.EQ.DGO TO 300 




START 

44 

c 






START 

45 



DU 3 N«1»N1 




START 

46 



XIiiSjNl-Xd.N) 




START 

47 



Y( 1»N)«0. 




START 

46 


3 

Y(MS,N)-.01R 




START 

49 



OCJ <» n«25»3A 




START 

50 


4 

Y(1#N)«Y(1,N-1)*DX^TANC2 




START 

51 



00 5 N«35»A3 




START 

52 


5 

Y(1»N)-Y(1,N-1)-0X^TANC2 




START 

53 



DO 6 N-17,3A 




START 

54 


6 

Y(MS»N)-Y(MS»N-1 )-0X*TANC2 




START 

55 



DO 7 N«35>A6 




START 

56 


7 

Y(MS»N)-Y(MS>N-1 )+0X+TANC3 




START 

57 



OY-IYIMS* A6)-.013)/.0225 




START 

58 



00 8 N-A7»N1 




START 

59 


8 

Y((1S/N)«Y<MS»M-1 )-OX^OY 




START 

6C 

c 


«*«»««#*«♦♦♦♦♦♦♦♦♦* 




START 

61 



lF(f1rt.EQ. DIRETURN 




START 

62 

c 


*««♦*♦♦♦♦♦ ♦♦♦♦♦♦♦ 




START 

63 


300 

CONTINUE 




start 

64 



Y((U»1)-.0A 




START 

65 



Y(M1/2)«.0A 




START 

66 



00 10 N»1>N1 




START 

67 



X(M1,n)-X ( 1»N) 




START 

68 



X(f1S»N)«X( 1»N) 




START 

69 


10 

Y(MS>N)-.018 




START 

70 



00 11 N«3»N1 




START 

71 


11 

Y(Ml»N)«Y(m,N-H-(X(Ml»N)-X(Ml>N-l) I^TANCl 




START 

72 



DO 12 N-37.A6 




START 

73 


12 

Y(MS»N)«Y(MS»N-1)-DX%TANC1 




START 

74 



0Y-( Y(MS»46)-.013)/.0225 




START 

75 



00 13 N-A7,N1 




START 

76 


13 

YIMS»N)-Y(MS»N-1 I-DX^uY 




START 

77 

c 






START 

78 



RETURN 




START 

79 



END 




START 

80 



iOBROoTINE TURB 




TURB 

1 



C0MNQN/Fl/NXM»NXm»NXM2/NXM3»NXMA»Nl»Ml»Nll»Kll 

# NS 1# NS2# nS# 

MS»HS1 

TURB 

2 



l>NC.<L/NCWLl>NCi<L2»NC.(t.M»NC«aP»NSTRUT 




TURB 

3 



C0MN0N/FAl/m2»M13#N12»N13f NXM5/NXM6#«i2#MS3>«SA/rtSPl#MSP2/MSP3 

TURB 

4 



l»NSNl>NSlPl,NS2Hl»NSlrtl>NS2Pl>NNS^NSlM2>NXM7 




TURB 

5 



Call vort 




TURB 

6 

c 






TURB 

7 



8 

c 

IF 

NFLDw-2 WAS SPECIFIED# EDDY VISCOSITY IS CALCULATED 

USING CALLS TO 

TURB 

9 

c 

SUBROUTINE EDDY. 




TURB 

10 

c 






TURB 

11 

c 

SUBROUTINE EDDY HAS SIX PARAMETERS! NNX 1,nNX 2/ N Yl> 

NY2 

/NrfALL/LNSTR. 

TURB 

12 

c 


NNX1>NNX2 and NY1>NY2 ARE THE GRID NUMBERS IN THE 

X 

AND Y 0IRECTI0NSTUR6 

13 

c 


RESPECTIVELY FOR THE FLOW DOMAIN BEING CONSIDERED. 


TURB 

14 

c 


NWALL-1 WHEN THERE IS 0>ILY ONE WALL (I.E.# A WALL 

AT 

NY2) 


TURB 

15 

c 


NWALL-2 WHEN THERE ARE TWO WALLS (I.E.# WALLS AT 

BOTH NYl AND 

NY2) . 

TURB 

16 

c 


NSTR-i FOR THE REGION IN WHICH THE LOWER SURFACE 

OF 

THE STRUT 

FORMS 

TUK8 

17 


17 
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C THE UPPER WALL OF THE FLuW DOMAIN (NY2-MS1). TUR3 

C NSTR-0 OTHERWISE TORB 

C TURB 

C USE AS MANY CALLS TD SUBROUTINE EDDY AS NECESSARY TO DEFINE ENTIRE KEGTURB 

C SEE examples IN THE USER'S MANUAL FOR FURTHER CLARIFICATION. TURB 

C TURB 

CALL E0DY(2,15,l.Ml»l,O) TURB 

CALL ED0Y(16»A6»MS>M1.2>0) TURB 

CALL EDDY(16,23>l>hSl>l>l) TURB 

CALL EDDY( 2A, A6, 1»MS1»2>1 ) TURB 

Call eody(a7#ni»i»mi>i/0) turb 

return turb 

END TURu 

SUBROUTINE XY( ETAMN»ETaMX>BETA) XY 

CtlMM0N/Fl/NXM»NXMl»NXH2»NXM3#NXMA>Nl>nl>Nll# Mll/NS 1/NS2#NS» MS# MSI XY 
l#NCWL#NCWLl#NCWL2#NCWLM#NCwLP»NSTRUT XY 

COMrtGN/FAl/K12#M13#N12#Ni3#NXM5»NXM6#MS2# MS 3 # MS A # M SPi > MS P2# MS P3 XY 
1 #NjM1#NS1P1# NS2M1»nS1M1#NS2P1»NNS#NS 1M2#NXM7 XY 

C0MM0N/F7/X(61#55)#Y(61»55»#XXI(61#55)#YXI (61# 55)# XY 

1 XETA(61#55)#YETA(61#55)#AJ(61#55) XY 

COMMON/FB/XS (55)# YS(55)#XSX1 (55)# YSXI(55)#XSETA(55 )# XY 

1 YSETA(55),AJS(55) XY 

DIMENSION ETA8(61)#ETAD(61) XY 

INTEGER ETAMX#ETAMN XY 

B-ALQG((BETA+1.)/(BETA-1.)) XY 

NETA-tTAMX-FTAMN XY 

DO 10 M-ETAMN# ETAMX XY 

A«(2.+(K-ETAMN)/NETA-1. )*B XY 

EXPA-EXP(A) XY 

ETAY«( (BETA+1. )*EXPA-3ETA+1. )/(l.+EXFA) XY 

ETA8 (M)-ETAY*NETA/2.+ETAMN XY 

10 ETA0(M)«2.*B*8ETA*£XPA/(1.+EXPA)**2 XY 

00 15 N-1»N1 XY 

DO 15 M-ETAMN#ETAMX XY 

X(M>N)»(X(ETAMX>N)*(ETAB(M)-ETAMN)+X(ETAMN#N)*(ETAMX-ETAB(M)) )/ XY 
INETA XY 

Y(i1»N)-( Y(ETAMX#N)+(ETAB(M)-ETAMN)+Y (ETAMN#N)7(ETAMX-ETAB(N) ) )/ XY 
INETA XY 

XETA (M#N)- (X(ETAMX#N)-X(ETAMN#N) )*ETAO(M)/NETA XY 

15 Y£TA(M#N)» (Y(ETA.MX#N)-Y(ETAMN#N))*ETAD(M)/NETA XY 

DO 20 N-2.N11 XY 

00 20 M-ETAMN»£TAMX XY 

XXI(M»N)«(X(M#N + 1)-X(H#N-1)) /2. XY 

20 YXI(M#N)-(Y(M,N+l)-Y(M#N-l))/2. XY 

DO 30 M«ETAMN#ETAMX XY 

XaI(M#1)«X(M#2)-X(M#1) XY 

YXI(M#1)«Y(K#2)-Y(M#1) XY 

XXI(M#N1)»X(M>N1)-X(M#N11) XY 

30 YXI (M»N1)«Y(M#N1)-Y(M#N11) XY 

RETURN XY 

END XY 

SUBROUTINE V IS C OUS ( S WE E P ) VISC 

BIT Bi Vise 

CCMMDN/F 1/NX M#NXM1#NXM2#NXM3#NXMA#M#M1# Nil# Ml 1# NS 1# NS2# NS# MS#MS 1 VI SC 

i#ncwl#ncwli#ncwl2»ncwlm#ncwlp#nstrut Vise 

C0MM0N/FA1/M12# M13#N12#N13#NXM5#NXM6#MS2#MS3#nSA#rtSPl#MSP2# MSP3 VISC 
1#NSM1»NS1P1#NS2M1#NS1 m1#NS2P1#NNS#NS1M2#NXM7 VISC 

COMMON/F2/RO(61,55)#U(61» 55 ) # V ( 61 # 55 ) # P ( 61 # 5 5) , T ( 61# 55) # VISC 

1 SH(61#55)#H(61#55)#R0S(55)#US(55)#VS(55)# VISC 

2 PS(55)#TS(55)>SHS(55)#HS(55) VISC 

C0MMQN/F3/VISL(51#55)# VIST(61#55)#VIS(61#55) #VISS( 55) VISC 

C£JMMaN/FA/SIGX(61#55)#SIGY(61,55)#TAUXY(61#55)#QX(61#55)# VISC 

1 UY(61#55)#SIGXS(55)#SIGYS(55)#TAUXYS(55)#C1XS(55)#(jYS(55) VISC 

COMMON/F5/AU(61»55#A)» AM(61,55#A) #AN(6i#55# A )# AU1(61#55#A) VISC 

CUMM0N/F6/AUS( 55# A)#AMS(55#A )#ANS (55#A)#AUS1 (55# A) VISC 

C0MM0N/F7/X(61»55)#Y(61#55)#XXI(61#55)#YXI(61# 55)# VISC 

1 XETA(61#55)#YETA(61#55),AJ(61#55) VISC 


TURB 

18 

TURB 

19 

TURB 

20 

iTURB 

21 

TURB 

22 

► TURB 

23 

TURB 


TURB 

25 

TURB 

26 

TURB 

27 

TURB 

28 

TURB 

29 

TURB 

30 

TURb 

31 

XT 

1 

xr 

2 

XY 

3 

XY 


XY 

5 

XY 

6 

XY 

7 

XY 

8 

XY 

9 

XY 

10 

XY 

11 

XY 

12 

XY 

13 

XY 

14 

XY 

15 

XY 

16 

XY 

17 

XY 

16 

XY 

19 

XY 

20 

XY 

21 

XY 

22 

XY 

23 

XY 

24 

XY 

25 

XY 

26 

XY 

27 

XY 

28 

XY 

29 

XY 

30 

XY 

31 

XY 

32 

XY 

33 

XY 

34 

XY 

35 

XY 

36 

XY 

37 

XY 

38 

vi:>cous 

1 

VISCOUS 

2 

VISCOUS 

3 

VISCOUS 

4 

VISCOUS 

5 

VISCOUS 

6 

VISCOUS 

7 

VISCOUS 

8 

VISCOUS 

9 

VISCOUS 

10 

VISCOUS 

11 

VISCOUS 

12 

VISCOUS 

13 

VISCOUS 

14 

VISCOUS 

15 

VISCOUS 

16 


18 
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C0f1f1LlN/FS/XS(55)»yS(55)»XSXl (55)>YSXI(i)5)»XSETA(55)> 

VISCOUS 

17 


i YjETA(55),AJS(55) 

VISCOUS 

18 


C0/1i1IjN/F9/INT( 55)»TE1(55)/T£2(55)/TE3(55)/TE4(55) 

VISCOUS 

19 


CUrtrt0N/F10/DT(61» 55),0TM 

VISCOUS 

20 


C0.1f1QN/Fll/ERR3(61>55),Bl (61,55) 

VISCOUS 

21 


C0Mf1QN/F12/LSY(1»MFUJH» 1. FI, LFO, ILT, UMAX, LW,L EXIT, CR IT, FDT, FOIL, 

VISCOUS 

22 


1 TIrtE,CCP,CCT,L,LK,TIi’1EL,CRITAVG 

VISCOUS 

23 


CUi1^1UN/F13/GAMA,RGAS,CP,PR, PRT,CV1S, Fn, PF, TF,UF, ROF, SHF, HF 

VISCOUS 

24 


CUfHaN/FlA/DAU(51,55, A),9DAU(55, A ),6DAUS (55, A),T1 (55,5), 

VISCOUS 

25 


lCCIMPY,Gt11,GM2,6M3,GhA,G^15 

VISCOUS 

26 



VISCOUS 

27 


INITIAL GUESS 

VISCOUS 

28 



VISCOUS 

29 


U(l, 1;NXM)-UF 

VISCOUS 

30 


V(i,l;NXH)-0. 

VISCOUS 

31 


T(l,ljNXM)-TF 

VISCOUS 

32 


sh(i,i;nxm)-shf 

VISCOUS 

33 


H(1,1)NXM)«HF 

VISCOUS 

34 


P(l,liNXM)»PF 

VISCOUS 

35 


RCi( l,llNXM)«P(l,l;NXrt)/RGAS/T(l,lINXh) 

VISCOUS 

36 


IFdtT.EQ.DGO TO 2 

VISCOUS 

37 


REA0(LFI,561) L,TIME 

VISCOUS 

38 


DO 3 N-1,N1 

VISCOUS 

39 


R£A0(LFI, 560) (X(M,N),Y(M,N),U(M,N), V(t1»N),P (h,N),T(rt,N),RO(M,N) , 

VISCOUS 

40 


lSH(rt,N),H(M,N),VIST(M,N),M»l,Ml) 

VISCOUS 

41 

3 

CONTINUE 

VISCOUS 

42 


LHAX-LMAX+L 

VISCOUS 

43 


iF(NSTRUT.EQ.O)GO TO 2 

VISCOUS 

44 


READ(LFI,560)(XS(N),YS(N),US(N),VS(N),PS(N),TS(N),RJS(N),SHS(N), 

VISCOUS 

45 


1HS(N),VISS(N),N«1,NS) 

VISCOUS 

46 

2 

CONTINUE 

VISCOUS 

47 


L-L + 1 

VISCOUS 

40 


CALL BOUND 

VISCOUS 

49 

1 

CONTINUE 

VISCOUS 

50 


LSTEP-1 

VISCOUS 

51 


£RR0(1,1INXN )«RO(l,iINXM) 

VISCOUS 

52 


QX(2,2JNXM4)«GAMA*P(2,2}NXHA)/R0(2,2JNXMA) 

VISCOUS 

53 


uX(2,2}NXMA)-VSQRT(ax(2,2;NXMA);QX(2,2jNXMA) ) 

VISCOUS 

54 


DT(2,2)NXMA)«FDT*YETA(2,2;NXMA) /( VABS(V(2,2jNX«4); 

VISCOUS 

55 


1QY(2,2)NXNA) )+0X(2,2}NXMA)) 

VISCOUS 

56 


DTM«08SMIN(DT(2,2;NXM4)) 

VISCOUS 

57 


DTNL«FDTL*DTM 

VISCOUS 

58 


B1(2,2)NXM4)«DT(2,2JNXM4) .GT.OTML 

VISCOUS 

59 


DT(2,2)NXN4)»08YCTRL( JTML,6I(2,2)NXM4))DT(2, 2JNXM4)) 

VISCOUS 

60 


time-time+dtm 

VISCOUS 

61 


call dviscp 

VISCOUS 

62 


CALL VECl 

VISCOUS 

63 


CALL DANPP 

VISCOUS 

64 



VISCOUS 

65 


PREDICTOR STEP WITH FORWARD DIFFERENCING. 

VISCOUS 

66 



VISCOUS 

67 


1F(FDT.GT.1. )60 TO 10 

VISCOUS 

68 


DO 100 I»l,4 

VISCOUS 

69 


Aol(2,2,I;NXM4)-AU(2,2,I;NXH4>-0T(2,2;NXf14)* ((AM(2,3,I;NX«4)-AM(2, VISCOUS 

70 


12,IJNXM4 ) ) + ( AN( 3,2,1 1 NXM4)-AN( 2,2,1 J NXM4) ) ) 

VISCOUS 

71 


IKNSTRUT.E0,0)G0 TO 100 

VISCOUS 

72 


DO 105 N»1,NS 

VISCOUS 

73 


NN-N+NSl-1 

VISCOUS 

74 


AUl(hSl,NN»I )-AU(MSl,NN, D-DT (MSI, NN)* (AN (MSI, NN + l,I)-An(nSl,NN,I) VISCOUS 

75 


1+AN5(N#I)-AS (MS1#NN,I) ) 

VISCOUS 

76 

105 

CONTINUE 

VISCOUS 

77 

100 

CONTINUE 

VISCOUS 

78 


CALL VEC2 

VISCOUS 

79 


CALL aauNo 

VISCOUS 

80 


GO TO 20 

VISCOUS 

81 

10 

CONTINUE 

VISCOUS 

82 


DO 15 I-IM 

VISCOUS 

63 


DAU(2^2>IiNXM4)--DT(2,2;NXrt4)+(Ah(2# 3/ 1 J NX M4 )-Ah ( 2/ 2# I j NXN4 ) + 

VISCOUS 

84 


19 
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16 

15 


17 

C 


20 


IbO 


C 

C 

C 


210 


30 


35 


36 


C 


‘rO 


c 

C 

C 


1AN(3>2> I JNXh4)-AN( 2»2» IjMXM'. ) ) /AJ (2> 21NXM4) 

IMNSTRUTtFO.O)GO TG 15 
OCj 16 N»1,NS 
NN"N+NS1“1 

0AJ(I1S1» I )--DT(M51/NN) ♦(AM(MSl»NN + l#i )-AM (MS1/NN>I)+ANS(N#I)- 

lAN(rtSl»NN> I ) )/AJ (MS1»NN) 

CONTINUE 

continue 

lAOO«l 

CALL IhPY(lADO) 

00 17 I-l/'f 

AU1(2>2»I JNXMA)«AU(2»2>1;NXM4)+OAO(2>2»I{NXM4)*AJ(2>2jNXN4) 
SUBROUTINE BOUND IS CALLED IN VECIM2. 

CALL VECIN2(IA00) 
continue 

1MLK.FQ.2)RETURN 
DO 160 I»1,A 

AU1(2»2, I}NXMA)»AU(2»2>I}NXf1A) +AU1(2>2»I}NXMA) 

continue 

CALL DVISCC 

CALL VECl 

Call damp: 

corrector step WITH BACKwORD DIFFERENCING 

IF (FOT .GT. 1. )Gl3 TO 30 
DO 210 I-1«A 

A01(2^2,I}NXMA)-D.5+(AU1(2#2,I|NXM'»)-0T(2»2J NXM4 )*(( AM ( 2 #2» I iNXrt f 
1-Art( 2>1»I JNXMA ) ) +( AN<2» 2» I J HXM4 )-AN ( 1 #2# I JNX ML) ) ) ) 

CONTINUE 
CALL VEC2 
CALL BOUND 
GO TO AO 

continue 

00 35 I-1,A 

0AU(2,2,I )NXMA)--DT(2,2JNXMA)*<AM(2,2>I;NXMA )-AM ( 2 , 1> I } NXHA ) + 
1AN(2/2>I;NXMA)-AN(1,2»IjNXMA) ) /AJ (2/ 2JNXMA) 

CONTINUE 

IA00«0 

CALL IMPY(IaDD) 

DO 36 1-1,4 
AUl ( 2, 2, IINXM4 ) - . 

IAOD-1 

SUBROUTINE BOUND 
CAlL VECIM2(IADD) 

CONTINUE 

1F(LR.EQ.2 IRETURN 
L-L + 1 
LSIEP-2 

JX{2,2jNXM4)-GAMA*P(2»2}NXM4)/R0(2»2iNXM4) 
DX(2,2;NXM4)-VSQRT(QX(2,2INXM4);QX(2,2;NXM4) ) 
OT(2,2}NXM4)-F0r*YETA(2,2jNXM4) / ( VABS ( V (2, 2 ;NXM4 ) j 

1UY(2,2JNXM4) ) + 3X (2,2;NXM4) ) 

DTM-08SHIN(0T{2,2}NXM4) ) 

OTML-FDTL+DTM 

t)l(2>2|NXM4)-DT(2,2jNXM4) .GT.DTML 

0T(2,2}NXM4) -38VCTRL(DTMl,o1(2>2;NXM4);0T(2, 2jNXM4)) 

TIME-TIME+OTM 
CALL DVISCC 
call VECl 
CALL DAMPC 

PREDICTOR STEP WITH BACKWARD DIFFERENCING. 

IF(FDT.GT.1. )GJ TO 50 
DO 110 1-1,4 

AU1(2,2,I JNXM4 )-AU( 2,2,1; NXM4)-DT( 2,2 ;NXM4)* ((AN (2,2,1} NXM4)-; 
11,I}NXM4) > + (AN(2,2, II nXM4)-aN( 1,2,1 ;NXh4)) ) 


5*(AJ1(2,2, I,'NXN4)+DAU(2,2, I;NXM4 ) *A J ( 2, 2; NXi14 ) ) 
IS Called in vecim2. 


VIISCUUS 

B5 

VISCOUS 

66 

VISCOUS 

87 

VISCOUS 

88 

VISCOUS 

89 

VISCOUS 

90 

VISCOUS 

91 

VISCOUS 

92 

VISCOUS 

93 

VISCOUS 

94 

VISCOUS 

95 

VISCOUS 

96 

VISCOUS 

97 

VISCOUS 

98 

VISCOUS 

99 

VISCOUS 

100 

VISCOUS 

101 

VISCOUS 

102 

VISCOUS 

103 

VISCOUS 

104 

VISCOUS 

105 

VISCOUS 

106 

VISCOUS 

107 

VISCOUS 

108 

VISCOUS 

109 

VISCOUS 

110 

VISCOUS 

111 

IVISCOUS 

112 

VISCOUS 

113 

VISCOUS 

114 

VISCOUS 

115 

VISCOUS 

116 

VISCOUS 

117 

VISCOUS 

lie 

VISCOUS 

119 

VISCOUS 

120 

VISCOUS 

121 

VISCOUS 

122 

VISCOUS 

123 

VISCOUS 

124 

VISCOUS 

125 

VISCOUS 

126 

VISCOUS 

127 

VISCOUS 

12b 

VISCOUS 

129 

VISCOUS 

130 

VISCOUS 

131 

VISCOUS 

132 

VISCOUS 

133 

VISCOUS 

134 

VISCOUS 

135 

VISCOUS 

136 

VISCOUS 

137 

VISCOUS 

138 

VISCOUS 

139 

viscous 

140 

VISCOUS 

141 

VISCOUS 

142 

VISCOUS 

143 

VISCOUS 

144 

VISCOUS 

145 

VISCOUS 

146 

VISCOUS 

147 

VISCOUS 

146 

VISCOUS 

149 

VISCOUS 

150 

VISCOUS 

151 

VISCOUS 

152 


20 



o o o 


APPENDIX A 


110 

Continue 

VISCOUS 

153 


CALL VEC2 

VISCOUS 

154 


CALL 30UND 

VISCOUS 

155 


Gu TO 60 

VISCOUS 

156 

bO 

CONTINUE 

VISCOUS 

157 


00 55 

VISCOUS 

158 


0aU( 2/2> I;NXMA)«-DT(2»2;NXNA)*(AM (2/ 2>1;NXMA )-Ah(2/l» i;nxma » + 

VISCOUS 

159 


1AN(2/2,I JNXMA)-AN(1^2, IjNXKA) )./AJ (2>2JNXMA) 

VISCOUS 

160 

bb 

CONTINUE 

VISCGU:) 

161 


lAOO-0 

VISCOUS 

162 


CALL IMPr(IADD) 

VISCOUS 

163 


JO 56 I»1M 

VISCOUS 

164 


AU1(2»2, I JNXMA )«AJ(2>2^ IjNXrtA) +DAU(2 »2»I}NXMA)*AJ ( 2,2}NXf1A) 

VISCOUS 

165 

56 

continue 

VISCOUS 

166 


SU3K0UTINE BOUND IS CALLED IN VECIM2. 

VISCOUS 

167 


Call ve:in2(iaod) 

VISCOUS 

168 

60 

CONTINUE 

VISCOUS 

169 


IE(LK.EQ.2)RETJRN 

VISCOUS 

170 


DO 165 I«1»A 

VISCOUS 

171 


AUl(2#2,I}NXNA»»AU(2>2iI{NXNA)+AUl(2»2/IjNXMA) 

VISCOUS 

172 


IF(NSTRUT.EQ.O)GQ TO 165 

VISCOUS 

173 


00 150 N«1,NS 

VISCOUS 

174 


NN«N+NS1-1 

VISCOUS 

175 

150 

AUS1(N>I)-AU1{NS1»NN»I ) 

VISCOUS 

176 

165 

CONTINUE 

VISCOUS 

177 


Call oviscf 

VISCOUS 

178 


CALL VECl 

VISCOUS 

179 


call dampp 

VISCOUS 

180 



VISCOUS 

181 


CORRECTOR STEP WITH FORWARD DIFFERENCING 

VISCOUS 

182 



VISCOUS 

183 


IF(FOT.GT.l. )G0 TO 70 

VISCOUS 

184 


DO 215 I-1,A 

VISCOUS 

185 


AU1(2»2,I)NXMA)-0.5»(AU1(2>2>I}NXH4)-OT(2>2;NXM4)*( (AH(2>3# I}NXN4)VISC0US 

186 


1-AM( 2>2j IINXMA ) )+ (AN(3>2^ I)NXNA)-AN( 2>2> I ;NXMA) ) ) ) 

VISCOUS 

187 


IF(NSTRUT.EQ.O)GQ TO 215 

VISCOUS 

188 


DU 120 N«1>NS 

VISCOUS 

189 


NN-N+NSl-1 

VISCOUS 

190 


AU1(HS1>NN,I )«3.5*(AUS1(N>I)-0T(MS1> NN) *( AM< hS 1» NN + 1, I )-AM( MS 1> NN, VISCOUS 

191 


1I) + ANS(N»I)-AN(MS1»NN> I))) 

VISCOUS 

192 

120 

CONTINUE 

VISCOUS 

193 

215 

CONTINUE 

VISCOUS 

194 


CALL VEC2 

VISCOU;> 

195 


CALL BOUND 

VISCOUS 

196 


GO TO 80 

VISCOUS 

197 

70 

CONTINUE 

VISCOUS 

198 


DO 75 I-liA 

VISCOUS 

199 


0AU(2j 2>I iNXMA )--DT (2» 2JNXMA) ♦( AM (2» 3»l;NXMA )-AM(2,2, DNXMA ) + 

VISCOUS 

200 


1AN( 3/2» IINXMA )-AN(2#2» I J NXM 4 ) ) / A J ( 2» 2iNXM4 ) 

VISCOUS 

201 


IF(NSTRUT.EQ.O)GO TO 75 

VISCOUS 

202 


00 76 N-1,NS 

VISCOUS 

203 


NN-N+NSl-1 

VISCOUS 

204 


0AJ(MS1»NN»I )--DT(MSl>NN)*(AM(hSl/NN + l>I )-AM IMS1>NN,I) + 

VISCOUS 

205 


1ANS(N/I)-AN(MS1>NN,I) )/AJ(MS1/NN) 

VISCOUS 

206 

76 

CONTINUE 

VISCOUS 

207 

75 

CONTINUE 

VISCOUS 

208 


1 A0D«1 

VISCOUS 

209 


CALL IMPY(IADD) 

VISCOUS 

210 


DO 77 I-1/4 

VISCOUS 

211 

77 

AU1(2/2»I;NXM4)-.5*(AJ1(2#2» IJNXM4)+DAU(2»2> 1 ; NXM4 ) * A J ( 2 a2 I NXM4 ) ) 

VISCOUS 

212 


lAOD-0 

VISCOUS 

213 


SUBROUTINE BOUND IS CALLED iN VECIM2. 

VISCOUS 

214 


Call VECIM2(IA0D) 

VISCOUS 

215 

80 

CONTINUE 

VISCOUS 

216 


IF (LK.EQ.2 1RETURN 

VISCOUS 

217 


IF(NFLQW.NE. 2)G0 TO 216 

VISCOUS 

218 


LL-L/20 

VISCOUS 

219 


LL-LL+20 

VISCOUS 

220 
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1F(LL.SE.L)G0 TO 216 

\^ISCOUS 

221 


CAUL TU?3 

VISCOUS 

222 

216 

CONTINUE 

VISCOUS 

223 


ERROdjl }NXr,)«(RJ( 1JNXM)-ERRD(1>1; NXM) )/ERR0(l#UNXf1) / FOT/LSTEP 

VISCOUS 

224 


EKROdf 1 JNXM )«VAtJS IERR0(1»UNXM>} ERROdi ijNXM) ) 

VISCOUS 

225 


ERR.'lX-08S‘1AX(ERRa(l»ljNXM)) 

VISCOUS 

226 


1F(£RRMX.LE.CRIT)G0 TU 217 

VISCOUS 

227 


IFaiME.GT.TlMEDGj TO 217 

VISCOUS 

228 


LL-L/50 

VISCOUS 

229 


LL-LL*50 

VISCOUS 

230 


IF(LL.SE.L)GD TO 220 

VISCOUS 

231 


SIGX(l»l}NXf1)-ERR0(l»ljNXM)*ERR0(l/l JNXM) 

VISCOUS 

232 


EPRAVG»Q8bSUH(SIGX(l»UNXM) ) 

VISCOUS 

233 


EKKAVG»SQRT(ERRAVG)/NXM 

VISCOUS 

234 


1F(EKRAVG. LE .CRITAVG) GO TO 217 

VISCOUS 

235 


4PITE(6»510)L>ERR,iX»ERRAV6 

VISCOUS 

236 

610 

FOR.iAT(/dOX, 'L-',15,10X, • MAX . ERROR- ' >E 12 . 10X> 

VISCOUS 

237 


I'MEAN SORT ERROR- SE12.5) 

VISCOUS 

238 


LL-L/LW 

VISCOUS 

239 


LL-LL*LW 

VISCOUS 

240 


IFdL.NE.DGO TO 220 

VISCOUS 

241 

217 

continue 

VISCOUS 

242 


Call print 

VISCOUS 

243 


CALL SPILL(SWEEP) 

VISCOUS 

244 


IF (ERRAVG.LE .CRITAVG) i<P I T E 1 6» 600 ) 

VISCOUS 

245 


IK EkkAVG.LE. CRITAVG) GO TO 218 

VISCOUS 

246 


IF (TlME.GT.TIMEL ) hRITE(6/610) 

VISCOUS 

247 


Ir (TIME.GT.TIMEL ) GO TO 218 

VISCOUS 

248 


IF(ERRMX.LE.CRIT) WR1TE(6»620) 

VISCOUS 

249 


1F(ERRMX.LE.CRIT)GQ TO 218 

VISCOUS 

250 

220 

continue 

VISCOUS 

251 


L-L + 1 

VISCOUS 

252 


IF (L.lE.LMAX )G0 TO 1 

VISCOUS 

253 


L-L-1 

VISCOUS 

254 


WRITE (6,630 ) 

VISCOUS 

255 

218 

CONTINUE 

VISCOUS 

256 


REWIND LED 

VISCOUS 

257 


WkIT£(LF0,561) L,TIME 

VISCOUS 

258 


DO 300 N-1,N1 

VISCOUS 

259 


WRITE(LF 0,560) (X(M,N),Y(li,N),U(M,N),V(M,N),P(M,N),T(M,N),R0(.1,N), 

VISCOUS 

260 


ISH (M,N),H( M, N) , VIST(M,N),M-l,iil ) 

VISCOUS 

261 

300 

CONTINUE 

VISCOUS 

262 


IF(NSTRUT.EO.O)GJ TO 310 

VISCOUS 

263 


WkITE(LFQ,560)(XS(N)»YS(N)»uS(N),VS(N),PS(N),TS(N),R0S(N),SHS(N), 

VISCOUS 

264 


1HS(N),VISS(N),N-1,NS) 

VISCOUS 

265 

310 

CONTINUE 

VISCOUS 

266 

560 

FORMAT (IX, 2F9.6,2F12.5,F1A.5, F13.5, F11.5,F1A.5,F15.6,F10.6) 

VISCOUS 

267 

561 

F0RMaT(1X,I6,2X,£15.7) 

VISCOUS 

266 

600 

FOR:1AT(//,20X, 'PROGRAM TERMINATED ON MEAN SORT ERROR CRITERION') 

VISCOUS 

269 

610 

F0RMAT(//,20X, 'PROGRAM TERMINATED ON PHYSICAL TIME ', 

VISCOUS 

270 


1 'CONVERGENCE CRITERION*) 

VISCOUS 

271 

620 

F0RMAT(//,2C'X, 'PROGRAM TERMINATED ON LOCAL ERROR CRITERION') 

VISCOUS 

272 

630 

FORMaK //,20X, 'PROGRAM TERMINATED ON UMAX*) 

VISCOUS 

273 


RETURN 

VISCOUS 

274 


END 

VISCOUS 

275 


SU8RUUTINE INVICIO(SWEEF) 

INVICID 

1 


BIT 81 

INVICIO 

2 


C0MMUN/F1/NXM,NXM1,NXM2,NXM3,NXMA,N1,M1,N11, M11,NS1,NS2,NS,MS,MS1 

INVICID 

3 


1,NCWL,NCWL1,NCWL2,NCWlM,NCwLP,NSTRUT 

INVICID 

4 


CUMMON/F Al/M12,M13,N12,N13,NXH5,NXM6,MS2,MS3#MSA,MSPi,MSP2, MSP3 

INVICID 

5 


1»NSM1,NS1P1, NS2Ml,NSli11,NS2Pl,NNS,NSlM2,NXM7 

INVICIO 

6 


CUMMUN/F2/R0(61, 55 ) , U ( 61, 55 ) , V (61, 55 ), P ( 61, 5 5 ) , T ( 6 1, 55 ) , 

INVICID 

7 


1 SH(6l,55),H(61,55),P0S(55),US(55),VS(55), 

INVICID 

8 


2 PS(55),TS(55),SHS(55),H3(55) 

INVICIO 

9 


C0MMQN/F3/VI SL ( 61, 55 ) , VI S T ( 61, 55 ) , VI S ( 61, 55) , V I SS ( 55 ) 

INVICIO 

10 


C0MM0N/FA/SIGX(51,55),SIGY(61,55),TAUXY(61,55),QX(61,55), 

INVICID 

11 


1 QY(61, 55),SIGXS(55),SIGYS(55),TAUXYS(55),(lXS(55),QYS(55) 

INVICIO 

12 


C0MM0N/F5/AU(61, 55, A ) , AM( 61, 55, A ) , AN ( 61, 55, A ), AUl ( 61, 55, A ) 

INVICID 

13 
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CuHrtLJN/F6/AUS(5 5,4)>A«S(5 5M)>ANS(55/4)#AUSl (55>4) 

INVICID 

14 

CQMt10N/F7/X(61»55)>Y(61#55)»XXl(61»5 5)»YXI (61,55), 

INVICID 

15 

1 XEfA(61, 55) ,YETA(61,55),AJ( 61,55) 

INVICID 

16 

C0(1M0N/F8/XS(55),YS(55),XiXI(55), YSXI <55),XSETA(55), 

INVICID 

17 

1 YSETA(55),AJS(55) 

INVICID 

16 

CCrtilOH/FP/INK 55),T£1(55),TE2(55),TE3(55),TEA(55) 

INVICID 

19 

Cul1rt0N/F10/DT(61,55),0T(1 

INVICID 

20 

COrtrlON/F 11 /ERRO( 61,55 ),B1 (61,55) 

INVICID 

21 

CDMMt3N/F12/LSYM,NFLOW,LFl,LFO,lLT,LMAX,LW,LEXIT,CRIT,FDT,HOTL, 

INVICID 

22 

1 TirtE,CCP,CCT,L,LK,TlrlEL,CRITAVG 

INVICID 

23 

COM,10N/F13/GAMA,RGAS,CP,PR,PRT,CVIS, F H, PF , TF , UF, RQF, SHF , HF 

INVICID 

24 


INVICID 

25 

INITIAL GUESS 

INVICID 

26 


INVICID 

27 

V(l,l)NX>1)-YXI(l,lJNXrl)/XXI(l,l;NXM)*UF 

INVICID 

28 

U(l,l;NXM)-UF*UF-Y(l,ljNXM)+V(l,l;NXM) 

INVICID 

29 

U(1,1>NXM)-VSQRT(J(1, IJNxh) iU{l,llNXh)) 

INVICID 

30 

1 (l,ljNX'1)»TF 

INVICID 

31 

p(i,i;nxh)«pf 

INVICID 

32 

ko(i,i;nxm)-rqf 

INVICID 

33 

SH(l,l)NX'1)»SHF 

INVICID 

34 

HU,liNXM)«HF 

INVICID 

35 

IF(NSTRUT.EQ.O) GU to 6 

INVICID 

36 

VS(liNS)«YSXI(l;NS)/XSXI(l}NS)*UF 

INVICID 

37 

USdlNS )»UF*UF-VS( IJNS )*VS(l;NS) 

INVICID 

38 

USdiNS )»VSQRT(US (1;NS)TUS(1 ;nS )) 

INVICID 

39 

PSdJNS)«PF 

INVICID 

40 

TSd}NS)-TF 

INVICID 

41 

SHSdjNS )»SHF 

INVICID 

42 

R(JS(UNS)«R0F 

INVICID 

43 

hs(i;ns)«hf 

INVICID 

44 

CONTINUE 

INVICID 

45 

1F(ILT.EQ.1)G3 to 2 

INVICID 

46 

REAOdFI, 561 ) L,TIME 

INVICID 

47 

00 3 N«1,N1 

INVICID 

46 

R£AD(LFI,560)(X(M,N),Y((1,N),U(M,N),V(n,N),P(h,N),T(rt,N),R0(M,N), 

INVICID 

49 

lSH(t1,N),H(M,N), YIST(M,N),M«1,M1) 

INVICID 

50 

CONTINUE 

la^icio 

51 

LNAX-LMAX+L 

INVICID 

52 

IF{NSTRUT.£0.0)GO to 2 

INVICID 

53 

READ(LFI,56C)(XS(N),YS(N),US{N),vS(N), PS ( N ) , TS ( N ) , RJS ( N ) , S HS ( N ) , 

INVICID 

54 

IHS (N),VISS (N), N-1,NS) 

INVICID 

55 

CONTINUE 

INVICID 

56 

L«L + 1 

INVICID 

57 

CALL IBDUMD 

INVICID 

58 

CONTINUE 

INVICID 

59 

LSTEP-1 

INVICID 

60 

EkRO(l,ljNXP, )-RO(l,UNXh) 

INVICID 

61 

JX(2,2)NXttA)»GA1A*P(2,2;NXhA)/R0(2,2iNXMA ) 

INVICID 

62 

OX(2,2}NXhA) -VSQRT(QX(2,2jNXM^);wX(2,2iNXNA) ) 

INVICID 

63 

0T(2,2JNX^14)•FDT*YtTA(2,^|NXM6) /( VABS( V (2,2|NXMA ) j 

INVICID 

64 

1QY(2,2JNXMA ) ) ♦OX ( 2, 2 J NXMA ) ) 

INVICID 

65 

0TI1»g8SMIN(DT(2,2jNXM4) ) 

INVICID 

66 

01/1L»FDTL*DTM 

INVICID 

67 

B1(2,2}NXM4)-DT(2,2|NXM4) ♦GT.DTML 

INVICID 

66 

0T(2,2)NX«4)«a8\/CTRL(0TML,61(2,2{NXM4)jDT(2, 2)NXM4)) 

INVICID 

69 

T1M£«TIME+DTM 

INVICID 

70 

CALL IVECl 

INVICID 

71 

CALL OAfIPP 

INVICID 

72 


INVICID 

73 

PREDICTOR STEP WITH FORWARD DIFFERENCING. 

INVICID 

74 


INVICID 

75 

00 100 1*1,4 

INVICID 

76 

AUK 2/ 2/ I ;NXM4)»AJ(2/2/I;nXM 4)-DT (2/ 2 jNXM4 )★ ( ( Art ( 2/ 3# I ; NXI14 )-A« ( 2# INVIC ID 

77 

12,IINXM4) ) + (AN(3,2,dNXM4)-AN(2,2,IjNXtl4) ) ) 

INVICID 

78 

IF(NSTRUT.EQ.O)GO TO 100 

INVICID 

79 

DU 105 N-1,NS 

INVICID 

80 

NN«N+NS1-1 

INVICID 

81 
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105 

100 


160 


ZiO 


110 


150 

165 


120 

215 


AUi(r,Sl» I )»AU(iSl>NN/l )-DT(MSl>NN )*JAM (MS 1» NN + 1»D- 

-AM(MS1#NN>I)I.SVICID 

62 

1 + ANS(N#I)-AN(MS1»NN»I) ) 

INVICID 

83 

CONFiNUE 

INVICID 

84 

continue 

INVICIO 

85 

call VEC2 

INVICID 

86 

CALL IBOUND 

INVICIO 

87 

IF(LK.EQ.2)kETURN 

INVICID 

86 

DO 160 I»1>A 

INVICID 

69 

A01(2>2> UNXIiA)«A'J(2»2»IiNXHA)+AUl(2>2^I»NXMA) 

INVICIO 

90 

CONTINUE 

INVICID 

91 

CALL IVECl 

INVICID 

92 

CALL OAMPC 

INVICID 

93 


INVICID 

94 

CUKRECTOR STEP WITH BACKWORO DIFFERENCING 

INVICIO 

95 


INVICID 

96 

00 210 I-l/A 

INVICID 

97 

AUK 2>2> I ;NXMA )-0.5*(AUl( 2>2> l;MXMA)-UT(2>2; NXMA)* ( (AM( 2>2» UNXMA) INVIC ID 

98 

l-AiK 2»I» IJNXMA ) ) + ( AN(2/2> iiNXMA)-AN( 1>2K )NX MA ) ) ) ) 

INVICIO 

99 

CONTINUE 

INVICID 

100 

CALL VEC2 

INVICID 

101 

CALL IBOUND 

INVICIO 

102 

IF(LK.Ea.2)RETJRN 

INVICID 

103 

L-L + 1 

INVICID 

104 

LSFEP-2 

INVICID 

105 

WX(2>2JNXMD -GANA*P(2»2;NXM'»)/RU(2/2iNXMA) 

INVICID 

106 

QX(2>2;NXMA)-VSQRT(aX(2>2iNXMA)iaX(2>2;NXMA) ) 

INVICID 

107 

DT(2^ 2;NXiiA )-FDT*YETA(2»2iNXMA) / ( VAB S ( V (2» 2} NXMA ) i INVICIO 

108 

1WY(2>2JNXHA) ) + OX < 2» 2i NXMA ) ) 

INVICIO 

109 

DTM-a8SMIN(DT(2»2jNXMA) ) 

INVICID 

110 

DTML-FDTL+DTM 

INVICIO 

111 

B1(2>2}NXMA)-DT(2»2JNXMA) .GT.OTML 

INVICIO 

112 

DT (2i2jNXMA) »Q8YCTRL(0TML»B1( 2>2}NXM A) IDT (2» 2JNXMA )) 

INVICID 

113 

TIME-TIME+OTM 

INVICID 

114 

CALL IVECl 

INVICID 

115 

call dampc 

INVICID 

116 


INVICIO 

117 

PREDICTOR STEP WITH BACKWARD DIFFERENCING. 

INVICIO 

118 


INVICIO 

119 

DO 110 I«1,A 

INVICID 

120 

AUl(2>2»UNXMA)-AU(2/2/liNXhA)-01 (2>2jNXMA)* ( ( AM ( 2 / 2, 1 J NXMA )-AM( 2» INVICI D 

121 

11> IINXMA) )+( AN{2»2jIJNXMA)-AN(1»2/1 JNXMA) ) ) 

INVICIO 

122 

CONFiNUE 

INVICIO 

123 

CALL VEC2 

INVICID 

124 

Call IBOUND 

INVICID 

125 

1F(LK.E3.2)PETJRN 

INVICID 

126 

DO lt>5 I-1»A 

INVICID 

127 

AUl(2>2>nNXMA)«AU(2»2K;NXMA) +AU1(2>2>I;NXHA) 

INVICIO 

128 

IE (NSTRUT.EQ.O)GO TO 165 

INVICID 

129 

DO 150 N«1>nS 

INVICIO 

130 

NN-N+NSl-1 

INVICID 

131 

AUS1(N>I )-AUl(MSl»NN,I ) 

INVICID 

132 

CONTINUE 

INVICID 

133 

call IVECl 

INVICID 

134 

CALL DAMPP 

INVICID 

135 


INVICID 

136 

CORRECTOR STEP WITH FORWARD DIFFERENCING 

INVICID 

137 


INVICID 

138 

DO 215 I-1,A 

INVICID 

139 

AUl(2>2>l;NXMA)«O.5*(AUl(2»2»nNXMA)-OT(2#2lNXMA)*((AM(2/3»IjNXMA)INVlCI0 

140 

1-AM(2>2»I JNXMA))+(AN(3>2> I; NXMA )-AN ( 2 > 2> 1 1 NX MA ) ) ) ) 

INVICID 

141 

1F(NSTRUT.EQ.0)G0 TD 215 

INVICIO 

142 

DO 120 N«1^NS 

INVICID 

143 

NN-N+NSl-1 

INVICIO 

144 

AU1(MS1»NNK )»D.5*(AUS1(N,I)-DT(MS1>NNK( AM( NS1»NN+1» 

I)-AM(MSi>NN>INVICID 

145 

ID+ANS (N» I )-AN(M5l»NN> liJ) 

INVICIO 

146 

continue 

INVICID 

147 

CONTINUE 

INVICID 

148 

CALL VEC2 

INVICID 

149 
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CALL IBDJND 

INVICID 

150 


IF«LK.£0.2)RETURN 

INVXCIO 

151 


ER«J(1>1 JNXM )«(R0(l»l|NXi)-ER«0(l>l jNXh) )/EftR0(l#l jNXH)/F0T/LSTEP 

INVICID 

152 


EkRa(l»l}NXM)«VABS(ERRO(l>ljNXM)JERRD(l>l;NXH) ) 

INVICID 

153 


ERRMX-08SMAX (ERRQ( liNXrt) ) 

INVICID 

154 


If (ERRMX.LE.CRIT)GO to 217 

INVICID 

155 


IF (TiME.GT.inEDGO TJ 217 

INVICID 

156 


LL-L/50 

INVICID 

157 


LL-LL*50 

INVICID 

158 


IF(LL.NE.L )G0 TO 220 

INVICID 

159 


ilGX{l/l}«fX« )»ERRa(l»UNXrt)*ERRO( I^IINXM) 

INVICID 

160 


tkRAVG«Q8SSUH(SIGX(l»l)NXf1) ) 

INVICID 

161 


EkRAVG«iQRT(ERRA\/G)/NXM 

INVICID 

162 


IF(EKKAVG.LE.CRITAVG) go to 217 

INVICID 

163 


WRITE{6>510)L>ERRn/EKRAtfG 

INVICID 

164 

5^10 

FORMATC /»10X/ 'L- '»I5»10X» 'MAX. ERROR- '/El 2.5 >10X> 

INVICID 

165 


I'F.EAN SORT ERROR- SE12.5) 

INVICID 

166 


LL-L/LR 

INVICID 

167 


lL-LL+LW 

INVICID 

166 


IFUL.NE.DGO to 220 

INVICID 

169 

217 

CONTINUE 

INVICID 

170 


Call PRINT 

INVICID 

171 


call SPILL(SWEEP) 

INVICID 

172 


IFCERkAVG.LE.CRITAVG) «RITE( 6/600) 

INVICID 

173 


IFCERkAVG.LE .CRITAVG) GO TO 218 

INVICID 

174 


IFITIME.GT. TIMED WRirE( 6 / 610 ) 

INVICID 

175 


IFCTIflE.GT.TIMEL ) GO TO 218 

INVICID 

176 


IF(ERRMX.LE.CRIT) rR1TE(6/620) 

INVICID 

177 


1M£RRMX.LE.CRIT)G0 TO 216 

INVICID 

170 

22C 

CONTINUE 

INVICID 

179 


L-L + 1 

INVICID 

180 


IF(L.LE.LMAX )G0 TO 1 

INVICID 

181 


L-L-1 

INVICID 

182 


WRIT6(6»630) 

INVICID 

183 

216 

CONTINUE 

INVICID 

184 


RERiND LFD 

INVICID 

185 


rRITE(LF0/561) L/TIME 

INVICID 

186 


00 300 N-l/Nl 

INVICID 

187 


RkIT£(LF0/560) (X (M/N), T(f1/N)/U(M/N) / V(M/N)/P(M/N), T(M»N)/R0(M»N)/ 

INVICID 

188 


1SH(N,N) ,H(M> N), VIST(M>N)»f1-l>Ml) 

INVICID 

189 

300 

continue 

INVICID 

190 


IF(NSTRUT.E0.0)G0 to 310 

INVICID 

191 


RkIT6(LFO,560)(XS(N)»YS(N)»uS(N)/VS(N)/PS(N)/TS(N)/ROS(N)/SHS(N)/ 

INVICID 

192 


1HS(N)/VISS(N),N»1,NS) 

INVICID 

193 

310 

continue 

INVICID 

194 

ito 

FORMAT (IX/ 2F9.6/2F12.5>F1A.5>F13.5/F11.5/F14.5/F15.6,F10.6) 

INVICID 

195 

561 

F0RNAT(1X>I6/2X»£15.7) 

INVICID 

196 

600 

FDRflAT(//,20X, 'PROGRAM TERMINATED ON MEAN SORT ERROR CRITERION') 

INVICID 

197 

610 

F0RMaT(//>20X> 'PROGRAM TERMINATED ON PHYSICAL TIME '/ 

INVICID 

198 


1 'CONVERGENCE CRITERION') 

INVICID 

199 

620 

FORMAT! ///20X> 'PROGRAM TERMINATED ON LOCAL ERROR CRITERION') 

INVICID 

200 

630 

F0RMAn//»20X» 'PROGRAM TERMINATED ON LMAX ' ) 

INVICID 

201 


RETURN 

INVICID 

202 


END 

INVICID 

203 


SUBROUTINE DVISCP 

DVISCP 

1 


COMrtON/Fl/NXM/NXMl/NXM2»NXM3/NXMA/Nl,Ml/Nll/ Ml 1/ NS 1/ NS2/ NS» MS/ MSI 

DVISCP 

2 


1/NCRL/NCWL l»NCRL2/NCkLM/NCWLP/NSTRUT 

DVISCP 

3 


CUM.MuN/FA1/M12»M13/N12/N13/NXM5/NXM6/MS2/MS3/MSA/MSP1/MSP2/MS?3 

DVISCP 

4 


1/NSM1/NS1P1/NS2M1/NS1M1/NS2P1/NNS/NS1M2/NXM7 

DVISCP 

5 


CUrtM 0 N/F 2 /Rn( 61 > 55 ) / U ( 61/ 55 ) / V ( 61 / 55 ) / P ( 61/ 5 5 ) / T ( 6 1/ 55 ) / 

DVISCP 

6 


1 SH(61/55)/H(61> 55)/ ROS( 55) /US (55 )/VS(55)/ 

DVISCP 

7 


2 PS(55)/TS(55)/SHS(55.)/HS(55) 

DVISCP 

8 


C0MMUN/F3/VI SL ( 61/ 55 ) / VIST < 61/ 55 ) / VI S (61/ 55 ) / V I SS ( 55 ) 

DVISCP 

9 


C0MMaN/FA/SIGX(Sl,55)/SIGY(61/55)/TAoXY(61/ 5 5)/QX(61/55)/ 

DVISCP 

10 


1 OY(61/55)/SIGXS(55)/SIGYS(55)/TAUXYS(55)/QXS(55)/0YS(55) 

DVISCP 

11 


C0MM0N/F5/AU(61>55/A)/ AM(61/55/4)/AN(61/55/4 )/AUl(61/55/A) 

DVISCP 

12 


C0MM0N/F6/AUS( 55/4 )/AMS(55/ 4)/ ANS (55/4)/ AUSl (55/ 4) 

DVISCP 

13 


C0MM0N/F7/X(6l/55)/Y(6l/55)/XXI(61/55)/YXI (61/55)/ 

DVISCP 

14 


25 
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10 

15 


20 


30 


1 XETA(61^ 55)>YET4(61»55)»AJ(61/55) DVISCP 15 

Cbrtf1ON/F0/XS<55),YS(55)/A5XI(55)»YSXI{55)>XSETA(55)> OYISCP 16 

1 YSETA(55)f AJS (55) DVISCP 17 

CC«M0N/F9/INT(55)>TE1(55)»TE2(55)>TE3(55)#TEA(55) DVISCP 18 

CUt1.iQN/F13/GAMA,RGAS>CP^PR>PRT>CVIS> FM>PF>TFjUF>ROF,SHF,HF DVISCP 19 

V13(1>1JNXM)-VISL {l>lINXM)+VlST(l/lJNXh) DVISCP 20 

NXKrl-NXM-Ml DVISCP 21 

AU(l>2»ljNXMM)»U(l»2|NXM'1)-U(l/ljNXMM) DVISCP 22 

AU(l»lfljMl)-U(l»2}Ml)-U(l»l}Ml) DVISCP 23 

AU(l#2»2}NXMM)«V(l>2|NXhM)-V(l/lJNXM(1) DVISCP 2A 

AU(l>l>2jMl)-V(l»2jMl)-V(l>l}fU) DVISCP 25 

Ah<l>2»ljNXMM)-SH(l»2jNXMM)-SH(l>liNXMM) DVISCP 26 

AM(l/l»lJMl)«SH(l,2|Mi)-SH(l/ljMl) DVISCP 27 

AU(2>l»3;NXMl)«J(2»iINXP.l)-0(l>i;NXMl) DVISCP 28 

rtl(ljNl)«Q8VGArHR(U(l»l}NXrt)»lNT{l;Nl);TEl(l{Nl)) DVISCP 29 

TE2(UN1)«Q8VGATHR( U ( 2/ UNXM ) » INT ( 1 1 N1 ) J TE2 ( 1 J Nl ) ) DVISCP 30 

TE3(1JN1 )"TE2(ljNl)-Tei(ljNl) DVISCP 31 

AU(i>l>3 5NXM)»Q8VSCATR (TE3(1}ND# INT (1;N1 ) J AU( 1>1, 3ji'(XM) ) DVISCP 32 

AU(2>l#4}NXMl).V(2»lINXhl)-V(l>ljNXMl) DVISCP 33 

TEl(im)-38VGATHR(V(i»l}NXi1)> INTdjNDiTEK IJNI) ) DVISCP 3A 

TE2(1JN1)-Q6VGATHR( V ( 2> 1; NX M) > INT ( 1 1 N1 ) j TE2 (1 JNl ) ) DVISCP 35 

TE3(liNl)-TE2(ljNl)-TEl(lJNi) DVISCP 36 

AU(1»1» AJNXM )»38VSCArR(TE3(l JNl)^ INT (IJNI ) J AU(1# 1* 4JNXM) ) DVISCP 37 

Ah(2>l>2)NXHl)»SH(2/l;NXNl)-SH(l>l;NXrU) DVISCP 38 

TEl(liSl)»Q8VGATHR(SH(l»lJNXM)>INT(lJNl)jTEl(lJNl) ) DVISCP 39 

TE2( l}Nl)"0eVGATHR(SH(2»i;NXM)>INT(UNl) JTE2 (IJNl) ) DVISCP 40 

TE3(1;N1)-TE2(1;N1)-T£1(1JN1) DVISCP 41 

AM(ld»2)NXM)-38VSCATMTE3(liNl)»lNT (liNl ) )AP,(1/1»2JNXM ) ) DVISCP 42 

IE(NSTRUT.E3.0)Ga TD 15 DVISCP 43 

DL 10 N-NS1/NS2 DVISCP 44 

AL(f1S»N,3)-U(MSPl>N)-J(MS>N) DVISCP 45 

At(.’iS^N»4)«V(MSPl>N)-V(MS>N) DVISCP 46 

Ah(HS,N,2)-SH(iSPl»N)-SH(MS#N) DVISCP 47 

CONTINUE DVISCP 48 

AM(1>1»3;NXM )-( YETA(1» IjNXn) ♦AM{i#l/ 1JNXM)-YXI ( 1/ 1 1 NXM ) ♦ AM ( 1^ 1» 2i DVISCP 49 

INXM) )/AJ (1,1 JNXM) DVISCP 50 

Ah(l,l,4JNXM)-(XXI(l»lJNXM)*Art(l,l,2jNXI1)-XETA(l,l JNXM)*AM(l,l»l; DVISCP 51 

INXrt) )/AJ(1,1;NXM) DVISCP 52 

AN(1,1,1JNXM)-(YETA(1,1JNXM)*AU(1,1, ljNXM)-YXI ( 1, Ij NXM ) ♦ AU( I, 1, 3 1 DVISCP 53 

INXrt) )/AJ (1,1}NXM> DVISCP 54 

AN(l,l,2;NXM)«(XXI(l»l;NXM)*AU(l#l,3jNXM)-XETA(l#liNXM)*AU(l#l#li DVISCP 55 

INXM) )/AJ(1,1;NXM) DVISCP 56 

AN(i,l»3jNXM)»(YETA(l»ljNXr.)*AU{l,l,2iNXM)-YXI(l,lJNXM)*AU(l,l»4J DVISCP 57 

1NXM))/AJ(1,1}NXM) DVISCP 58 

AN(l,l,4;NXH)-(XXI(l,i;NXh)*AU(l,l/4;NXM)-XETA(l,l JNXM)*AU(1,1, 2; DVISCP 59 

INXM) )/AJ(l»llNXM) DVISCP 60 

SIGX(l,l)NXM)-P(l,ljNXMH-2.*VIS(l,l)Nxn)/3.*(AN(l, 1, 4JNXM)-2 . • OVISCP 61 

1AN(1,1,1JNXM )) DVISCP 62 

SlGY(l>l}NXM)-P(l,l;NXM)+2.*VlS(i>l} NXM) / 3 . * ( AN( 1, l>l;NXM)-2.* DVISCP 63 

1AN(1,1,4}NXM) ) DVISCP 64 

TAuXY(l,ljNXM)--VIS(l,lJNXM)*(AN(l,l>2lNXM)+AN(l,l,3iNXM)) DVISCP 65 

QX(l,ljNXi)--( VISL(l>lJNXM)/PR+VIST( 1,1JNXM) /PRT )*AM(1,1,3|NXM) DVISCP 66 

QY(1,1 JNXM)--( VISE (1,1; NXM)/PR+VI ST (1x1 iNXM) /PRT)*AM( 1,1,4} NXM) DVISCP 67 

If (NSTRUT.EO.O)RETURN OVISCP 68 

DO 20 N-1,NS OVISCP 69 

NN-N-l+NSl OVISCP 70 

A0S(N,3)-US(N)-U(MS1,NN) DVISCP 71 

AUS(N,4)-VS(N)-V(MS1,NN) DVISCP 72 

AMS(N,2)»SHS (N)-SH(nSl,NN) OVISCP 73 

DC 30 N»2,NS DVISCP 74 

A(^S(N,1)«US(N)-JS(N-1) DVISCP 75 

AUS(N,2)»VS(N)-VS(N-1) DVISCP 76 

AfiS(N,l)-SHS(N)-SHS(N-l) DVISCP 77 

AUS(l,l)-US(l)-U((1S,NSlfU) DVISCP 78 

AUS(1,2)-VS(1)-V(MS,NS1M1) DVISCP 79 

AMS(1,1)-SHS(1)-SH(MS,NS1M1) DVISCP 80 

ans(i,3;ns)« (yseta(i;ns)*ams(i,i;ns)-ysxi (1}NS)*ams(i,2;ns) )/ dviscp si 

lAJS(ljNS) OVISCP 82 
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10 

15 


AI'S(1#4}MS)» (XSXKll Mi )*AMS (l>2jNS)-XSETA(l}NS)*A«S(i#llNS» )/ OVISCP 8 3 

lAJSdlNS) DVISCP 84 

ANS(l>ljNS)-(rSETA(li'<S)*AUS(l»liKS)-ySXI(l)NS)*AUS(l»3}NS) )/ DVISCP 8 5 

lAJS(llMS) DVISCP 86 

ANS(1/2JMS)« (XSXKDNS )*ADS( l/35NS)-XSETA(l)NS) + AJS(ldJNS) )/ DVISCP 87 

lAJSdlNS) DVISCP 88 

ANSd>3}MS)« (rSETA(lJNS)*AUSd»2jNS)-riXHliHS)*AUSd/4}NS) )/ DVISCP 89 

lAJSdJNS) DVISCP 90 

ANS (lj4;NS>»(XSXId}NS )*AUS ( 1/ 4; N S ) - XSETA ( 11 NS ) *ADS d> 2 1 NS ) )/ DVISCP 91 

lAJSdlNS) DVISCP 92 

SIGXS (IJNS )-PS ( 1} NS )+2.*VlSS dlNS )/ 3 . *{ ANS ( 1 >41 NS ) -2 . *ANSd/ 1 J NS ) )DVISCP 9 3 

SlGYSdjNS)«PS(llNS) + 2.*VlS j(llNS)/3.*(ANS(l>llNS J-2.*ANS(l»4lNi))0VlSCP 94 

TAJXYSd)NS)--VISS(i;NS)*(ANSd>2}NS)+ANSd»3lNS)) DVISCP 95 

aXSdlNS)- -VISSd}NS)/PR*AMSd»3iNS) DVISCP 96 

DYSdlNS)- -VI3Sd;NS)/Pk*AriS(l>41NS) DVISCP 97 

RETURN DVISCP 98 

END DVISCP 99 

SUBROUTINE OVISCC DVISCC 1 

CUMi 1UN/F1/NXM> NXM1> NXn2>NXK3>NXP4>Nl>Ml>Nll» Ml 1> NS 1> NS 2> NS> MS> M S 1 DVISCC 2 

1>nCWL>NCWL1>NCVL2>NC(<lM»NCwLP>NSTRUT DVISCC 3 

CCiN(10N/FAl/M12>M13>N12>N13,NXM5>NXtl6>MS2>MS3>MS4>MSPl>MSP2»MSP3 OVISCC 4 

l>NSrtl>NSlPl>NS2Ml»NSlMl»NS2Pl>NNS>NSlM2>NXK7 DVISCC 5 

CUf1rtON/F2/RQ (61> 55)>U(61>55)»V(61>55)>P(61>55»>T(61>55)> OVISCC 6 

1 SH(61>55)»H(61»55)/R0S(55)/US(55)>VS(55)> DVISCC 7 

2 PS(55),TS(55)>SHS(55)>HS(55) DVISCC 8 

C0i1i10N/F3/VlSL(61»55)>VlST(61>55)>VIS(61»55)»VISS(55) OVISCC 9 

CDrt'10N/c4/SIGX(61»55>> S1GY(61> 55 ) >T AUXY ( 61> 5 5) > OX ( 61> 55 ) > OVISCC 10 

1 CiY<61»55)>SIGXS(55)>SIGYS<55)>TAUXYS(55)>QXS<55)>0YS(55) OVISCC 11 

CON.NON/F5/AU(61» 55> 4 ) > AM ( 61/ 55> 4 ) > AN ('61> 5 5> 4 )> AUl ( 6l> 55> 4 ) DVISCC 12 

COrtliON/F6/AII5(55/4)/AMS(55/4)/ANS(55/4)/AUSl(55/4) DVISCC 13 

CuMrt0N/F7/X(61>55)/Y(61/55)/XXH61/55)/YXI(61/55J/ OVISCC 14 

1 XtTA(61/55)/YETA<61/P5)/AJ(6i/55) OVISCC 15 

CaH.1GN/F8/XS (55), YS(55 )/XSXI (55)>YSX1(55)/XSETA(55 ), OVISCC 16 

1 YSETA(55)/AJS(55) OVISCC 17 

C0.Nrt0N/F9/INT(55),TEl(55),Te2{55)/TE3(55)/TE4l55) DVISCC 18 

CUM.1DN/F13/GAMA,RGAS,CP,PR,PRT/CVIS/FM/PF/TF/UF/R0F/SHF/HF DVISCC 19 

VISd/llNXM)-VISLd/llNXH)+vIST(l/llNXM) OVISCC 20 

NX(1(1»NXM-H1 DVISCC 21 

AUd/l/llNXMM)-Ud/21NXMN)-Ud/llNXMM) OVISCC 22 

AUd/Nl/llMl )-Ud,Nllhl)-Ud/Nll}Ml) DVISCC 23 

Acd/l/2lNXMM)-Vd/2;NXMh)-Vd/llNXMM) DVISCC 24 

AUd/Nl/2lMl)«Vd,NlJMl)“Vd#NlilMl) OVISCC 25 

AMd/Nl/llhl)«SH(l/NllMl)-SHd/Nll)Ml) DVISCC 26 

Ahd/1/llNXMH) ■SHd/2lNXMM)-SHd/llNXMM) OVISCC 27 

AUd/l/3)NXMl)-J(2/llNXMl)-Ud/l|NXt11) OVISCC 28 

TEldlND-aeVGATHRCUdil/llNXlD/INTdlNDlTEldJNin DVISCC 29 

TE2dlNl)-Q8VGATHR( U( Ml , 1 1 NXM ) , INT dlNl ) 1 T E2 dlN 1 ) ) DVISCC 30 

TE3(llNl)«TE2d)Nl)-TEldlNl) OVISCC 31 

AU(Ml,l/3lNXM)«Q8VSCATR(TE3(llNl)/INT(llNl)lAU(Ml/l/3)NXrt)) DVISCC 32 

A0d/i/4lNXMl)-V(2,ilNXMl )-Vd/l)NXMl) OVISCC 33 

TEldlNn-QeVGArHR(V(Mli,lINXM)/INTd)Nl)lTEldlNl)) OVISCC 34 

TE2dlNl)«Q8VGATHR ( V( Ml / 11 NXM) / INT d INI ) ) T E2 (llNl ) ) OVISCC 35 

TE3dlNl)-TE2d;Nl)-TEl(llM) DVISCC 36 

AL(Ml/l/4}NXM)-Q8VSCATR(IE3dlNl)/lNT(llNl)l AU(Ml/l/4iNXM)) DVISCC 37 

AM(l/l>2}NXMl)»SH(2/llNXMl)-5Hd/llNXMl) DVISCC 38 

TEldlNl)-Q8VGATHR(SH(Mll,llNXh) / INT (llNl )lTEldlNl) ) DVISCC 39 

TE2(1)N1)-Q8VGaTHR(SH( Ml, 1 1 NXM )/ INT dlNl ) iT E2 (1 |N1 ) ) OVISCC 40 

TE3(llNl)-TE2d)Nl)-TEl(l)Nl) DVISCC 41 

Ah(ru,l,2;NXM)-O0VSCATR(TE3(llM),INT(llNl)lAM(Ml,l,2lNXM)) OVISCC 42 

IF(NSTRUT.EQ.0)G0 to 15 DVISCC 43 

0C( 10 N-NS1/NS2 DVISCC 44 

NN-N-NSl+1 DVISCC 45 

AU(1S1/N,3)«US(NN)-U(MS1,N) DVISCC 46 

AU(MSl/N/4)-V3 (NN)-V(MS1,N) DVISCC 47 

AM(MS1/N,2)«SHS(NN)-Sri (MSl/N) OVISCC 48 

continue DVISCC 49 

AMd,l/31NXH)-( YETAd/l}NXM)*AM(l/l/ 11NXM)-YX1 d/llNXM)*AM( 1/1/21 OVISCC 50 

1NXN))/AJ(1,1)SXM) OVISCC 51 
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AhCl#l>4;NXM)-(XXI (1*1 JNXi1)»AM(1,1,2}NXM)-XETA(1>1 DVISCC 

lNXrt))/AJ(l»l;NXM> OVISCC 

An( 1/1,1;NXM)-( YETA(l»l|NXh)*AU(l>l# liNX«)-YXl (1>1 jNXM)*AU( OVISCC 
INXfl) )/AJ (1/ljNXN) OVISCC 

AN(iil#2}NXM)»(XXni>I}NXf1)*AU(l>l/3}NXM)-XETA(l>liNXh)*AU(l>l/lJ OVISCC 
IMXfl))/AJ(l,l}NXM) OVISCC 

AN(l>l,3)NXH)»(YETA<l»liNXMI»AU(l/l>2iNXM)-YXI(l>l jNXh)»AU(l/l» A) OVISCC 
iNXiDj/AJCl^UNXM) OVISCC 

AK(l»l*A5MXM)-{XXI(l*UNX«)*AU(l>l/A;NXM)-XETA(l/liNXM)*AU(l»l»2j OVISCC 
1NXM))/AJ(1,1;NXM) OVISCC 

SlGXa^ 1}NXM)»P( l»lJNXM) + 2.*VIS(l>ljNXM)/3.* (AN(1# l/A|NXrt)-2.* OVISCC 

1AK(1#1»1JNXM ) ) OVISCC 

S1GY(1»1 jNXM)-Pa/l;NXh»+2.*VlS(l/liNXM)/3.* (AN(1/ l/l;NXM)-2.* OVISCC 

1AN(1,1,A}NXM)) OVISCC 

TAJXY(1»1}NXM)--VIS(1»1JNXM)*(AN(1>1#2JNXM)+AN(1#1/3;NXM)) OVISCC 

QX(l>ljNXM)--( VISL (l»ljNXM)/PR+VIST( 1>1JNXM) /PRT )♦ AM ( 1/ 3J NXH ) OVISCC 

QY(l^l}NXM)«-( VISL (l/l;NXrt)/PR+VlST< l»l;NXM) /PRT ) *AM( 1>A; NXM) OVISCC 

IF(KSTRUT.EO.O)RETURN OVISCC 

00 20 N-l^NS OVISCC 

NK-S-1+NSI OVISCC 

AUS(N»3)-US(N)-U(MS1»nN) OVISCC 

AU5<N>A)-VS(N)-V(MS1>NN) OVISCC 

20 AhS(N>2)-SHS (N)-SH<MS1/NN) OVISCC 

00 30 N»1,NSM1 OVISCC 

A0S(N,1)-US{N+1)-US(N) OVISCC 

AUi(N»2)-VS(N+l)-VS(N) OVISCC 

30 AMS<N»1)-SHS(S+1)-SHS(N) OVISCC 

A0S(NS»1)-U(HS>NS2P1)-US(NS) OVISCC 

A0S(NS>2)-V(MS>NS2P1)-VS(NS) OVISCC 

AMS < NS » 1 )-SH(MS>SS2Pl)-SHS(NS) OVISCC 

AMS{1,3;NS)- (YSETA(1;NS)*AMS(1>1JNS)-YSXI ajNS)*ArtS(l>2jNS) )/ OVISCC 

lAJS(liNS) OVISCC 

AhS(l^A}NS)«(XSXI(lJNS ) ♦ AMS ( 2} NS )-XSET A ( 1 } NS ) ♦ AMS ( 1/ 1 J NS ) )/ OVISCC 

lAJS(ijNS) OVISCC 

ANS(1>1JNS )■ (YSETA( 1JNS)*AUS(1>1JNS)-Y5XI (l; NS)*AUS(i>3 JNS) )/ OVISCC 

lAJS(liNS) OVISCC 

ANi(l/2}NS)-(XSXI(ljNS )*AUS ( 1# 3JNS )-XSETA (1 J NS ) *AJS ( 1» UNS ) )/ OVISCC 

1AJS(1;NS) OVISCC 

AN5(1> 3jNS )«(YSETA(l}NS)i‘AUS(l>2;NS)-YSXiaiNS)'t'AUS(l>AjNS» )/ OVISCC 

1AJS(1;NS) OVISCC 

ANS< A;*JS )• (XSXI ( liNS )*AuS( 1, AIN5)-XSETA(1J NS ) ♦ AU S ( 2 > NS ) ) / OVISCC 

1AJS<1;NS) OVISCC 

SlGXS(l;MS )-PS ( 1;NS)+2.^VI:.:) (liNS )/3.*(ANS(l#4|NS)-2.+ANS(lil)NS ) )DVISCC 
SlGrsaiMS)«PS (1 jNS)+2.*VlS:,(liNS )/3.^(ANS (1#1JNS)-2.*ANS(1#4;N5) )DVISCC 
TAUXYS(l;MS)-“VISS(liNS )*(ANS(1/2;NS )+ANS(l# 3;NS) ) OVISCC 

QXS(l;SS)« -VI.SS (ljNS)/PR*AMS(l>3jNS ) OVISCC 

QYS(ljNS)- -VISS (ljNS)/PR*AMS(1^4JNS ) OVISCC 

RETURN OVISCC 

END OVISCC 

Subroutine veci veci 

CUrli10N/Fl/NXM^NXil^NXrt2/NXn3/NXh4f Ml 1 # NS NS2^ NS# MS# MSI VECI 

1#NCRL#NCWL1#NCWL2#NCWLM#NCWLP#N$TRUT VECI 

COMMON /F 2/ RO (61, 55 ) , U ( 61 # 55 ) , V ( 61 , 55 ) # P ( 61 , 5 5 ) # T ( 6 1# 55 ) # VECI 

1 SH{61#55)#H(61#55)#RG5{55)#US(55)#VS(55)# VECI 

2 PS(55),TS(55),SHS(55),HS(55) VECI 

COMiON/FA/SIGX(61# 55 ) # S I G Y ( 6 1# 5 5 ) # TA UXY ( 61# 5 5 ) # OX ( 61# 55 ) # VECI 

1 QYC61#55)#SIG)(S( 55)#sIGYS(55)#TAUXrS(55)#QXS{ 55)#QYS(55) VECI 

C0«M0N/F5/AU(6l, 55 , 4 ) , AM ( 61 , 5 5, 4 ) , AN ( 61# 5 5# 4 )# AU 1 ( 61# 5 5# 4 ) VECI 

CCiMMGN/F6/AUS( 5 5#4)#AMS(55#4),ANS (55#4)#AUS1 (55# 4) VECI 

CurtM0N/F6/AUS(55,4),AMSt55,4)#ANS (55#4)#AUS1 (55#4) VECI 

CaMfl0N/F7/X(61#5 5),Y(61#55)#XXl ( 6 1 # 5 5 )# YX I( 6 1# 5 5 ) # VECI 

1 XETA(61, 55) ,YETA(61,55), AJ( 61#55) VECI 

CGMM0N/F8/XS (55)#YS(55)#XSXI ( 5 5 ) # YS X I ( 55 ) # XS ETA ( 55 )# VECI 

1 YSETA(55),AJS(55) VECI 

AN(l#l#4jNXM)-R0(l#liNXh)*H( l#l;NXM)-P(i#ljNXM) VECI 

AU(l#l#ljNXh )-R3( 1#1;NXM)*U( 1#15NXM) VECI 

AU(1#1#2;NXM)-AU(1#1# l;NXM)*U(l# l;NXM)+SIGX (1#1;NXM) VECI 

Au(i,i, 3 ;nxm)»a'J( 1,1, i;nx«)+v(i, unxm)+taux y( i# i;nxm) veci 


52 

53 

54 

55 

56 
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71 
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95 

96 

97 

98 

99 
100 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

10 

11 

12 

13 

14 

15 

16 

17 

18 


28 



APPENDIX A 


JNXM)»(AN(1, 1> ^)NX)i)+SIGX(l^l>NXM) ) ♦U(ljl;NXM)+V(l>l|NXr1)VECl 


l^TAJXYd# 15NXM)+3X (1# i;NXrt) \JEC1 

A^(i#l#l^NX^< )-R3(l/l;^IXH)+V( 1#1;NXK) VECl 

AN(l,i^2jNXM)*AN(l,lj i;NXM)AU(ld5NXM)+TAuXY(id5NXi) VECl 

AN(l>i,3 JNXM )«AS( l>ld jNXM) »V( 1#1 JNXM) + SIGY( liNXM) VECl 

AN(ljl#^;SXM )■ (A^(l,d4;NXh)-^SIGY(ld)NXM) )^V(ldiNXf1)+U(l#l5NXM) VECl 
1+TAJXY( 1# 1 ;NXM)+QY(li 1;NXM) VECl 

uO 20 I«l>4 VECl 

AKddd )SXK )-YETA(l,l)NXn)+ AU ( 1 dd JNXM ) -X ET A C 1 d ; NXM ) +AN ( Id d i V EC 1 
INXM) VECl 

ANd,ld ;NXM)--YXI(1/1 jNXM)* AUddd;NXM)+XXIdd;NXM)+AN( Idd ; VECl 
INXrl) VECl 

20 CONTINUE VECl 

Atd^l>l;NXM )-RDddiNXK)»AJ d#ljNXM ) VECl 

ALd/ldJNX^' )«AUd#ld;NXrt)»Ud>i;NXM) VECl 

AUddd;NXK)»AJ( ldd;NXH)*VddjNXh) VECl 

AU(ldd;NXM )■ (ROdd;NXM)*H(ld|NXN )-P (Id i NXM ) ) *A J (Id ;NXM ) VECl 

IF(NSTRUT,EQ.O)RETURN VECl 

DO 30 N-l/NS VECl 

ANS(Nd) -RDS (N)tHS(N)-PS(N) VECl 

AOS(Nd)»R0i (N)*US(N) VECl 

AUS(Nd)-AUS (Nd)^US(S )+SIGXS(N) VECl 

AO$(Nd)*AUS (Nd)*VS(N)+TAUXYS(N) VECl 

AOS(Nd )• ( ANS(N#4) ♦SIGXS (N) )*US (N )♦ V S (N ) AU XYS ( N ) +QXi> ( N ) VECl 

AN:^(Nd)-ROS(N)^VS(N) VECl 

ANS(Nd)«ANS (Nd )*US (N)+TAUXYS (N) VECl 

ANS(N>3)-ANS (Md )^VSiN)+bIGYS(N) VECl 

ANS(Nd)- (ANS('(d)4SlGYii (N) )^VS(N)-^US(N)tTAUXYS(N)+QYS(N) VECl 

30 CONTINUE VECl 

Du AO I«1^4 VECl 

Du AO N«1^NS VECl 

Ah:>(Nd)«YSETA(N )*AUS(N>1 )“XSETA(N)»ANS(N^I ) VECl 

AN$(Nd ) --YSXI (N)^AUS(Nd ) + XSXI (N)*ANS(N> I ) VECl 

AO CONTINUE VECl 

00 50 N*1#NS VECl 

AbS(Nd)«R3S(N)<'AJS(N) VECl 

AOS(Nd )-AUS (Nd )*US (N) VECl 

AUG(Nd)-AUS (Nd)+Vb(N) VECl 

AUi(Nd)- (P0S(N)*HS(N)-PS (N) )*AJS (N) VECl 

to CONTINUE VECl 

RETURN VECl 

ENO VECl 

jUdROuTiNE IVECl IVECl 

CuiiNuN/Fl/NXMdXMdNXNa# NXh3 ^ NXMA#N1 /rtl^ N 11# Mil# NS 1#NS2 #nS# MS# MSI IVECl 
1 #nCWL#NCWL 1# NCWL2#NCWLM#NCWL P#NSTRUT IVECl 

C0.1rtON/F2/RO(61# 55)#U(61#55)#V(6l#55)#P(61#55)#T(61#55)# IVECl 

1 SH(6l#35)#H(61#55i#R0S(55)#uS (55)# VS (55)# IVECl 

2 PS(55)#TS( 55)#SHS(55)#HS(5p) IVECl 

C0NNDn/F 5/AU(61#55#A)# Ah(61# 55# A ) # AN ( 61# 5 5# A )# AUl ( 61# 5 5# A ) IVECl 

C0M.i0N/F6/ AUS(55# A )# AiS ( 5 5# A ) # ANS ( 55 # A) # AUS 1 ( 55# A ) IVECl 

CUrtM0N/F7/X( 61#55)# Y(61#55)#XX1(61#55)#YXI(61#55)# IVECl 

1 XETA(61#55) ,YETA(61#55)#AJ (61#55 ) IVECl 

COMiON/FB/XS (55)# YS(55 )#XSX: ( 5 5 ) #YSX I ( 55 ) # XS ET A ( 55 ) # IVECl 

1 YSETA(55)#AJS(55) IVECl 

AO(l#l#l;NXM )*ROd, l;NXM) *U( 1#1;.NXM) IVECl 

AOd#l, 2 JNXM)-AU(1#1# ljNXM)^Ud#liNXM) + P(l# UNXM) IVECl 

AU( 1#1# 3;NXM)-AU(1#1# 1;NXM)^V(1# DNXM) IVECl 

AJd#l#A;NXM)-AJ(l#l#l;NXM)^Hd#liNXM) IVECl 

ANd#l#UNXM)*Ra(l#l;NXH)+V(l#lJNXM) IVECl 

AN(l#l,2nxM )-AS(1#1#1;NXM)*U( 1 # 1 )NXM) IVECl 

AN(l#l#35NXM)-ANd#ld;NXM)*Vd#l;NXM)+P(l#l JNXM) IVECl 

AN(1#1# A)NXM)-AN(l#l#liNXrt)^Hd#liNXM) IVECl 

DO 20 I-1#A IVECl 

AKd#l#I ;NXM )»YETAd# i;fiXM)^ AU d # 1 # I )NXM ) -X ET A ( 1 # 1| NXfl) +AN ( 1# 1 # I )I VEC 1 
INXM) IVECl 

AN(1#1#I JNXM )»-YXI (1#1JNXM)^ AU (1 # 1# I )NXM ) +X XI ( 1 # 1 j NXM ) ♦ AN ( 1# 1# I ; IVECl 
INXM) IVECl 

20 CONTINUE IVECl 
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30 


40 


50 


AL(i>l»i;NXh)«RO(l/l;NXH)*AJ(i,ljNXM) 

Ab(l>l,2jMXM)»AU(l>l>l;NXHH‘U(l>ljNXM) 

Al(l>l»3jNXM)«AU(l#l#l;hxM)*V(l»ljNXM) 

AU(1^1>4;NXM)-(R3{ 1> liNXM)+H(l>liNXr NXM) ) ♦ A J ( 1/ 1 JNX M ) 
n (NirRUr.EQ.OjRETURN 
DO 30 N»1,NS 
AUi(N#l)«ROS(N)*US(N) 

AUa(N>2)»AUS(N>1) *03 (:S ) + PS (N ) 

AUj(N, 3)«AUS (N>1)*VS(N) 

AI<5(N>4)«AUS(N>1)*HS(N> 

AN3(N,1)«R0S(N)*VS(N) 

AN3(N<2)«ANS (N/1)*0S(N) 

ANs(N>3)-ANS (N# 1 )*VS(N)+PS(H) 

ANS(N<4)»ANS (N,1)*HS(N) 

continue 

OCI 40 I«l,4 
DC 40 N»1>NS 

AMi(N^I)»YSETA(N) *AUS (N> I)-xSETA(N)*ANS(N>I) 

AKS (N» I )«-rSXl (N)*AUS(N^1 )+XSXl (N )fANSlN> I ) 


IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 

IVECl 


CONTINUE IVECl 

DO 50 N"1»NS IVECl 

AUS(N,1)»R0S (N)*AJS(N) IVECl 

AUS(N/2)*4US(N>1)*US(N) IVECl 

AU5(N,3)»AUS(N,U*VS(N) IVECl 

AU5(N>4)»(R0S(N)*HS(N)-PS(N) )*AJS(N) IVECl 

CONTINUE IVECl 

RETURN IVECl 

END IVECl 

SUBROUTINE DAMPP DAHPP 

C0MrtaN/Fl/NXM>HX?11^NXrt2>NXK3^NXM4f Nl#hi#Nll> h^ll^NS 1#NS2#NS# MS/ MSI OAMPP 
l#NCWL/NCWLl/NCWL2/NCWLM/NChlP/NSTRUT OAMPP 

C0MMDN/FA1/M12# M13/N12/N13/NXM5/NXM6/MS2/ MS3/MS4>jiSPi/MSP2/MSP3 OAMPP 
l/NSMl/NSlPl/NS2Ml/NSl.il/NS2Pl/ NNS/NS1M2/NXM7 OAMPP 

C(jMM0N/F2/R0(61, 55)/U(61, 55 ) / V ( 61 / 55 ) / P ( 61 / 5 5) / T ( 6 1/ 55 ) > OAMPP 

1 SH ( 61/ 55 ) / H ( 61/ 5 5 ) / R OS ( 5 5 ) / US { 55 )/ V S ( 55 } / OAMPP 

2 PS( 55)/TS (55) /SHS (55)/HS (55 ) OAMPP 

C0M.10N/F4/SIGX(61/55 ) / S IG Y ( 6 1/ 55 ) / TAUXT ( 61 / 5 5) / QX ( 61/ 55 ) / OAMPP 

1 QY(61/55)/SIGXS(55)/SIGYS(55)/TAUXYS(55)/QXS(55)/ QYS(55) OAMPP 

C OMMON/F 5/ AU (61/55/4)/ Ah l61/55/4)/AN( 61 /55/4)/AU 1(61/55/4) OAMPP 

CDMM0N/F7/X(61/55)/ Y (61/ 55)/ XXI (61/55 )/YXl (61/55)/ OAMPP 

1 XETA (61/ 55 ) /YETAC 61/55 )/AJ (61/55 ) OAMPP 

CUMrtUN/F10/DT( 61/ 55 )/ DTM OAMPP 

COMfiON/F12/LSYM/NFLQW/ LFI/LFU/ ILT/LMAX/L W/ L E XI T/ CR I T/ FOT/ FOIL/ OAMPP 

1 TIME/CCP/CCT/L/LK/TIhEL/CRITAVG OAMPP 

CuMMdN/F13/GAMA/RGAS/CP/PR/PRT/CVIS/FM/PF/TF/UF/R0F/5HF/HF OAMPP 

QX(i/ 1|NXM)»GAMA*RGAS+T(1/1 ;nXM) OAMPP 

QX(l/l;NXh)«VS3RT(UX(i/ UNXh ) ) CiX ( 1/ 1 ;NXM ) ) OAMPP 

SlGX(l/i;NXM)*(VA3S(U{l/l;NXM)/SIGY(l/lJNXM))+QX(l/ijNXM)) OAMPP 

SlGY(2/2|NXM4).p(2/3jHXh4)-2,*P(2/2|NXM4)+P( 2/1JNXM4) OAMPP 

SIGY(2/2;NX14)-VABS(SIGY( 2/2/NXM4);S 1GY(2/2J NXM4) ) OAMPP 

UY(2/2;NXM4)-T(2/3;nXM 4)-2.*T(2/2;NXM4)+T(2/ l;NXM4) OAMPP 

QY(2/2jNXM4)-VABS(QY(2/2iNXM4)|QY(2/2|NXM4)) OAMPP 

;ilGX (2/2)NXM4)-S IGX(2/ 2iNXM^)»(SIGY( 2/2INXM4 )*CCP/ OAMPP 

1( P(2/3,-NXM4) +2. + P(2/2;NXM4) + P( 2/l,*NXM4) )+QY( 2/2;NXM4)^CCT/ OAMPP 

l(T(2/3;NXM4)+2.+T(2/2iNXM4)+T(2/l)NXM4) ) ) OAMPP 

QX(l/l)NXM)«(VABS(V(l/l;NXM);SlGY(l/liNXM) )+GX(1/1;NXM) ) OAMPP 

SI6Y(2/2;NXM4)-P(3/2;NXM4)-2 .+P (2/ 2; NXM4 ) +P( 1/2;NXM4) OAMPP 

$IGY(2/2jNXM4)-VABS(SIGY(2/2;MXM4);SIGY(2/2; NXM4)) OAMPP 

QY(2/2iNXM4)-T(3/2lNXM4)-2.»T(2/2|NXM4)+T( 1/2JNXM4) OAMPP 

CiY(2/2|NXM4)-VABS (QY(2/2jNXM4) iUY (2/ 2)NXM4) ) OAMPP 

SIGY(2/2;NXM4)- GX ( 2/ 2;NXM4 ) + ( S I GY ( 2/ 2 ; NXM4 ) *00 P/ OAMPP 

1 (P(3/2jNXM4)+2.tP(2/2iNXM4)+P(l/2iNXM4))+0Y( 2/2iNXM4)*CCT/( OAMPP 

2 T(3/2)NXM4) +2.^T(2/2;NXrt4)+T(l/2;NXM4) ) } OAMPP 

IMNSTP.UT.EO.O)GO TO 35 OAMPP 

00 30 N-NS1/NS2 OAMPP 

NN-N-NSl+1 DAMPP 


SIGY(MS1/N)-QX(MS1/N)+ (ABS(PS(NN)-2,*P(MS1/N )+P(MS2/N) )»CCP/(PS(NNOAMPP 
l)+2,+P(MSl/N) + P(MS2/N) ) + A BS ( TS ( NN ) -2 . *T ( MS 1/ N ) +T ( M S2/N ) ) +CC T/ ( T S ( OAMPP 
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2NM + 2.*T(M51/N)+T(hS2,Nl) 1 

DAMPP 

41 

30 

CONTINUE 

OAMPP 

42 

35 

CONTINUE 

DAhPP 

43 


OL 100 I-1/4 

OAMPP 

44 


UX(2,2;NXM4 ) -SIGX(2,2;NXh41* ( AU ( 2 , 2 , I ; NX M4 ) - AU ( 2, 1, I JNXM4 ) ) 

DAMPP 

45 


aY(2/2;NXM4)»SIGY(2/2;NXM4)* ( AU (2, 2, 1 J NXM4 ) - AU (1/2/ 1 JNX«4 )) 

OAMPP 

46 


0(j 10 N«2,N11 

OAMPP 

47 


ax(i,N)«o. 

DAMPP 

48 


QX(rtl,N)«0. 

OAMPP 

49 


QY(1,N)»0. 

DAMPP 

50 

10 

aY(rtl,N) -0. 

OAMPP 

51 


IF(NSTRUT.EO.O)GO TO 25 

OAMPP 

52 


OU 20 N-MS1/NS2 

DAMPP 

53 


aX(rtS,M)»0. 

OAMPP 

54 

20 

0Y(rtS,N)»0. 

DAMP? 

55 

25 

continue 

OAMPP 

56 


Ah(2/2,I ;NXM4) -Art (2/2, IjNXhA )-QX( 2/2 ;NXM4) /XXI (2/2iNXM4) 

DAMPP 

57 


AN(2,2/I;NXh4)-AN(^/2/ IjNXhA )-QY( 2/2 ;NXM4)/YETA(2/ 2 ;NXM4) 

OAMPP 

58 

100 

CONTINUE 

OAMPP 

59 


RETURN 

OAMPP 

60 


ENO 

OAMPP 

61 


SUBROUTINE OAMPC 

OAMPC 

1 


C0rtrt0N/Fl/NXM,NX'U/NXrt2/NXh3/NXrt4/Nl/Ml/Nll/ Ml 1/ Nil/ NS2 , NS, MS, MS 1 

OAMPC 

2 


l»NCWL>NCWLl,NC^L2,NCWLM>MCwLP»NSTRUT OAHPC 

CGrt;l0N/FAl/M12»M13>N12>M3>NXK5/NXM6i«S2> MS3 #MSA».'1bPl»MSP2» MSP3 OAhPC 
l»NSii>NSlPl»NS2'U»NSlnl/(4S2Pl#NNi>NSlM2>,SXh7 OAhPC 

C0rtrt0N/F2/R0(61>55)>U(61»55)> V(61>55 )/P(61#55)»T(61#55)> OAMPC 

1 SH(6i»55)>H(61>55)>ROS(55),US(55),VS(55), OAhPC 

2 PS(55)»TS(55)>SHS(55)/HS(55) OAhPC 

C0rth0.^/FA/SIGX(61»55)»SIGY(61#55J#TAUXY(61/5 5)#QX(61#55)» OAhPC 

i QY(61>55)»SIGXS(55)>SIGYS(55),TAUXYS(55),QXS(55»»QrS(55) OAhPC 

CCJhhuN/FS/AU (61» 55> A)» Ah( 61/ 55>A) >AN (61»5 5, A )» AUl(61>i)5»4) OAhPC 

CljhrtON/F6/AUS( 55/ A)/ Arts (5 5/4)/ AhS (55/41/ AUSl (55/4) OAhPC 

C0rtrt0N/F7/X(61/55)/Y(61/55)/XXI(61/55)/YXl(61/55)/ OAhPC 

1 x£TA(bl/55)/YETA(61/55)/AJ(61/55) OAhPC 

C(Jrtrt0N/F10/DT(61/55)/0Th OAhPC 

C0hrtaN/F12/LSYrt,rtFL0W/ LF 1 / LF 0/ ILT /L h AX , L W/ LE XIT/ CR IT, FDT »FDTL/ OAHPC 

1 TIrtE,CCP/CCT/L/LK/TIrtEL,CRlTAVG OAhPC 

C DrtrtOrt /F 13/GAMA, RGAS/ CP/ Pk/PRT/C VIS, FM/PF/TF /OF /ROF/ SHF /HF OAhPC 

QX(l/ljMXM)»GAMA*RGAS*T(l,l;NXH) OAhPC 

QXd/ijNXM) «VS3RT(0X(1, IJ nX. 1) ; 0X( 1,1 ;NXrt) ) OAhPC 

SlGXd, liMXh 1- ( V4BS(U( l,liNXM) |SIGY( l,ljNXh) )♦ QX <1 ,1 J NXM ) ) OAhPC 

S16Y(2,2jNXM4)-P(2,3)NXM4)-2.*P(2,2;NXh4)+P( 2/ IjNXhA) OAhPC 

SIGY(2,2;NXM4)-YABS(S1GY(2,2JNXM4) jSIGY(2,2jNXM4)) OAhPC 

aY(2,2)MXM4)«T{2/3|NXh4)-2.t‘T{2/2)NXh4)+T( 2, IJhXhA) OaMPC 

QY(2,25NXM4)-VABS (UY (2,2; NXH4) J0Y(2, 2;NXM4 ) ) OAhPC 

S16X(2/2;NXh4)-SIGX(2,2;NXh‘t)*(SIGY( 2,2;NXh4 )*CCP/ OAhPC 

l(F(2,3jMXM4>+2.+P(2,2iNXh4)+P(2,lJNXM4))+JY(2, 2iNXh4)*CCT/ OAHPC 

l(T(2,3;NXrt4)+2.*T(2/2;NXh4)-t-T(2/l;NXM4) 1 1 OAhPC 

axd/i;NXrti - (VA8s(v{i,i;Hxh) ;siGYd, i;nxm) )+cx( i,i,*NXrt) ) oahpc 

SIGY(2,2;NXh4)«P(3,2;rtXM4)-2.*P(2/2;NXn4)+P( l/2jNXh4) OAHPC 

SlGY(2,2;NXK4)«VAdS(SIGY(2,2jHXH4) jSIGY(2,2;NXh4)) OAHPC 

aY(2,2;rtXH4)«T(3,2;NXM4)-2.*T(2/2JhXh4)+T( 1,2 ;NxH4) OAHPC 

QY(2/2 jNXM4)«VABS(tiY(2,2;NXh4) ;0Y(2, 2 JNXM4)) OAHPC 

SlGY(2,2;rtXM4)- QX ( 2, 2JNXH4 ) ♦ (S I GY( 2/ 2 )rtXh4 )*CCP/ OAHPC 

1(P(3/2JSX14)+2.»P( 2/2iNXM4)+P( l,2)NXrt4) )+QY( 2,2)NXh4)*CCT/( OAHPC 

2 T(3/2;NXrt4)+2.*T(2/2;NXh4;+T(l/2;NXh4))) OAHPC 

IF(rt5TRUT.E0.0)G0 TO 35 OAMPC 

Ou 30 N-NSl, NS2 OAMPC 

■NN-N-NSl + 1 OAMPC 

SIGY(MSl,M)«QX(rtSl,N)*(A3S(PS(NN)-2.*P(rtSl/N)+P(hS2/N) )*CCP/(PS(NNOAhPC 

l)+2.*P(hSl,N )+P(rtS2,N) )+ABS(TS(NN)-2 .♦KrtSl, N)+T(rtS2,rtn*CCT/(TS ( OAMPC 
2Nh)+2.*T(MSl,N)+T(MS2,N)) ) OAHPC 

CCirtTiNUE OAMPC 

CONTINUE OAMPC 

DO 100 I»l/4 OAMPC 

aX(2,2;NXM4) -SIGX(2,2;NXM4)*(AU(2,3, 1 ; NXM4 ) - AU ( 2, 2 , I ;NXM 4) ) OAMPC 

aY(2/2;NXM4)-SIGY(2/2jNXrt4)*(AU(3,2,I;NXrt4)-AU(2/2/IiNXM4)) OAMPC 


31 
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10 


20 

25 


100 


DO 10 N»2,N11 

UX(1>N)«0. 

ax(m/N)«o. 

QY(l>N)-0. 

UY(i'll^N)»0, 

IF (NSTRUT.EO.O)GO TO 25 
00 20 N«NS1>NS2 
NN-N-NSl+1 

aYI.’ISl^Nj-SIGYiMSljfO’-CAUSINN^D-AUlhSl/N^ 1) ) 

0X(rlS»N)-0. 

0Y(rli>fJ)-0. 

CONTINUE 

AM(2»2>I jNXF'.A)«Af1(2>2»IJNXh4)-0X(2»2iNXf1A>/XXI (2/2 JNXM4) 
aN( 2»2>I }NXMA)-AN(2>2/IiNXF!4 )-QY( 2/ 2 J NXMA ) / Y ET A ( 2/ 2) NXH4 ) 

CCNTlNUF 

RETURN 

END 

SU3R0UTINE VEC2 

C0f1,10N/Fl/NXM,NXMl»NX«2>NXri3/NXM4/Nl>Nl/Nll/ hll/NSl#NS2,NS# NS/rtSl 
1 >NC(iL»NCWL1# NCi<L2^NChLM>NCWLPiNSTRUT 
C0fl!iUN/F2/RD(61» 55)/U(61»55)>V(6l#55)>P(61^55)>T(61»55)» 

1 SH(61,55) >H(61, 55)»RJS(55)#US (55)>V5(55)> 

2 Ps(55)/TS(55),SHS(55)>HS(55) 
CDNrtCN/FA/iIGX(61,55),SlGY(61>55)/TAUXY(61/55)»0X(61,55)» 

1 QY(61»55>>SIGXS(55)»SIOYS(55),TAUXYi(55)#QXS(55)* QYS(55) 
CONMaN/F5/AU(61,55><»)»AM(61,55>A),AN(61>55»A)>AUl(61»55»A) 
C0r1fluN/F7/X( 61»55)# Y(bl,55)/XXI (61/55)/YXI (61»55)» 

1 X£TA(61,55),Y£TA(61,55),AJ(61,55) 
COrliaN/F13/6AMA,RGAS»CP^PR»PRT>CVIS> FM/ PF» T F »U F> RO F» SHF» HF 
C0Mii0N/FlA/DAU(61»55» A)» dOAU (55><f )>BDAUS( 55» A)> T1 ( 55>5)> 
lCClMPY>GMl>GM2»GM3>GN<»»Gn5 
R0(2# 2JNXM4) »AUl(2f 2»1|NXMA) /AJ(2/2JNXMA) 

U(2>2JNXM4) -AU1(2»2»2JNXM<» ) /AUl ( 2^ 2 # U NXM4 ) 
V(2#2;NXM4)-AU1(2>2»3;NXM4)/AU1(2>2 a 1JNXM4) 
0X(2>2jNXM4)-(J(2>2JNXM4)*U(2»2iNXM4)+V<2/2;NXH4»*V(2/2iNXM4))/2. 
SH(2»2jNXM4)»(AUl(2>2» 4jNXf14)/AUl (2> 2»1JNXM4 )-QX ( 2> 2} NXM4 ) ) t-GAMA 
HI 2> 2JNXH4 )-SH(2»2iNXM4) + QX( 2/2JNXH4 ) 

P(2,2iNXH4)-GM2*R3(2,2 JNXH4)+i>H(2/2l NXh4) 

T(2^2JNXM4 )-SH(2»2jNXM4)/CP 

RETURN 

END 

SUBRUUTINE BOUND 

CUrtNUN/Fl/NXM, NXtU»NXrt2/NXM3/NXft4,Nl^m#Nll> Ml 1# NS 1^ NS2, NS/ MS# Ni 1 
1»NCNL/NCWL1>NCWL2»NCWLM/NCWLP/N5TRUT 
CtJl1N0N/FAl/M12,'113,N12/N13/NXf15/NXM6/MS2»MS3/MS4/rtSPl»NSP2»MSP3 
l/NSrll/NSlPl/NS2Ml»NSlrtl»NS2Pl/NNS/NSlM2/NXM7 
COnNON/F2/RO(61»55 )»U(61/55) /V(bl/55),P(6i,55)/T(61»55)> 

1 SH(61»55)/H(61»55)/R0S(55)/US(55)/VS(55)> 

2 PS(55)»TS(55)/SHS<55)/HS(55) 

C0rtNUN/F3/VISL<61» 55)/VlST(61/55)/VIS(61#55)/VISS(55) 
C0MN0N/F7/X(61»55)#Y(61»55)/XXI (61/55)/YXI (61/55)/ 

1 X£TA(61,55)>YETA(6l/55)/Aj(61/55) 

C0flf10N/F8/XS(55)/ YS(55 )/ XSXI { 55 ) / YSXl ( 55 ) / XS ET A ( 55 )/ 

1 YSETa(55)/AJS(55) 

C0MMUN/Fg/INT(55)/TE1(55)/TE2(55)/TE3(55)/T£4(55) 
C0MM0N/F10/0T(61/55)/ JTM 

COMrtON/F12/L SYM/NFLOW/ LFI /LFO/ILI/LMAX/LW/ LEXIT/CRIT/FDT/FOTL/ 

1 TIME/CCP/CCT/L/LK/TIMEL/CRITAVG 

C0Mf10N/F13/GAMA,RGAS/0P/PR/PRT/CVlS/FM/PF/TF/UF/R0F/SHF/HF 
TEl(25N12)»Q8VGATHR{T(Mll/2»NXM3)/INT(ljN12) )TEl(2jN12) ) 
T£2(2iN12)-QeVGATHk(P(hll/2;NXM3)/INTaiN12) ;TE 2 ( 2 lN 12 ) ) 

TE3(2)N12 )«TE2(2;N12)/RGAi/TEl(2)N12 ) 

TE4(2JN12)-CP^TE1(2JN12) 

T( i11/2;NXM3 )«3eVSCATR(TEl(2;N12) / INT (1;N12) i T ( M 1 / 2;NXM3 ) ) 
P(11/2JNXM3)»Q8VSCATR(TE2(2/N12)/1NT (1/N12); F(M1/2|NXM3) ) 
R0(fll/2JNXM3 )-38VSCATR{TE3(2jN12) /INT(liN12) JR0(H1/2|NXM3)) 
SH(H1/2;NXK3 )«Q8VSCATR (TE 4( 2 ; N12) / INT (1 ; N12 ) iSH ( Hl/2 JNXM3 ) ) 
TEl(2jN12)-0. 


DAMPC 

46 

DAMPC 

49 

DAfIPC 

50 

DAMPC 

51 

DAMPC 

52 

DAMPC 

53 

DAMPC 

54 

DAMPC 

55 

DAMPC 

56 

DAMPC 

57 

DAMPC 

58 

DAMPC 

59 

DAMPC 

60 

DAMPC 

61 

DAMPC 

62 

DAMPC 

63 

DAMPC 

64 

VEC2 

1 

VEC2 

2 

VEC2 

3 

VEC2 

4 

VEC2 

5 

VEC2 

6 

VEC2 

7 

VEC2 

6 

VEC2 

9 

VEC2 

10 

VEC2 

11 

VEC2 

12 

VEC2 

13 

VEC2 

14 

VEC2 

15 

VEC2 

16 

VEC2 

17 

VEC2 

18 

VEC2 

19 

VEC2 

20 

VEC2 

21 

VEC2 

22 

VEC2 

23 

VEC2 

24 

BOUND 

1 

BOUND 

2 

BOUND 

3 

BOUND 

4 

BOUND 

5 

BOUND 

6 

BOUND 

7 

BOUND 

8 

BOUND 

9 

BOUND 

10 

BOUND 

11 

BOUND 

12 

BOUND 

13 

BOUND 

14 

BOUND 

15 

BOUND 

16 

BOUND 

17 

BOUND 

18 

BOUND 

19 

BOUND 

20 

BOUND 

21 

BOUND 

22 

BOUND 

23 

BOUND 

24 

BOUND 

25 

BOUND 

26 

BOUND 

27 
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dOUND 

U(^U^^5NXM3)«08\/SCATRtTEl(2;N12J^ INT(liN12);U(f11>2;NXM3) ) B3UN0 

V(Mi/2jNXM3 )«Q3 VSCaTR(TE 2(2;N12)^ INT (1;N12) j V(Ml/2 jNXM3) ) SOUND 

TE3(2;N12)-TEA(2;N12)MTEl(2|Nl2)*Te:li2)N12)+TE2(2 jN12)^TE2(2iN12)BQUND 
l)/2. BOUND 

Htil>2;NXf13)-08\/S:ATR(TE3(2;N12)>INT(ljN12) JH(hli2iNXM3) ) BOUND 

IF (NSTRUT.EQ.O)Ga TO 15 BOUND 

DO 10 N*NS1#NS2 BOUND 

Ul?1S>N)-0* BOUND 

V(.iS#N)-0. BOUND 

P(rtSiN)»P(MSPl,M) BOUND 

T(fli>N)«T(MSPl>N) BOUND 

SH(rtS^N)-CP+T(MS>N) BOUND 

R0(f1S,N)-P (NS# N)/RGAS/T(NS>N) BOUND 

10 H(MS#N)-SH(MSP1#N) BOUND 

DO 20 N»1#NS BOUND 

NN-N^NSINI BOUND 

US(N)-0. BOUND 

VS(N)-0. BOUND 

PS(N)*P(MS1#NN) BOUND 

rS(N)-T(MSl#NN) BOUND 

SHS(N)*CP+TS (N) BOUND 

ROS (N)«PS (S ) /RGAS/TS(N) BOUND 

20 HS(N)*SHS(N) BOUND 

15 CONTINUE BOUND 

DO 25 N*2#N11 BOUND 

T(1#N)-T(2#N ) BOUND 

P(1#N)«P(2#N) BOUND 

SH(I#N)-CP*T(1#N) BOUND 

RD(1#N)-RD(2#N) BOUND 

UU#N)»U(2^N ) BOUND 

V(1>N)*VC2#N) BOUND 

25 H(1#N)-H(2#N) BOUND 

IF(LSYN.EQ.O)GO TO 55 BOUND 

DO 60 N*2#N11 BOUND 

V(1#N)*0. BOUND 

60 H( 1#N)-SH( 1#N)+U(1#N)*UU#N) /2. BOUND 

55 CONTINUE BOUND 

C COWL PLATE OR CENTER STRUT oTAkTS FROM N-NCWLl AND ENOS AT NCWL2. BOUND 

IF(NCWL.EQ.O)Ga TO 35 BOUND 

DO 30 N-NCkLl# NCWL2 BOUND 

H(1#N)»SH{1#N) BOUND 

U(i#N)«0. BOUND 

30 V(1#N)*C. BOUND 

35 CONTINUE BOUND 

IF(LEXIT.EQ*0)G3 TO 36 BOUND 

U(i#Nl;Ml )-2.*U( l#Nll;hlJ -U(1#M2;M1 ) BOUND 

V(l#Ni)Ml)-2.’^\/(l#Nll;Ml)-V(l#ra2iMl ) BOUND 

p(i#NUin)-2.*‘ p(1#nii;md- p(i#ni2;mi) bound 

SH(1#N1;M1 )-2.*SH( 1#N11;M1)-SH(1#N12;M1) BOUND 

GO TO 37 BOUND 

36 CONTINUE SOUND 

u(i#numi)-u(i#nii;mi) bound 

V(i#Ni;Ml).V{l,Nll;hl) BOUND 

P(1#M|M1)-F(1#N11JM1) BOUND 

SH(l#NljMl)-SH(l#Nil;Ml) BOUND 

37 CONTINUE BOUND 

T(i#NijMl)«SH(l#Ni;Ml)/CP SOUND 

R0( 1#N11M1 )«P( 1#N1;M1) /(RGAS^T (1# Nli Ml) ) BOUND 

H( i#NlJMl)-5H(l#Nl5Ml) + ( J(l#NliMl)*U(l#Nl;Ml ) + V(l#Nl;Ml) ♦ BOUND 

1V(I#M;M1) )/Z. BOUND 

1I-Q8SLT(T(1#1;NXM)#C.5) BOUND 

IF(II.EQ.NXM)GO TO 50 BOUND 

i^FXTE(6# 500)L BOUND 

500 F0RMAT(/#2X# 'NEGATIVE TEMPERATURE IN THE FIEL0S15X# BOUND 

1 'NO. OF ITERATIONS*'# 15) BOUND 

DO 40 N.ljNl BOUND 

WRI r£(6f 550)N#X(1#N) BOUND 
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00 40 

BJUND 

96 



aOUNO 

97 


WPIT£(6, 560)Y(d»M)#U(.’1»N)»V(M/N),P(i1»N)»T(M>N)>H(rt>M)>VIST(M/N) 

30UN0 

98 


CONTINUE 

SOUND 

99 


IRNSTRUT.EQ.0)G3 TO 45 

BOUND 

100 


WRITE(6,570) 

SOUND 

101 


mRITE( 6» 560) (YS(M)»Ui(N)»VS(N)>PS(N)/TJ,(N)>HS(N)RISS(N)>N-l»NS) 

SOUND 

102 


CONTINUE 

SOUND 

103 


LK«2 

BOUND 

104 


RETURN 

SOUND 

105 

50 

CONTINUE 

BOUND 

106 

550 

F0RrtAT(//8X, 'Y'» 10A» 'U'»10X> 'WllX/ ‘P'/llX# 'TS12X#' H'#10X/ 

BOUND 

107 


I'VIiTSlOXR BOOT STATI0N«'»I3>5X» 'X- SF10.5) 

BOUND 

108 

560 

FORrlAT(2X,F12.5>F12.5,F12.5»F13.4>F12.5»2F15.5) 

BOUND 

109 

570 

F0Kf1AT(/>8X» *YS'/PX> 'USS9X/ •VS*/1CX#‘PS'/10XRTS'/12X#'HS'» 

BOUND 

110 


110X» 'VISS'f/) 

BOUND 

111 


IJNXM )-VSaRT(T(l»liNXM)iVISL( 1»1JNXM) ) 

BOUND 

112 


VIiL(l,i;NXH)-C\/IS*VlSL(lR;NXh)*T(l>l{NXM)/ (T ( 1 iNXM ) ♦ 110 . ) 

BOUND 

113 


IKNiTRUT.EQ.O)RETURN 

BOUND 

114 


ViS3(lJNS)»VSQRT(Ti>(lJNS.) JVISS (1;NS) ) 

BOUND 

115 


vissdjNS )-c vis*vis5(i;N3)*rs(i;NS) / (TidiNS )+iio. ) 

BOUND 

116 


return 

BOUND 

117 


END 

BOUND 

110 


3U3RQUTINE IBOUNO 

IBOUND 

1 


CtjrtrtCiN/Fl/NXMRXMl,NX?12,NXh3>NXM4»Nl>m>Nll> hi 1 » NS 1> NS2» NS* MS# MS 1 

IBQUND 

2 


l#NC(a#NCWLl#NCWL2»NCRLM#NCwLP#NSTRUT 

IBOUND 

3 


CD1rtuN/FAl/M12*M13»N12»M3*NXM5»NXM6#MS2»MS3*hS4*MSPl»MSP2#MSP3 

IBOUND 

4 


l#NSf1i#NSlPl# NS2M1»NS1I11#NS2P1>NNS>NS1M2#NXM7 

IBOUND 

5 


CUf1M0.'t/F2/RD(61» 55 )»U(6l»55)#V(61>55)#P(61#55)/T (61/55) » 

IBOUND 

6 


1 SrU61»55),H(61#55),R0S{55)#US(55)#VS(55)# 

IBOUND 

7 


2 PS(55)»TS(55)#SH3(55)»HS(55) 

IBOUND 

8 


C0n.10h/F4/SIGX(61*55)»SIGY(61#55)#TA0XY(61#55)#QX(61»55)# 

IBOUND 

9 


1 Qr(61*55),SIGXS(55),iI6YS(55)#TAUXYS(55)#QXS(55)#UYS(55) 

IBOUND 

10 


COMMON/FS/ AU(61» 55# 4 ) , AM ( 61# 55# 4) » AN ( 61, 5 5 » 4 ) # AU 1 ( 61* 55# 4 ) 

IBOUND 

11 


C0,1rt0N/F6/AUS(55#4)#AMS(55#4)#ANS(55#4)#AUSl (55#4) 

IBOUND 

12 


C0flnGN/F7/X(61*55 )*Y(61*55)#XX1(61*55)*YXI(6 1*55)* 

IBOUND 

13 


1 XETA(61* 55)#YETA(61#55)#AJ (61#55) 

IBOUND 

14 


CuMMON/FB/XS (55) * YS<55)#XiXl ( 55 )# YSX I (55 ) # XS ETA(55 )# 

IBOUNO 

15 


1 YjETA(55)*AJS(55) 

IBOUND 

16 


COrt.lUN/FlE/LSYM/MFLOW# LFI# LF 0# I LT# LM AX* LW# LE XIT# CR I T# FOT#FOTL* 

IBOUND 

17 


1 TIHE*CCP*CCT*L* LK*TIMEL*CRITAVG 

IBOUND 

18 


C0MMDN/F13/GAHA* RGAS* C P# P R# P RT# CV IS# FM# PF#TF#UF# ROF#SHF#HF 

IBOUND 

19 


OL) 5 N«2#N11 

IBOUND 

20 


R0(1*N)-RQ(2#N) 

IBOUND 

21 


Pd*N)- P(2*N) 

IBOUND 

22 


7(1*N)- T(2#N) 

IBOUND 

23 


Hd*N)» H(2#N) 

IBQUND 

24 


3H(1#N)»SH(2*N) 

IBQUND 

25 


CDNl»YXId#N )/XXI d#N) 

IBOUND 

26 


CON2-1.+CON1*COM1 

IBOUND 

27 


Ud*N)«(H(l*N)-SH(l#N) )*2./C0N2 

IBOUNO 

26 


Ud#N)>SORT(U(l*N) ) 

IBQUND 

29 

5 

Vd#N)«CDNl*U( 1*N) 

IBOUNO 

30 


00 20 N-2#N11 

IBQUND 

31 


RD(M1*N)>R0(M11*N) 

IBOUND 

32 


SH(«1*N)«SH(M11#N) 

IBOUND 

33 


P(M1*N)« P(M11*N) 

IBOUND 

34 


T(M1#N)» T(M11#N) 

IBOUNO 

35 


H(M1*N)- H(M11*N) 

IBOUNO 

36 


CONl-YXI (M1*N)/XXI (M1*N) 

IBOUND 

37 


CON2-1.+CON1*CQN1 

IBOUND 

38 


U(Mi#N)-(H(Ml#N)-SH(Mi#N))*2./CQN2 

IBOUNO 

39 


U(M1*N)«SQRT(U(M1*N)) 

IBOUND 

40 

20 

V(rtl*N)»CONl*U(Ml#N) 

IBOUND 

41 


IFUSYM.EO.DGO TO 60 

IBOUND 

42 


IF(NCRL.EQ.O)GJ TO 62 

IBOUND 

43 


IF(NCWL1.LE.2.ANO.NCWL2.6E.N11)GO TO 60 

IBOUND 

44 


00 50 N«1*NCWLM 

IBOUNO 

45 
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JU>N)«U(2>N) 

I30UN0 

46 

50 

V(1»N)«V(2»M) 

I BOUND 

47 


IF(?<CWL 2 .GE.Nll)Ga TO 60 

I BOUND 

48 


00 55 N-NCWLP>N11 

IBOUND 

49 


U(1>N)>U(2>N) 

IBOUNO 

50 

55 

V(1>N)-V(2»N) 

IBOUND 

51 


60 TO 60 

IBOUND 

52 

62 

Continue 

IBOUND 

53 


00 56 N-2»N11 

IBOUND 

54 


U(1/N)>U(2»N) 

IBOUND 

55 

56 

V(i»N)-V(2»N) 

IBOUND 

56 

60 

continue 

IBOUND 

57 


IF(NSTRUT.EQ.0)63 TO 6 

IBOUND 

56 


00 10 N»1>NS 

IBOUNO 

59 


NN«N+NS1-1 

IBOUNO 

60 


R0(MS>NN)-R0(NSP1»NN) 

IBOUNO 

61 


SH(NS>NK)-SH(MSP1>NN) 

IBOUNO 

62 


P(MS>NN)- P(MSP1»NN) 

IBOUND 

63 


T(Mi>NN). T(MSP1»NN) 

IBOUND 

64 


H(rtS,NN)- H(MSP1»KN) 

IBOUND 

65 


CONl-TXl (f1S»NN)/XXI(F!i^NN» 

IBOUND 

66 


CCN2-1.+C0N1*C0N1 

IBOUND 

67 


U(NS>NN)-(H(MS»NN)-SH(hS>NN) )*2./CDN2 

IBOUNO 

68 


0(MS»NN)»SQRT( J(MS»NN) ) 

IBOUNO 

69 


V(NS»NN)«C0N1*U(MS»NN) 

IBOUNO 

70 


RDS(N)«R0(MS1> NN) 

IBOUNO 

71 


PS(NJ- P(MS1»NN) 

IBOUND 

72 


T5(N)» T(MS1»NN) 

IBOUND 

73 


HS(N)- H(MS1,NN) 

IBOUNO 

74 


SHS(N)«SH(MS1>NN) 

IBOUNO 

75 


C0N1«YSXI(N)/XSXI(N) 

IBOUNO 

76 


C0N2-1.+C0N1+C0N1 

IBOUND 

77 


US(N)-(HS(N)-SHS(N))*2./C0N2 

IBOUND 

78 


US(N)-SORT(US(N) ) 

IBOUND 

79 

10 

VS(N)-C0N1*US(N) 

IBOUNO 

80 

6 

continue 

IBOUNO 

81 


IF (LExIT.EQ.O) 63 TO 36 

IBOUNO 

82 


UU>Ni;MH«2.*U{l»NllJMl)-U( l>N12jhU 

IBOUNO 

63 


V<1>MJM1)«2.*V(1»N11;m 1)-V( l>N12;r.l ) 

IBOUNO 

84 


Sh(l,NljMl)«2.*SH(l,NllJMl)-SH{l»N12iMl) 

IBOUNO 

85 


P(i#Nl}rtl)-2.* P(1»N11JM1)- P(i»N12jMl) 

IBOUND 

86 


GO TO 37 

IBOUNO 

87 

36 

CONTINUE 

IBOUND 

86 


U(l»NlIhl)«U(l»NlUni) 

IBOUNO 

89 


V(l/Nijm»«V(l»Nll}f11) 

IBOUND 

90 


jH{i»NljMll-SH{l»NllJfHI 

IBOUNO 

91 


P( i,NUKl)«P(l,NH}hl) 

IBOUNO 

92 

37 

CONTINUE 

IBOUNO 

93 


T(1/N1}N1)-SH(1,N1}N1)/CP 

IBOUNO 

94 


RO(l»Nl;Ml)»P(l»NlJrtl)/(RGAi*T(l»MjNin 

IBOUNO 

95 


H(l»NljMl)«SH(l,Ni;Ml» + (U(l,NIiMl)*U(l/Nl;t11) + V(l»NUrtU* 

IBOUNO 

96 


lV(I»NijMl) )/2. 

IBOUNO 

97 


Il-08SLT(T(l,l;NXt1)#0.5) 

IBOUNO 

98 


IF(II.EQ.NXM)RETURN 

IBOUNO 

99 


WRITE(6,500)L 

IBOUND 

100 

500 

F0RMAT(/»2X» 'NEGATIVE TEMPERATURE IN THE FIELD'#15X» 

IBOUND 

101 


1 'NO. OF ITERATIONS -*/I5) 

IBOUND 

102 


DO 40 N«1,N1 

IBOUND 

103 


RPITE<6»550)N,X(I.N) 

IBOUNO 

104 


00 40 MM-1,M1 

IBOUNO 

105 


M«rtI-MM+l 

IBOUND 

106 


WRITE (6, 560) Y(M»N)»U(M»N»/V(M>N),P(M#N)»T(M#N)»H(.*1»N» 

IBOUNO 

107 


CONTINUE 

IBOUNO 

108 


IF(NSTRUT.EQ.0)G0 TO 45 

IBOUNO 

109 


WRITE(6»570) 

IBOUNO 

110 


WRITE (6. 560) (YS(N)»US(N)» VS ( N ) / PS (N ) » TS ( N ) » H S< N ) , N-1* NS ) 

IBOUND 

111 

45 

CONTINUE 

IBOUNO 

il2 

550 

F0RMAT(/#ex»'Y'»10Xf 'U'»10X, ' v • » 1 1X» ' P ' » 11X» 'T ' » 12 X> • H*,10X» 

IBOUNO 

113 
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STATION* 13/ 5X# 'X» S Flo. 5) I BOUND 

560 F0RnAT<2X/F12.5/F12.5/F12.5/F13.4/F12.5/F15. 5) I BOUND 

5 70 FCjR.iAT(//8X/ «YS'/9X/ •US'/9X/ • VS MOX / ' PS • / 10 X/ ' T5 ' / 1 2X/ ' HS • / / ) I BOUND 
LK-2 IBUUNO 

ftETJRN IBQUND 

END IBOUND 

SUdKGUTINE VORT VQRT 

CDnMJN/Fl/NXM/NXMl/NXn2/NXN3/NXh^/Nl/hl/fUl/ M 1/ N5 1/ NS 2/ NS / MS/NSl VORT 
l/NCwL/NCWLl/NCWL2/ NC it L 11/ NCk L P/ NS T RUT VQRT 

CGN.nDN/FM/N12/M13/N12/Nl3/NXN5/NXh6/f1S2/;1S3/hSA/NSPl/MSP2/MSP3 VORT 
l/Nif11/NSlPl/NS2Ml/NSlrU,NS2Pl/NNj/NSlM2/NXh7 VQRT 

CUNM0N/F2/R0(6I/ 55 ) / U ( 61/ 55 ) / V ( 61 / 55 ) /P ( 61 / 5 5 ) / T ( 6 1/ 55 ) / VQRT 

1 SH (61/55)/H (61/ 55) /RGS (55 )/US (55 )/ VS (55)/ VORT 

2 Pi(55)/T5 (55)/SHS(55)/HS(55) VORT 

C0f1r1QN/FA/SIGX(61/55)/SlGY(6 1/55)/TAUXY(61/5 5)/ 0X( 61/ 55)/ VORT 

1 UY(61/55)/SIGXS(55)/SlGYS(55)/TAUXYS(55)/CXS(55)/UYS(55) VORT 

C(JrtNUN/F7/X(61,55)/Y(61,55)/XXI(ol/55)/YXl(6 1/55)/ VORT 

1 XETa(61/55)/YETA(61/55)/ AJ(61/55) VQRT 

CUM.iUN/F8/XS (55)/ YS(5:>)/XSXX (55)/YSXI (55)/XStTA(55)/ VORT 

1 YSETA(55)/AJS(55) VORT 

SlGX(2/l;NXM2)-(U( 3/l;NXrt2)-U(l/i;NXM2))/2. VORT 

^1GY(2/1;NXM2)-( V(3/i;NXrt2)-V(l/l/NXM2))/2. VQRT 

DU 10 N-l/Nl VORT 

S1GX(1/N)-U(2/N)-U(1/N) VQRT 

SlGYd/N) .V(2/N)-V(l/N) VORT 

SlGY(i/N)»V(2/N)-V(l/N) VQRT 

SlGX(hl/N)«U(Ml/N)“U(Mll/N) VORT 

10 SIGY(hl/N)-V(Ml/N)«V(Hll/N) VORT 

UX(1/2;NXM3) *(U(1/3;NXM3)-U(1/1JNXM3 ) )/2. VORT 

QY(i/2jNXM3)-( V(1/3|NXM3)-V( 1/ liNXM3 ) )/2. VORT 

DO 20 M*l/11 VORT 

MX(f1/l)«U(iM,2)-U(M/ 1) VQRT 

QY(H/l)*V(M/2.)-V(M/l) VQRT 

OX(.i/Nl)*U(P/Nl)-U(M/Nll) VORT 

20 QY(rt/Nl)*V(M/Nl )-V(M/Nll) VQRT 

IF(NSTRUT.EQ.O)GQ TO 50 VORT 

DO 30 N-NS1/NS2 VORT 

NN-N-NSl+1 VQRT 

SIGX(MS/N)*U(MSP1/N)-J(MS/N) VORT 

SlGY(MS/N)*V(MSPl/N)-V(rtS/N) VORT 

S1GX(MS1/N)- (US(NS)-U(MS2/N) )/2. VORT 

30 SIGY(MS1/N)» (VS(NN)-V(MS2/N) )/2- VORT 

DO 40 N»2/NSM1 VORT 

40 QXS(N)«(US(N+1)-US(N-1) )/2. VORT 

QXS(l) -(US (2 )-U(nS/NSlNl) ) /2. VQRT 

QXS(NS)-(U(MS/NS2Pl)-JS(NSMl))/2. VQRT 

DU 60 S-l/NS VORT 

NN-N-l+NSl VQRT 

60 S1GXS(N)«US(N)-U(MS1/NN) VORT 

QXS(i;NS)-(XSXI(UNS)*SIGXS( l;NS)“XS£TA(liNS )♦ QX S ( U NS ) ) / A J S ( 1 ; VQRT 
INS) VQRT 

50 CUNTINJE VQRT 

QX(1/1)NXM)-(XXI (l/liNXM)*SIGX(l/ 1 jNX M )-X £ T A ( 1 / 1 jNXM ) ♦ OX(l/li VORT 

INXM) )/AJ(l/lrSXM) VORT 

QYd/ljNXM) - (YETA(1/1;NXM)+ Q Y ( 1 / 1 ; NxM ) -YX I (l/l;NXM)^SIGY( l/l; VORT 

INXN) )/AJ(1/1;NXM) VORT 

TAUXY(l/l|NXM)-QX(l/l;NXM)-QY(l/l|NXM) VORT 

TAJXY(l/l;NXM)-VA8S(TAUXY(l/ 1 ; NXM ) i T AUXY ( 1/ 1 j nXM ) ) VQRT 

RETURN VORT 

END VORT 

SUBROUTINE EODY( NNX 1/ NNX2/ N Y 1/ NY2 / Nirf ALL/ NS TR ) EDDY 

C ELDY VISCOSITY MODEL TAKEN FROM AIAA PAPER NO. 78-257 BY BALDWIN EDDY 

C AND LOMAX. NWALL ■ 1 FOR ONE WALL AND ■ 2 FOR TWO WALLS. EDDY 

DIMENSION YD(61)/VISTI (61)/ VIST0(61)/FY(61) / UDIFK61) EDDY 

CUrtMDN/Fl/NXM/NXMl/NXM2/NXM3/NXM4/Nl/Mi/Nli/ Mi 1/ NS 1/ NS2/ NS/ MS/ MSI EDDY 
1/NCWL/NCWLl/ NCWL2/NCWLM/NCWLP/NSTRUT EDDY 

C0MM0N/F2/R0(61/ 55 ) / U ( 61 / 55 ) / V ( 61 / 55 ) / P ( 61 / 5 5 ) / T ( 6 1/ 55 ) / EDDY 

1 SH(61/55)/H(61/55)/R0S(55)/US (55 )/VS(55)/ EDDY 
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2 PS(55),TS(55),SHS(55)^HS(55) 



EDDY 

9 


CCJ.i?10N/F3/VISL (61,55), VIST(61#55)fVIS(6i^55) 

>VISS(55) 


EDDY 

10 


CUf1rt0N/F4/$IGX(61,55),SIGY(61,55),TAUXY(61,5 5),QX(61,55), 


EDDY 

11 


i GY(61,55),SIGXS(55),SIG»'5( 55),TAi;XYS(55),0XS(55),0YS(55) 


EDDY 

12 


CDMrt0N/F7/X(61,55),Y(6l,55),XXl(61,55),YXI(61,55), 


EDDY 

13 


1 XETA(61,55),YETA(6i,55), AJ(61,55) 



EDDY 

14 


C0rt.1jN/F8/XS(55)» YS(55 )»XiXi ( 5 5 ) » YSXI ( 55 ) # XS ET A ( 55 »» 


EDDY 

15 


1 Yi£TA(55),AJS(55J 



EDDY 

16 


real k»kk 



EDDY 

17 


DATA AP»CCP»CKL£8»CWK»K>KK/26.>1.6»0.3»0.25»0.4» .0168/ 


EDDY 

18 


KSTR EQUAL TO 1 FOR THE REGION WHERE LOWER 

SURFACE GF THE 

STRUT 

EDDY 

19 


FL'Rrii THE UPPER WALL OF THE DOMAIN IN WHICH 

EDDY VISCOSITY 

IS 

EDDY 

20 


BEING CALCULATED. IT IS ZERO OTHERWISE. 



EDDY 

21 


nNY-NY2-NY1+1 



EDDY 

22 


NNYl-NYl 



EDDY 

23 


NNY2-NY2 



EDDY 

24 


IF (NWALL.EQ.DGO TO 5 



EDDY 

25 


NNY2-NY2 



EDDY 

26 


NNYI-NYI+NSY/2 



EDDY 

27 


continue 



EDDY 

28 


NNY2M1«NNY2-1 



EDDY 

29 


NMYlPl-NNYl+1 



EDDY 

30 


NNY-NNY2-NNY1+1 



EDDY 

31 


DO 100 N-NNX1.NNX2 



EDDY 

32 


iKNSTR.EQ.DGG TO 10 



EDDY 

33 


TAOW«ABS( VIS(NNY2.N)*QX(NNY2,N) )*R0(NNY2»N) 



EDDY 

34 


YPI*SQRT(TAUH)/VIS(NNY2#N) 



EDDY 

35 


OU 20 M-NNY1.NNY2 



EDDY 

36 

20 

YL(r1)»Y{NNY2»N)-Y(M»N) 



EDDY 

37 


GO TO 15 



EDDY 

38 

10 

NN-N-NSl+1 



EDDY 

39 


TAJW"ABS(VISS(NN)AQXS(NN) )*ROS (nN) 



EDDY 

40 


YPi»SQRT(TAUW)/WISS (NN) 



EDDY 

41 


DO 30 M«NNY1»NNY2 



EDDY 

42 

30 

Y0(.D«YS(NN)-Y((1>N) 



EDDY 

43 

15 

CONTINUE 



EDDY 

44 


DO 40 M-NNY1.NNY2 



EDDY 

45 


YP-YPI + YDH) 



EDDY 

46 


RL— YP/AP 



EDDY 

47 


RL-EXP(RL) 



EDDY 

46 


RL«YD(M)*(1.-RL) 



EDDY 

49 


V15TI (M)-RO(M»N)*K*KWRL*RL*TAUXY(M»N) 



EDDY 

50 


UDIF1(M)-U(M»N)+U(M»N) +V(M.N )*V(M.N) 



EDDY 

51 


UDIF1(M)-SQRT(UDIF1(M) ) 



EDDY 

52 

40 

FY(rt)-RL+TAUXY(M.N) 



EDDY 

53 





EDDY 

54 


SOMETIMES CALCULATIONS OF FYMAX MAY REQUIRE 

SOME attention 

DUE 

EDDY 

55 


TO MORE THAN ONE MAXIMA OCCURING IN FUNCTION 

FY. 


EDDY 

56 





EDDY 

57 


FY.1AX-FY (NNY2M1) 



EDDY 

58 


Y.1AX«Y0(NNY2M1) 



EDDY 

59 


DO 50 MM-NNY1P1.NNY2M1 



EDDY 

60 


M-NNY2-MM+NNY1 



EDDY 

61 


IF(FY{M-1).GT.FY(M))GJ TO 50 



EDDY 

62 


FYMAX-FY(M) 



EDDY 

63 


YMAX-YO ( M) 



EDDY 

64 


GO TO 51 



EDDY 

65 

50 

CONTINUE 



EDDY 

66 


FYflAX-FY(NNYlPl) 



EDDY 

67 


YMAX-YO(NNYIPI) 



EDDY 

66 

51 

CONTINUE 



EDDY 

69 


ULIF-QeSMAX(UOIFl(NNYIJNNY) ) 



EDDY 

70 


FWAKEl-YHAX+FYMAX 



EDDY 

71 


FWAKE2«CWK*YMAX*U0IF*J OIF /FYMAX 



EDDY 

72 


FWAKE«AMIN1( FWAKE1.FWAKE2) 



EDDY 

73 


DO 60 M«NNY1»NNY2 



EDDY 

74 


FKuEB-l./ (l.+5.5't'(CKLEa*Y0(M )/YMAX)**6) 



EDDY 

75 

60 

VISTQ(M)-KK*CCP*RO(M»N)*FWAKE*FKL£B 



EDDY 

76 
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70 

100 


210 


220 

221 

230 


2^0 

200 

500 


IGJT-O 

EDDY 

77 

DU 70 MM»NNY1»NNY2 

EDDY 

7b 

M-.>iNY2-hM + NNYl 

EDDY 

79 

DVIST*VISTO(M)-VISTI (M) 

EDDY 

80 

Ih (OVIST. LE.O* )I0JT-1 

EDDY 

81 

VI;>T(MfN)«VISTr(M) 

EDDY 

82 

IF( iuUT.EQ.l )VIST(h>N) -VISTG(M) 

EDDY 

83 

CONTINUE 

EDDY 

8A 

CONTINUE 

EDDY 

85 

IE(N«ALL.E0.1)G0 TU 500 

EDDY 

8b 

NKY2-NNY1-1 

EDDY 

87 

NNY1*NY1 

EDDY 

88 

NNr*NNY2-NNYl+l 

EDDY 

89 

NNYlPl-NNYl+1 

EDDY 

90 

NNY1P2«NNY1+2 

EDDY 

91 

DD 200 N«NNX1,NNX2 

EDDY 

92 

TAUirf-ABS(VIS (NNY1>N) + JX(NNY1/N) )*P0( NNU^N) 

EDDY 

93 

YPI«S0RT(TAUw)/VIS(NNY1^N) 

EDDY 

94 

DO 210 M-NNYl^NNY2 

EDDY 

95 

YD(N)«Y(M,N)-Y(NNY1^N) 

EDDY 

96 

YP»YPI+YO(M) 

EDDY 

97 

KL--YP/AP 

EDDY 

9b 

RL-EXP(PL) 

EDDY 

99 

RL-Y0(M)*(1.-RL) 

EDDY 

100 

VISTI (M )«R0( M^N) *K*K+RL*KL*TAUXY(h/N ) 

EDDY 

101 

UDIFl(N)«U(M/N><‘U(n^N)+V(M>N)^V(ri^N) 

EDDY 

102 

UD1F1{M)-S0RT(UDIF1(M) ) 

EDDY 

103 

FY(ii)«RL*TAUXY(M#N) 

EDDY 

104 


EDDY 

105 

S0N6TIKFS CALCULATIONS OF FYhAX MAY REQUIRE SOME ATTENTION DUE 

EDDY 

106 

TO MORE THAN ONE MAXIMA OCCURING IN FUNCTION FY. 

EDDY 

107 


EDDY 

108 

FYMAX-FYCNNYIPI) 

EDDY 

109 

YMAX«YD(NNY1P1) 

EDDY 

110 

DO 220 M«NNY1P2> NNY2 

EDDY 

111 

IF(FY(M).GT.FY(M-1))GU TO 220 

EDDY 

112 

FYMaX«FY(M-1) 

EDDY 

113 

YMAX-YO(M-l) 

EDDY 

114 

GO TO 221 

EDDY 

115 

CONTINUE 

EDDY 

116 

F YMAX«FY(NNY2) 

EDDY 

117 

YMaX-YD (NNY2) 

EDDY 

118 

CONTINUE 

EDDY 

119 

UDIF«QBSMAX(UDIF1(NNY1;NNY) ) 

EDDY 

120 

FWAKEl-YMAX^FYMAX 

EDDY 

121 

FWAKE2*CWK + YMAXtJ0IF*'JOlF /F YMaX 

EDDY 

122 

FWAKE-AMINK FWA<E1^FWAKE2) 

EDDY 

123 

DO 230 M-NNY1/NNY2 

EDDY 

124 

FKLEb-1./ (l.+5.5*(CKL£B*YD(fl)/YMAX)**6) 

EDDY 

125 

V1ST0(M)*(<K*CCP*R0(M^N) + FWAKE^FKLEB 

EDDY 

126 

IOUT-0 

EDDY 

127 

00 2A0 M«NNY1#NNY2 

EDDY 

128 

DVIST-VISTO(M)-VISTKM) 

EDDY 

129 

IF(DVIST.LE.O. )IQUT»1 

EDDY 

130 

VliT(M^N)-VISTI (M) 

EDDY 

131 

IF( rOUT.EQ.l)VIST(M>N)-VISTG(M) 

EDDY 

132 

CONTINUE 

EDDY 

133 

CONTINUE 

EDDY 

134 

CONTINUE 

EDDY 

135 

RETURN 

EDDY 

136 

END 

EDDY 

137 

SUbROUTINE IMPY(IAOD) 

IMPY 

1 

C0MM0N/F1/NXM>NXMI>NXM2,NXM3jNXM 4>N1> M1>N11> Ml NS 1> NS 2 /NS > MS > MS 1 

IrtPY 

2 

l^NCWL/NCWLl# NCWL2 ^ nCWLM>NCWL P>NSTRUT 

IMPY 

3 

C0MM0N/FA1/M12#M13^ Nl2iN13>NXM5#NXM6iMS2#MS3>MS4>MSPi>MSP2#MSP3 

IMPY 

4 

l>NSMl^NSlPlf NS2Ml>NSlrtl>NS2Pl^NNS^NSlM2>NXM7 

IMPY 

5 

CDMM0N/F2/RD(61^55)#U(61>55) >V(61>55),P(61,55)^T(61>55)^ 

IMPY 

6 

i SH(6i#55)#H(61^55)/RUS(55)#US(55)/VS(55)/ 

IMPY 

7 
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2 P$(55)#TSC55),SHS(55),H$(55) IMPY 

CLif1iiON/F3/VISL (61, 55), VIST(61,55) ,VIS(61, 55) ,VISS( 55) IMPY 

C0MM0N/F4/SIGX(61,55), SIGY(6i,55)>TAUXY(61/55)#QX(61#55)/ IMPY 

1 QY(61,55),$IGXS (55),SIGYS(55) , TA t XY S ( 55 ) , QX 5 ( 5 5 ) , QY3 ( 5 5 ) IMPY 

COMiON/F5/AU(61, 5 5> 4 ) , AM ( 61, 55 > ^ ) / AM ( 61> 5 5# 4 ), AUl ( 61, 5 5, 4 ) IMPY 

C0MMQN/F6/AUS(55,4),AMS(55f ^), AMS (55,4), AUbl (55,4) IMPY 

C0MM0N/F7/X( 61,55), Y(61,55),XX1(61,55),YXI(61,55), IMPY 

1 XETA(61,55),YETA(61,55),AJ(61,55) IMPY 

CDMM0M/F8/XS (55),YS(55 ),XSXI (55),YSXI (55) ,XSETA(55), IMPY 

i YSETA(55),AJS(55) IMPY 

CUrtM0N/F9/IMT(55),TEl(55),TE2(55),TE3(55),TE4(55) IMPY 

C0MM0N/F10/DT(61,55),t)Trt IMPY 

C0Mf1QN/F12/LSYM,NFL0W, LFI , LF 0, 1 LT, LM AX, L W, L£ XI T, IT, FOT, FOIL, IMPY 
1 TIM£,CCP,CCT, L, LK, TIMEL, CRITAYG IMPY 

CDMMQN/F13/GAMA,RGAS,CP,PR,PRT,CV1S,FM,PF,TF,UF,R0F,$HF,HF IMPY 

C0Mrt0M/F14/DAU(61,55,4),B0AJ(55,4 ),SDAU3(55, 4),TI(55,5), IMPY 

lCCI.iPY,GMl,GM2,GM3,GM4,GM5 IMPY 

gX(l,l5MXM)»GM5*T( 1,UMXM) IMPY 

gX(i,liNXM)»VSaRT(UX(i, ljNXM)|QX(l,l INXM) ) IMPY 

SIGX(1,UNXM )-.5^(U(l, ljNXh)+U(l, ljMXM) + V(l, IjNXM) (1, 1;NXM ) ) IMPY 
QY(l,l|NXM)«OTM/YETA(i,i;MXM) IMPY 

TAJXY(l,l;NXM)«GM4tVlS (1, 1; NXM ) 'I'OY ( 1 , 1 ;NXM ) / (RO ( 1, 1 )MXM ) ♦ IMPY 

lYETAU,liNXM>) IMPY 

SIGY(1,1)NXM)*VABS (V(l,l;MXf1 )IS16Y(1,1INXM) ) IMPY 

Ah(i,l,l)MXM)-SIGY(l,l; NXM)«QYa,l;NXM)-.5 IMPY 

6IGY(i,i;NXM)-(5IGY( 1, l;NXM)+QX(l,li MXM) )tCiY (l,ljNXM)-. 5 IMPY 

AM( i,i,2|NXM )-V( 1, l|NXh)+QX( 1, lINXM) IMPY 

Ah(i,l,2)NXM )-(\/A3S(AM(l,l,2>NXM) jAM(l,l,2;NXM) ) )^QY(1,1 jNXM)-.5 IMPY 
AM(1,1,3JNXM)»V( 1, 1;MXM)-*QX( l,liMXM) IMPY 

Ah(i,i,3lMXM )■ YABS(AM(1,1,3;NXM) jAM(l,l,3lNXM))^QY(i,l|NXM )-.5 IMPY 
MIM-0 IMPY 

NM-2 IMPY 

NN2-NS1M1 IMPY 

IF(r(STRUT. EQ •0)SM2-N11 IMPY 

15 CONTINUE IMPY 

NlM-NIM+1 IMPY 

DO 10 N-NN1,NN2 IMPY 

DO 10 M»2,M11 IMPY 

rtM»(l-IADD)*M+IADD*(Mil+2-M) IMPY 

IF(S1GYCMM,N).LT.0, )G0 TO 10 IMPY 

VV»7(MM,N) IMPY 

CC-gX(MM,N) IMPY 

ALP-SIGX(MM,N) IMPY 

RCCSQ«1./(CC^CC) IMPY 

UU«0(MM,N) IMPY 

RR-RO(MM,N) IMPY 

IF(IADD.EQ.0)G0 TO 30 IMPY 

IF(Mrt.NE,Mll)G3 TO 20 IMPY 

DO 110 I«l,4 IMPY 

110 DAJ(hll,N, I )«DAU(M11,N, I)+BDAU(N, I)tUY(Ml,N) IMPY 

GO TO 20 IMPY 

30 CONTINUE IMPY 

IF(MM.NE,2)G0 TO 20 IMPY 

DO 230 1-1,4 IMPY 

230 DAJ(2,N,I )«DAU(2,N, I)+BDAU(N,I )*0Y(1,N) IMPY 

20 CONTINUE IMPY 

0P«(AlP + DAU(MM,N,1)-UU*DAU(MM,N,2 )-VV*DAJ(MM,N,3) + DAU(MM,N, 4) )^GM1IMPY 

xxi-dauh«,n,i )-op*r::sq impy 

XX2- (DAU( MM, N, 2 )-UU^DAU(MM,N,l ) )/ RR IMPY 

XX3--VV^CC + DAU(MM,N,1)+CC*0AU(MM,N,3 )+0P IMPY 

XX4--XX3+2,*DP IMPY 

VIS1-TAUXY(MM,N)+ABS(XX1) /(GM1*RR/GAMA)*CCIMPY IMPY 

ODl-AMAXK (VlSlf AM(MM, N, 1) ), 0. ) IMPY 

0D3-AMAXK (VIS1^AM(MM,N,2)),0.) IMPY 

DD4-AMAX1 ( ( VIS1+AM(MM, N,3 ) ),0. ) IMPY 

YYi«XXl/(l.+DDl) IMPY 

YY2-XX2/(1,+D01) IMPY 

YY3-XX3/ (1,+D03) IMPY 
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YYA-XX4/ (1.+D04) IMPY 

OK-YYH*. 5^( YY3 + YY^)*RCCS0 IMPY 

OAJ(Mrt>NU )-DR IrtPY 

DAU(MM,N>2)-UU*DR+PR*YY2 IMPY 

DAJ(M«,N>3)* VV^DR^-.S^C YY3-YY4)/CC IMPY 

0AJ(MM#N,4)-ALP*0R+RR^UU + YY2 + .i?*( VV* ( YY3-YY4 ) / C C + ( YY 3 + Y Y ^ ) / GMl) IMPY 
IF(M.Ea.Mll)GO 10 IMPY 

Mj*MM-2tIAC)0 + l IMPY 

DD-YETA(MM,M )/YETA(MJ#N) IMPY 

ZZ1*UD^DD1+YY1 IMPY 

ZZ2-DD+DD1+YY2 IMPY 

ZZ3-DD*DD3^YY3 IMPY 

ZZ^-D0*00^*YY4 IMPY 

Ok-ZZl+.S* (ZZ3+ZZA)*RCCSa IMPY 

OAU(MJf N, l)«DAU(MJ^Nil )+DR IMPY 

DAJ( 2 )-DAU(MJ/N^2 ) + JU + 0R + ZZ2*RR IMPY 

0Au(MJ/N#3)-DA’J(MJ#N>3)>VV + JR+.5* (ZZ3-ZZ4)/CC IMPY 

OAJ(MJ#N^ A) -OAU(MJ^N>A )+ALF*DK+RR»UU^ZZ2+.5* ( VV*(ZZ3-ZZ4)/CC+ IMPY 

1(ZZ3‘»*ZZ4)/GM1) IMPY 

10 CGMTiNUE IMPY 

IF(SSTRUT,EQ,0)R6TURN IMPY 

NM«NS2P1 IMPY 

NK2-N11 IMPY 

IP(fUM.EQ.l )G0 TO 15 IMPY 

Ob 5 N«NS1,NS2 IMPY 

NN«N-NS1+1 IMPY 

5 QYS(N)-DTM/YSETA(MN) IMPY 

Ob AO N«NS1>NS2 IMPY 

Db A5 M«2>MS1 IMPY 

MM*( 1-IADD)*M+IADD+ (Mil+2-M) IMPY 

IF(SlGY(MM,m.LT.O. )G3 TU A5 IMPY 

(MM» M ) IMPY 

CC*4X(MM>N) IMPY 

ALP-SIGX(MM,N) IMPY 

RCCSQ-1./(CC*CC) IMPY 

UU*J(MM#N) IMPY 

Rk-RO(MM,N) IMPY 

IF(IADD,EO,0)GO IG 65 IMPY 

IF(MM.NE.MS1 )G3 TG 60 IMPY 

00 7C I«1^A IMPY 

70 OAj(MSl#S>I) -DAO( MSl^M/ I ) +&0 AOS ( I )^OY:>(M ) IMPY 

Gb rO 60 IMPY 

65 COMTINUE IMPY 

IF(MM,KE.2)G0 TO 60 IMPY 

00 75 I»lfA IMPY 

75 0AJ(2^N# I)»DAU(2#N/I)+B0AU(ron*QY(l/N) IMPY 

60 CONTINUE IMPY 

0P«{ ALP*OAU( MM#M / 1 )-UU^DAU( MMi N/ 2 )“ V V *UAU ( MM # N> 3 ) ♦OAU ( MM# N# A ) )>*'GMiIMPY 
XX1«DAU(MM#N#1 )-D?+RCCSQ IMPY 

XX2« COAU(MM#N# 2)-UU*DAU(MM#N#l ) )/Rk IMPY 

XX3--YV^CC^DAU(MM#N#i) +C:*uAU(MM#K#3 )-i'DP IMPY 

XXA»-XX3+2.^DP IMPY 

VISl-TAUXY (MM#N)^A6S(XX1) /(GM1*RR/GAMA)*CCIMPY IMPY 

ODl-AMAXK (VISl>AM(MM#N/i) ), 0. ) IMPY 

DD3-AMAX1 ( (VIS1+AM(MM#N#2) )#0. ) IMPY 

DLA-AMAXK (^ISi + AM(MM#N#3))#0. ) IMPY 

YY1-XX1/(1,+0D1) IMPY 

YY2-XX2/ (1,+D01> IMPY 

YY3-XX3/ (1.+DD3) IMPY 

YYA-XXA/(1.+DDA) IMPY 

DR-YY1+,5*(YY3+YYA)+RCCSU IMPY 

DAJ(MM#N#1)«DR IMPY 

OAU(MM#N# 2)-UU*DR + RR*YY2 IMPY 

0AJ(Mn#N#3)-VV*0R+.5*( YY3-YYA) /CC IMPY 

0AJ(Mf1# N# A )-ALP^DR + RR*UU*YY2 + .5*( VVt (YY3-YYA ) /CC ♦ ( YY3 + YY A ) / GMl ) IMPY 
IF(M.E0*MS1)GQ TJ A5 IMPY 

MJ-MM-2+IADD+1 IMPY 

OD»YETA(MM#N)/YErA(MJ>N) IMPY 


76 

77 

78 

79 

80 
81 
62 
83 
8A 
85 
66 

87 

88 

89 

90 

91 

92 

93 
9A 

95 

96 

97 

98 

99 
100 
101 
102 
103 
lOA 

105 

106 

107 

108 

109 

110 
111 
112 
113 
ilA 

115 

116 

117 

118 

119 

120 
121 
122 
123 
12A 

125 

126 

127 

128 

129 

130 

131 

132 

133 
13A 

135 

136 

137 

138 

139 
lAO 
lAl 
1A2 
1A3 


40 



APPENDIX A 


45 


90 

85 


95 

80 


50 

40 


ZZ1«DO*D01*YY1 

irtpy 

144 

ZZ2«D0*DD1*YY2 

irtpy 

145 

ZZ3-0U*0D3*YY3 

ifipy 

146 

ZZ^«DD*D0^+YY4 

iMpy 

147 

DR-ZZl+.5*(ZZ3 + ZZ«t)*RCCS3 

Mpy 

148 

DAU(MJ»N,l)«DAU(MJ»N,i)tDR 

iMpy 

149 

0AU(MJ»N»2)»DAU(MJ#N»2) +UU*0R+ZZ2+kR 

iMpy 

150 

DAU(Mj,N,3)«0AU(»1J>N»3)+y/V*DR+.5+tZZ3-ZZA)/CC 

irtPY 

151 

0AU(MJ,N»4)«0AU(MJ,N> A )+ALP*DR+RR*UU*ZZ2+.5*(VV*{ZZ3-ZZA)/CC+ 

IMpy 

152 

1(ZZ3+ZZA)/GM1) 

IMPY 

153 

CaNTINUE 

IMPY 

154 

OCi 50 M-MSP1»M11 

IMPY 

155 

MM-{ l-IA0D)+M>IA0D*(fUl + f1SPl-M) 

IMPY 

156 

IF(SIGY(MM»N).LT.O. )Ga TO 50 

IMPY 

157 

VV«V(rlM»M) 

IMPY 

158 

CC-JX(MH»N) 

IMPY 

159 

ALP-S1GX(MM»N) 

IMPY 

160 

RCCSQ-1./(CC*CC) 

IMPY 

161 

U0-0(MM»NJ 

IMPY 

162 

Rk-RO(MM»N) 

IMPY 

163 

IF(iAOD.EQ.O)GO TO 85 

IMPY 

164 

lF(rtf1.NF.Mll)G0 T3 80 

IMPY 

165 

00 90 I-1»A 

IMPY 

166 

OAJ(hll,N»I )-DAU(Mll,N, I)+-BOAO(N» I)*CY(M1»N) 

IMPY 

167 

GO TO 80 

IMPY 

168 

CO.'ITINUE 

IMPY 

169 

IFCflM.NE.MSPDGD TO 80 

IMPY 

17C 

00 95 1-1»A 

IMPY 

171 

DAU(MSP1»N»I )-0A0{h5Pl»N» I)+B0A0S (N, I)*QY{hS/N) 

IMPY 

172 

CONTINUE 

IMPY 

173 

DP-(ALP*DAU(MM>N>1 )-Uo*DAU(MM/N#2 )-VV*DAO{NM/N^ 3) + 0AU(MM,N» An*GmirtPY 

174 

XXl«DAU(‘1M»N/U-0P’<‘RCCSa 

IMPY 

175 

XX2-(0AU(HM»N, 2)-UU*0AU(Mrt>N»l))/RR 

IMPY 

176 

XX3--VV*CC*DAU(MM»N»l)+C:*0Au{hN»N»3 )+0P 

IMPY 

177 

XXA— XX3 + 2.*0P 

IMPY 

178 

VISl-TAUXY(MM>N)+ABS(XXl)/(GMl*RR/GAhA)*CCIMPY 

IMPY 

179 

0Dl-AMAXl((VISl + AM(MM»N»in»O. ) 

IMPY 

180 

0D3-AMAXK (VISl + AM(Mh»N/2) )>0. ) 

IMPY 

181 

DDA-AMAXK ( VIS1 + AM < MH> N> 3 ) ) > 0. ) 

IMPY 

182 

YYl-XXl/ (l.+DOl) 

IMPY 

183 

YY2-XX2/(1.+DD1) 

IMPY 

184 

YY3-XX3/(1.+DD3) 

IMPY 

185 

YYA-XXA/ (1,+DOA) 

IMPY 

166 

0k-YYl+.5*(YY3+^YYA)*RCCSQ 

IMPY 

187 

0AU(«M»N»1)-DR 

IMPY 

188 

0AU(MM»N»2)«UU*0R+RR*YY2 

IMPY 

189 

DAO(MM»N, 3)«VV*0R+.5*( YY3-YY A) /CC 

IMPY 

190 

DAJ{hH»N, A)«ALP*0R + RR*UU<‘YY2 + . 5*( VV* ( YY3-YYA )/ CC + ( YY 3 + Y Y A ) / GMl ) 

IMPY 

191 

IFtM.EQ.MlDGO TO 50 

IMPY 

192 

MJ-.1M-2*IAD0+1 

IMPY 

193 

DD-YETA(MM,N )/YETA(MJ»N) 

IMPY 

194 

ZZ1-0D*DD1*YY1 

IMPY 

195 

ZZ2-DD*0D1*YY2 

IMPY 

196 

ZZ3-0D*D03*YY3 

IMPY 

197 

ZZA-OD*OOA*YYA 

IMPY 

198 

Ok-ZZl+.5*(ZZ3+ZZ4»*RCC5Q 

IMPY 

199 

0AU(MJ»N,1)-0AJ(MJ>N,1)+0R 

IMPY 

200 

OAU(MJ»N,2)“DAU(MJ»N»2)+UO*OR+ZZ2*RR 

IMPY 

201 

OAU(MJ»N»3)-OAJ(MJ»N»3)>VV'*‘OR+.5'f(ZZ3-ZZA)/CC 

IMPY 

202 

OAJ(MJ»N» A)«DAU(MJ»N»A )+ALP + 0R+RR*UU’»ZZ2+.5* (V V* ( ZZ3-ZZA ) /CC + 

IMPY 

203 

1(ZZ3+ZZA)/GM1) 

IMPY 

204 

continue 

IMPY 

205 

CONTINUE 

IMPY 

206 

RETURN 

IMPY 

207 

END 

IMPY 

208 

SUBROUTINE VECIM2(IA0D) 

VECIM2 

1 

C0M,1UN/Fl/NXH,NXM1>NXM2»NXM3»NXMA»M/Mlf Nli> Ml NS NS2# NS» MS# MSI 

VECIM2 

2 

l>NCrfL»NCWLl»NCi<L2#NCWLM»NCWLP#NSTRUT 

V6CIM2 

3 
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C 


10 

20 


16 

15 


17 

30 


C 

C 


COflMClN/F Al/M12>m3>N12/N13» NXM5#.SXh6/fli)2» MS3>hSA>rtSPi»MSl'2»MSP3 
l»NSMl»NSlPl>NS2>11»NSlf11»Ni2Pl/NNS>NSlrt2>NXh7 
C0MM0N/F2/R0(61,55)>U(61,55)>V(61»55)>P(61>55),T(bl»55)» 

1 SH(61»55)/H(61»55)/ft£JS(55)/Ui(55)/VS(55)# 

2 Pi(55),TS(55),SHS(55)#HS(55) 

C0N.1QN/F3/VISL(61,55), V1ST(61,55),V1S(61,55) /VISS( 55) 
CUMrtaN/F4/SIGX(61,55)#SI6Y(61/55)>TAUXY(61»55)>QX(61#55)> 

1 QY(6i>55)>SIGXS(55)>SlGrS(55) » TAOXYS (55 ) / UX S ( 5 5 ) < QYS ( 55 ) 
CUf1M0N/F5/AU(61»55#4)»Ah<61/55/A)>AN(61»55/A )/AUl( 61/55/4) 
C(JMiiQN/F6/AUS(55>4)/Arli(5 5»4)/AN5(55/4J/AUSl (55/4) 
CilMrt0N/F7/X(61/55)/Y(61/55)>XXl(61/55)/YXI(61/55)/ 

1 XETA(61/55)/YETA(61/55)/AJ(61/55) 

CUMrtQN/F8/XS (5 5 ) , YS (55 ) » XSX 1 ( 55 ) > YSX I ( 55 ) / XS ET A ( 55 ) / 

1 YiETA(55)/AJS(55) 

C0MMQN/F')XINT(55)/TE1(55)/Tt2(55)/TE3(55)/T£4(55) 

C0Mrl0N/F10/DT(61/55)/DTM 

CUMiiOM/FlE/LSYM/NFLOW/LFl/LFO/lLT/LMAX/Lw/LEXIT/CRir/FOT/FOU/ 

1 TlflE/CCP/CCT/L/LK/TlMEL/CRlTAVG 
C0f1MUN/F13/GAHA/RGA5/CP/PR/?RT/CVlS/F«/PF/TF/UF/ ROF/SHF/HF 
C(jflM0N/F14/DAU(61» 55/ 4 ) / BDAU ( 55/4 )/BDAOS ( 55/ 4)/TI(55/5)/ 
1CCI.1PY/GM1/GM2/GM3/GM4/6M5 

IFdADD.EQ.l )G0 TO 10 
K-Hl 

GO TD 20 
K-1 

CONTINUE 

TI ( l/i;Nl ) - Q8 VGATHP(RlI(K/ 1/ NXM) /Ml/ N1 ;TI ( 1/ 1 /N1 ) ) 

T1 (1/2JN1)»Q8VGATHP(U(K/ U NXM ) / Ml/ Nl JTI ( 1/ 2 J N1 1 ) ) 

Tl(l/3;Nl)»Q8VGATHP(V(K/liNXM)/Hl/Nliri(l/3/M)) 

T1(1/4;M)»Q8VGATHP(T(K/1;NXM)/M1/)<1;TI(1/4;N1)) 

Tl ( 1/5)»U ) -Q8VGATHP( VIS (K/liNXM)/ Ml/ Nl) Tl ( 1/ 5 j Nl ) ) 
1 F(.nSTRUT.EQ.0)GJ TQ 30 

THE SUBROUTINE ASSUMEi U-V«0 CONDITIONS ON THE 
STRUT SURFACE. 

IFdADD.EQ.l )GJ TO 15 
DO 16 N-NS1/NS2 
NN-N-NSl+1 
ANS(N/1)«RDS (NN) 

ANS(N/2)-TS(NN) 

ANS(N/3)«VISS(NN) 

GO TO 30 

continue 

DO 17 N-NS1/NS2 
ANS(N/1)«R0(MS/N) 

ANS(N/2)»T(MS/N) 

ANS(N/3)«V1S(MS/N) 

CONTINUE 

RL(2/2)NXM4)»AUl{2/2/l;NXM4)/AJ(2/2/NXM4) 

U( 2/2)NXM4) ■AJl(2/2/2/NXM4)/AUl( 2/2/1/NXM4) 
V(2/2)NXH4)-AUK2/2/3)NXM4 )/a 01(2/2/1)NXM4) 
QX(2/2;NXM4)«(U(2/2)NXM4)*U(2/2)NXM4)+V(2/2)NXM4)*V(2/2)NXM4) )/2. 
SH(2/2;NXM4)-(AU1(2/2/4)NXM4)/AU1(2/2/1;NXM4 )-QX ( 2/ 2)NXM4 ) ) *GAMA 
H(2/2)NXM4 )»SH(2/2)NXil4)+QX(2/ 2)NXM4 ) 

P(2/2)NXH4 )-6M2FR0 (2/2)NXM4)*5H(2/2)NXM4) 
r (2/2)NXM4) ■SH(2/ 2)NXM4)/CF 

CALL bOUND 

00 40 N«2^N11 

ALP».5*(TI(N/2)*TI(N/2)+ri(N/3)*TI(N/3)) 

CC"SUkT(GM5*TI (N/4) ) 

AMS(N/1)-R0(K/N)-TI(N/1) 

AM3(N/2)-RO(K/ .'n*U(K/N)-TI(N/i)*TI(N/ 2 ) 
AMi(N/3)»R0(K/N)*V(K/N)-TI(i'l/l)*TI (N/3) 


\/EClM2 

4 

VECIM2 

5 

VEC1M2 

6 

VECIM2 

7 

VECIM2 

8 

VECIM2 

9 

VECIH2 

10 

VECIM2 

11 

VEC1M2 

12 

VECIM2 

13 

VECIM2 

14 

VECIM2 

15 

\/ECIh2 

16 

\/ECIM2 

17 

VECIM2 

18 

VECIM2 

19 

VECIM2 

20 

VECIM2 

21 

VECIM2 

22 

VEClrt2 

23 

VECIH2 

24 

VECIM2 

25 

VEC1M2 

26 

VECIM2 

27 

VECIfia 

28 

VECIM2 

29 

VECIM2 

30 

\^ECIM2 

31 

\/ECiH2 

32 

VECIM2 

33 

VECIM2 

34 

VECIH2 

35 

VECIM2 

36 

\/ECIM2 

37 

VECIrt2 

38 

VECIf12 

39 

VECIM2 

40 

VECIM2 

41 

VECIN2 

42 

VECIH2 

43 

VECIM2 

44 

VECIf12 

45 

V£Cih2 

46 

VEC1M2 

47 

VECIM2 

48 

VECIM2 

49 

\/EClM2 

50 

VECIM2 

51 

VECIM2 

52 

VEC1M2 

53 

VECIM2 

54 

VECIK2 

55 

VECIM2 

56 

VECIH2 

57 

VECIM2 

58 

VECIh2 

59 

VEClh2 

60 

\/EClM2 

61 

VECIM2 

62 

VEC1M2 

63 

\/ECIh2 

64 

VECIM2 

65 

VEC1M2 

66 

VECIM2 

67 

VECIM2 

68 

VEC1M2 

69 

VECIM2 

70 

VECIM2 

71 
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AMS (N> A)"RQ(Ki M)*(GM3*T(K>N)+.5*(U(K,»M)*U(K»N)+V(K/N)*V(K>N)))- 

in(rta)*(GM3*TI(N>4)+ALP) 

0P«GM1*(ALP*AMS(N»1)-TI (N>2»»AMS(K>2 )-Tl (N/3 )*AMS(N»3)+AMS(N>4) ) 
RCCS^J»1./<CC*CC) 

XXi«AMS(M»l)-OP»RCCSQ 

XX2-(AMS(N,2)-TI (N,2)*AKS(N^ 1) )/TI(N»l) 

XX3«CC*(AMS(N» 3)-TI (N> 3)*AMS (N»1) )+DP 
XX4»-XX3+2.*DP 

V1S1»GM4*T1 (N> 5) / (TI ( N» 1 ) *Y ETA { K» N ) ) -.5*YETA (K> N ) / DIM 

VAiJ;>l«AeS(TI (N,3) )+Vlsl 

VPCABS-ABS(TI(N, 3)+CC) ♦VlSl 

VMCaBS-ABS(TI (N>3)-CC)+VIS1 

0L/1*AMAX1(VABS1>0. ) 

D03-AMAX1(VPCA8S>0. ) 

D04»AMAX1(VMCA8S»0. ) 

YY1-D01*XX1 

YY2-DD1*XX2 

YY3«D03*XX3 

YY4«0D4*XX4 

AYl«(YY3*YY4)/2. 

AY2»(YY3-YY4)/(2.*CC) 

BDaJ{N»1)-YY1+AY1*RCCSQ 

BUAU(N»2»-TI (N>2)*B0AJ(N, 1)+TI(N,1)*YY2 
BDAU{N>3)»TI (N>3 )*BDAU (N> D + AY2 

BDAd(N»4)»ALP*aDAU(N>l)+TI(N»l)*Tl (N»2)*YY2 + TI (N»3 )*AY2 + AY1 /GMl 
CONTINUE 

IF(NSTRUT.EQ.O)RETURN 
IFdADO.EQ.DGQ TO 55 
DU 50 N«NS1>NS2 
NN-N-NSl+1 
TEKN)-ROS(NN) 

TE2(N)-TS(NN) 

TE3(N)«YSETA (NN) 

60 TO 60 

CONTINUE 

00 56 N-NS1>NS2 

TE1(N)»R0(MS>N) 

TE2(N)-T(HS»N) 

TE3(N)-YETA(MS>N) 

continue 

DO 65 N-MS1/NS2 
BDAJStN,l)«TEl(S )-ANS(N»l) 

B0AJS(N,4)«GM3*(TE1(N)*TE2(N)-ANS (N> 1 )*ANS (N>2 ) ) 

YlSl“GM4*ANS{N>3)/(ANS(N»l)*TE3(N))-.5*rc3(N)/DTM 

TE4(n)-S3RT(GM5*AS5(N/2)) 

TI(N>1)-AMAX1(\/1S1>0. ) 

TI(N>2)«AMAX1(V!S1+TE4(N)/0.) 

BDAUS(N» l)»TI(N»l)»aDAJS(N, 1 ) -GMl / ( T E4 (N )*TE 4( N ) )*(TI(N»1 )-TX(Nj 
1)*BDAUS(N,4) 

BDAJS(N»4)-TI(N» 2)*BDAUblN»4) 

> CONTINUE 


RETURN 

END 

SUBROUTINE 


SPILKSwEEP ) 


1*NCWL^NCWL1»NCWL2,NCWLM/NCWlP,NSTRUT 
C0MM0N/FAlYM12»M13>N12»N13>NXM5»NXM6>rtS2#f1S3/N54/rtSPl/MSP2/ MSP3 
1# NSMl^NSlPlj NS2M1> NGlfU»ri5 2Pl»NN5^NS IM2,NXM7 
CON(10N/F2/R0(61» 55)»U(61»55)>V(61»55)»P(61#55)^T(61/55)/ 

1 SH(61»55)»H(61»55J,kOS(55)»Ui(55)/VS(55), 

2 PS(55),TS(55),SHS(55)/HS(55) 
C0NN0N/F4/SIGX(61»55)>SIGY(61»55)/TAUXY161/55)/QX(61»55)^ 

1 0YI61/55) /S IGXS ( 55) >S IGYS ( 55) # TAUXYS {55)>QXS( 55)» QYi ( 55 ) 
OjMMDN/F7/X( 61» 55 )>Y (61,55 )» XXI (61»55)»YXI (61,55 )» 

1 XETA(61,55),YETA(61,55),AJ(61,55) 

C0i1MDN/F8/XS (55),YS(55 ),XSXI (55),Y5XI(55),XS£TA(55), 

1 YSETA(55),AJS(55) 


VECIH2 

72 

VEC1M2 

73 

VECIM2 

74 

VECIM2 

75 

VECIH2 

76 

VECIM2 

77 

VECIH2 

76 

VECIM2 

79 

VECIM2 

80 

V£CIh2 

01 

VECIM2 

82 

VECIM2 

83 

VEClrt2 

84 

VECIM2 

05 

VECIH2 

06 

VECIM2 

87 

VECIM2 

80 

VECIrt2 

89 

VECIM2 

90. 

VECIM2 

91 

VECIM2 

92 

VECIM2 

93 

VECIM2 

94 

VECIM2 

95 

VECIM2 

96 

VECIM2 

97 

VEC1M2 

98 

VECIM2 

99 

VECIM2 

100 

VECIM2 

101 

VECIH2 

102 

VECIH2 

103 

VECIM2 

104 

VECIrt2 

105 

VECIM2 

106 

VECIH2 

107 

VECI«2 

108 

\/ECIh2 

109 

VECIh2 

110 

VECIM2 

111 

VECIf12 

112 

VECIf12 

113 

VEC1W2 

114 

VEC1M2 

115 

VECIM2 

116 

VECIH2 

117 

VECIH2 

118 

2)VECIM2 

119 

VECIM2 

120 

VECIM2 

121 

VECIM2 

122 

\/ECI«2 

123 

VECIM2 

124 

VECIM2 

125 

SPILL 

1 

,1 SPILL 

2 

SPILL 

3 

SPILL 

4 

SPILL 

5 

SPILL 

6 

SPILL 

7 

SPILL 

8 

SPILL 

9 

SPILL 

10 

SPILL 

11 

SPILL 

12 

SPILL 

13 

SPILL 

14 
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30 

40 

to 

45 

55 

35 

275 

280 

255 


COMf10N/F12/LSYM»NFLOW»LFl/LFO/ILT/LnAX>LW>LEXIT>CRIT>FOT>FDTL# 

1 T I«E>CCP>CCT>L»LK>TlM£L,CftlTAVG 
C0i1.10N/F13/GAMA»ftGAS/CP»Pft/PkT#CVIS> FM/PF^ TF>UF>RaFtSHF>HF 
DIMENSION XP(61»55)/LiiP<61/55)>USTl<65)»USPl(55) 

DIMENSION FMASS(55)>UST(61,55)>FMAST (55) 

P1»4.*ATAN(1.) 

THESh«SWEEP*PI/180. 

CTHESW-CDS(THFSrf) 

STHcSi«-SIN(THES»0 
FMS» FM/CTHESW 

XP(i;ljNXM)«X(l,UNXM)/CrHESW 

VT"FNS*SQRT(GAMA*RGAS*TF)*STH£SW*CTHESW 

osn lINXM ) “VT 

USn(lJNS)«VT 

IF(NFLOW.EQ,0)GO to 35 

0X(1>1;NXM)»VABS(U(1/1}NXM)j0X(1»1|NXM)) 

IN Mli2»Ml-12/ 12 IS AN ARBITRARY NUMBER. WHAT IT MEANS IS THAT 
rIRST 12 POINTS LIE IN THE B.L. PROFILE. THIS IS REOUIREO TO PUT 
VISCuUS EFFECTS ON THE INYISCID COMPONENT VT . NUMBER 12 MAY 
NEED TO BE CHANGED FOR CERTAIN CALCUL ATIONb . 

M112-M1-12 
00 30 N«1/N1 
DU 30 M-M112iMl 

0ST(M>N)«VT+QX(M»N)/QX(M12»N) 

IF(NSTPUT.E0.0)G0 to 45 
0XS(1JNS)"VABS(US(1;NS )} OXS ( 1;N5 ) ) 

MSll-MS+11 
DO 40 N-NS1^NS2 
00 40 M«MS^MS11 

UST(M,N) -VT^QXCMjNl/QXtMSll^ N) 

MSlO-MS-10 

DO 50 N«NS1»NS2 

NN-N-NSl+1 

DO 50 M»MS10>MS1 

QYS(NN)»0X(MS10^N) 

UST(M^N)»VT*QX(M,N)/QX(MS10jN) 

USTl(lJNS)»YT+QXS<lINS)/OYS(liNS) 

USP1(1;NS)«UST1(1;nS)-US(1;nS)*STHESW 

CONTINUE 

IF(NC«L.EQ.0)G3 Tfl 35 
DO 55 N«NCWL1# NCWL2 
Ou 55 M«l,12 

U5T(M,N)«VT+QX(M,N)/0X(12^N) 

continue 

USP( 1>1 JNXM) «UST( 1» 15NXM)-U( 1» 1JNXM)*STHESW 
NSM-NSIMI 

IF(NSTRUT.EQ.0)NSM-N1 

fmass(i;ni)»o. 

DC 275 N«1>NSM 
00 275 M-2>M1 

FMASS(N)»FMASS{N)+(R0(M-1>N)*USP(M-1#N)+R0(M>N)*UsP(M>N) )*(Y(M»N) 

1Y(M-1»N) )/2. 

CONTINUE 

IF(NSTRUT.EO.O)GO TO 290 
Ou 280 N-NS2P1»N1 
DO 280 M«2>M1 

FMASS(N)"FMASS(N)+(R0(M-1»N)*USP(M-1»N)+R0(M/N)*USP(M»N) )*( Y(M/N)' 
1Y(M-I,N) )/2. 

CONTINUE 

DO 250 N-NS1#NS2 
00 255 M»2>MS1 

FMASS (N) -FMASS (N )+ (R0( M-1>N) *USP ( M-1 >N) + R0 ( M » N ) P ( M»N ) )♦ ( Y ( M» N )• 
1Y(M-1,N) )/2. 

continue 

DU 260 M"MSP1^M1 

FMASS(N)»FMASS(N)+(R0(M-1#N) ♦USP( M-1 >N) + R0{M>N ) ^Uip (M>N ) )*( Y( M» N)- 
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260 


250 

290 


285 


300 


805 


810 


20 

806 


1Y(M-1>N) )/2. 

CCNTlNUE 
MN-.S-NSl + l 

FHASS(N)-FhASS(N)+(ROS(NN)*USHl(NN)+ 

1RC(M31»N)*USP(MS1»N) )♦( YS (NN )-Y(MSl» N ))/2. 

CONTINUE 
CCNTlNUE 

FFAST(l)-FhASS(l) 

DO 285 N»2»N1 

FMA3T(N)«FMAST(S-l) + (FMASS(N)+FflA5S{N-l) ) ♦ ( X F( 1» N ) -XP ( 1# N-1 ) ) /2 . 
CONTINUE 

WRITE(6>300>(N>X?(1/N)»FMAS6 (N ) » F HAS T ( N » / N«1 >N1) 

FCRMATI 10X> ' BODY STATION NO . ■ ' > 1 4 » 3X » ' XP« ' / F 10 . 6/ 3 X» ' LUC AL MASS St 
1ILLE0«'»F10.5»3X, 'TOTAL MASS SPILLED UPTO XP«'>F10.5» 
OX(l»i}NXM)«P(l» 1;NXM)/K la ) 

WRlTEt6>805) 

RRITE(6/810)<N»XP(l>N)>QXll^N)#QX(Ml/N)>N»l/Nl) 

F0RMAT{///»5X»*'t'»16X» 'XP'»21Xi 'CENTER LINE PR. RATIO'# 

120X> 'SIDEWALL PR. RATIO'#/) 
F0RrtAT(4X»I3»7X»F13.5#20X#F15.A#20X#F15.4) 

IF(NSTRUT.EQ.0)60 TO 20 
0XS(NS1}NS»«PS(UNS)/P(1»1) 

WR1TE{6»806) 

WRITE(6#810) (N#XP<HS#N)#WX(hS#N)» UXS (N)#N«NSl#NS2) 

CONTINUE 

F0RMAT(///»5X#'M'»16X» 'XP'#16X/'STRUT UPPER SURFACE PR. RATIO'# 
116X»'STRUT LOWER SURFACE PR. RATIO'#/) 

RETURN 

End 

SUBROUTINE PRINT 

C0M.10N/Fl/NXM>NXm#NXM2>NXM3#NXf',4#Nl#Ml#Nll# M11#NS1#NS2#NS# rtS#NSl 
1#NC«L#NCWL1»NCWL2 »NCwLM#NCWLP#NSTRUT 
C0rtM0N/FAl/M12#M13#N12#N13#NXM5#NXM6#MS2#f1S3#hS4#MSPl#MSP2#MSP3 
1#NS(11,NS1P1#NS2M1»NS1M1#NS2P1#NNS#NS1M2#NXM7 
C0rtrt0N/F2/R0{61#55)#U(61#55) #V(61#55)#P(61#55)#T(61#55)» 

1 SH(61#55)#H(61# 55)#R0S(55)#US (55)#VS(55)# 

2 P3(55)»TS(55)#SHS155),HS(55) 

C0MM0N/F3/VISL(61»55)#VIST(61#55) #VIS(61#55) #VISS(55) 
C0Nrt0N/F7/X(61#55)#Y(6i# 55)#XXI(61#55)#YXI(61#55)# 

1 XETA(61#55)#YETA(61#55)# AJ(61#55) 

CDrt.10N/F8/XS(55)» YS(55 )#XSXI ( 5 5 ) # YSX I ( 55 ) # XS ET A ( 55 ) » 

1 YSETA(55)#AJS(55) 

C0MMUN/F11/ERR0(61#55)#B1 (6i#55) 

CartM0N/F12/LSYM»NFL0w» LFI#LFU#ILT#LHAX#LW#LEXIT#CRIT#FDT#FDTl# 

1 TIME#CCP#CCT#L#LK#TIf1EL#CRlTAVG 
C0f1NUN/F13/GAMA» RGA3»CP#PR#PRT#CVIS#FM#PF#TF#UF#R0F#SHF#HF 
i«KlTE{6,500)L»riME#TlMEL 
HRITE(6#505> CCP#CCT#F0T#FDTL 
WPITE(6#510) CRIT#CRITAVG#LEXIT 
DO 250 N-1#N1 
WR1TE(6» 545)N 
WRITE(6#550> 

DU 250 MM-1»M1 

rl-Ml-MH + 1 .. 

WRITE(6»560)X(M#N)#Y(M#N)#U(M#N)»V(M#N)#P(h#N)#T(i#N)»R0(rt»N)# 

lSH{f1#N)#VIST(M#N)#ERR0(M»N) 

250 CONTINUE 

IF(NSTRUT.EC.0)G0 TO 255 

WRITE(6»570) ^ . 

WRITE (6, 560 (XS(N)# YS<N)#US(N)# VS{N)#P5(N)#TS(N)#R0S(N)#SHS <N)# 
1HS(N)#VISS(N)#N-1»NS) 

255 CONTINUE 

ERR0(1#1#NxM)»R 0( l#l;NXM) *U( 1#1;NX7 i) 

DO 265 N-1#N1 
FMASS-0. 

IF(NSTRUT.EQ.0)G0 to 261 
IF (N.LT.NS1.0R.N.GT.NS2)G0 TO 261 
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I'iN-N-NSl + l 
DO 263 M«2>MS1 

263 F.1ASS«FMASS + ( ERRO(f1»N )+ERRU(M-l>N) ) / 2.* { Y ( M>N )-Y { h-l» N ) ) 

FMA3i«F7'ASS + (ERR0(MS1> K) + ROS (NN)*USt NN) )/2.* (YS (NN )-Y(HSl>N ) ) 
i)ti 262 K*MSP1>M1 

262 F!1ASS-F‘1ASS + (ERRD(rt, ■1)+ERRJ(M-1>N)) / 2.* (Y(rt>N)-Y( n-l>N) ) 

GO TO 26A 
261 CO.-ITINIJE 

DC 260 M-2»M1 

260 FiiA3i-F MASS + { ERR 0(rt,N ) + ERRO ) ) / 2. + { Y ( f1/N )-Y ( M-1,N ) ) 

26A continue 

wRITE(6^700)N>X(l»N),FMASS 
265 CONTINUE 

500 FURMaT(/> 2X^ 'NO. Ot ITERATIONS ■ ' > 15 » 2 a> ' TI M E •'»E15.7^ 

1 2X, 'TIME LIMIT-'#E15.7) 

505 F0RMAT( /,10X, 'CCP-',F8.A.10X>'CCT-'.Fe.A,10X^ 'FDT» '^FS.A^IOX* 

1 'FDTL-'»F8.AW) 

510 F0Rf1AT(/> 1X> 'LOCAL ERROR CR I TERION- ' #F10, 5/ 5 X# 

1 'AVERAGE ERROR C RITER ION- • / FI 1 . 7 /5X, 

1 'ORDER OF EXTRAPOLATION AT OUTFLOW *•, 12 ) 

5A5 FORMAT!/, 10X#'BODY STATION N0.-',I3,/) 

550 FORMAT ( /,6X, 'X ', 9X, 'Y',11X, 'U',12X/'V',12X,'P',12X,'T',10X,'R0», 
113X, 'SH',11X,'VIST',5X, 'ERRO') 

560 FORMAT (IX, 2F9 . 6, 2F12. 5, F 14, 5, F 13 .5, F 11 . 5, FI A. 5, F 1 5 . 6, F10.6 ) 

570 FORMAT!/, 5X, ' X S ' , 8 X, • Y5 • , lOX , • US • , 11 X, • VS ' , 1 lA, ' PS ' , IIX, *TS ' ,10X , 

1 'ROS ',11X,'SHS',11X,'HS ', 9X, 'VlSS ' ) 

700 FORMAT (/, lOX, • BODY STATION NO . - • , 1 3, lOX, ' X ! M ETER S ) ■ ' , F9 . 5, 1 OX, 
I'MASS FLOW RATE!KG/M/SEC)*',F10.5) 

RETURN 

END 


PRINT 

40 

PRINT 

41 

PRINT 

42 

PRINT 

43 

PRINT 

44 

PRINT 

45 

PRINT 

46 

PRINT 

47 

PRINT 

48 

PRINT 

49 

PRINT 

50 

PRINT 

51 

PRINT 

52 

PRINT 

53 

PRINT 

54 

PRINT 

55 

PRINT 

56 

PRINT 

57 

PRINT 

58 

PRINT 

59 

PRINT 

60 

PRINT 

61 

PRINT 

62 

PRINT 

63 

PRINT 

64 

PRINT 

65 

PRINT 

66 

PRINT 

67 

PRINT 

68 

PRINT 

69 


46 



APPENDIX B 


EXAMPLES OF GEOMETRIES AND ASSOCIATED PARAMETERS 
Nonsi^nmetric Problems 

Exan^les of geometries and associated parameters for nonsymmetric problems are 
as follows: 


Example 1 



A: 

N = 

1, 

M 

= 1 

Nl 

= 

61 

B: 

N = 

1. 

M 

= Ml 

Ml 

= 

51 

B* : 

N = 

2, 

M 

= Ml 

SWEEP 

= 

0 

C: 

N = 

Nl, 

M 

= Ml 

NSTRUT 

= 

0 

D: 

N = 

Nl, 

M 

= 1 

LSYM 

= 

0 






NCWL 

= 

0 

FOR 

NFLOW = 

2: 

CALL EDDY 

(2,N1,1,M1,1,0) 
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Example 2 



A; 

N 

= 

1, 

M 

= 

1 

B: 

N 

= 

1, 

M 

- 

Ml 

B' : 

N 

= 

2, 

M 

= 

Ml 

C: 

N 

= 

Nl, 

M 

= 

Ml 

D: 

N 

= 

Nl, 

M 

= 

1 

E; 

N 

= 

31, 

M 

= 

1 


N1 = 61 
Ml = 51 
SWEEP = 0 
NSTRUT = 0 
LSYM = 0 
NCWL = 1 
NCWLi = 31 
NCWL2 = Ni 


FOR NFLOW = 2 


CALL EDDY ( 2 , 30, 1 , Ml , 1 , 0 ) 
CALL EDDY ( 31 , Nl , 1 , Ml , 2, 0 ) 
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Example 3 

B B' 



A: 

N 

= 

1, 

M 

= 

1 

Nl 

= 

61 

B: 

N 

= 

1, 

M 

= 

Ml 

Ml 

= 

51 

B' ; 

N 

= 

2, 

M 

= 

Mi 

SWEEP 

= 

0 

C: 

N 

= 

Nl, 

M 

= 

Ml 

NSTRUT 

= 

0 

D: 

N 

= 

Nl, 

M 

= 

1 

LSYM 

= 

0 

E: 

N 

= 

41, 

M 

= 

1 

NCWL 

= 

1 

F: 

N 

= 

21, 

M 

= 

1 

NCWLl 

= 

21 


NCWL2 = 41 


FOR NFLOW = 2: CALL EDDY ( 2 , 20, 1 , Ml , 1 , 0 ) 

CALL EDDY ( 21 , 41 , 1 , Ml , 2 , 0 ) 
CALL EDDY (42, N1 , 1 , Ml , 1 , 0 ) 
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Example 4 

B B* 



A: N = 1, M = 1 

B: N = 1, M = Ml 

B ' : N = 2, M = Ml 

C: N = Nl, M = Ml 

D: N = Nl, M = 1 

P: N = 36, M = 1 

I: N = 1, M = 26 

E: N = 16, M = 26 

F: N = 31, M = 26 

G: N = 46, M = 26 

H: N = 46, M = 26 

J: N = Nl, M = 26 


Nl = 61 
Ml = 51 
SWEEP = 0 
NSTRUT = 1 
NSl = 16 
NS2 = 46 
MS = 26 
LSYM = 0 
NCWL = 1 
NCWLl =36 
NCWL2 = Nl 


FOR NFLOW = 2: CALL EDDY ( 2 , 1 5 , 1 , Ml , 1 , 0 ) 

CALL EDDY ( 16 , 46, MS, Ml , 2 , 0) 
CALL EDDY ( 47, Nl , 1 , Ml , 2 , 0 ) 
CALL EDDY ( 16 , 35 , 1 , MSI , 1 , 1 ) 
CALL EDDY ( 36 , 46 , 1 , MSI , 2 , 1 ) 
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Symmetric Problems 

The code solves only the upper half of the flow for symmetric problems. 
Examples of geometries and associated parameters for symmetric problems are as 
follows: 


Exaitple 1 






N1 = 61 
Ml = 51 
SWEEP = 0 
NSTRUT = 0 
LSYM = 1 
NCWL = 0 

FOR NFLOW = 2: CALL EDDY ( 2 , N1 , 1 , Ml , 1 , 0 ) 
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Example 2 




N1 = 61 
Ml = 51 
SWEEP = 0 
NSTRUT = 0 
LSYM = 1 
NCWL = 1 
NCWLl =31 
NCWL2 = 51 

FOR NFLOW = 2: CALL EDDY ( 2 , 30 , 1 , Ml , 1 , 0 ) 

CALL EDDY ( 31 , 51 , 1 , Ml , 2 , 0 ) 
CALL EDDY ( 52 , N1 , 1 , Ml , 1 , 0 ) 
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Example 3 



N1 = 61 
Ml = 51 
SWEEP = 0 
NSTRUT = 1 
NSl = 26 
NS2 = 51 
MS = 27 
LSYM = 1 
NCWL = 0 

FOR NFLOW = 2: CALL EDDY { 2 , 25 , 1 , Ml , 1 , 0 ) 

CALL EDDY ( 26 , 51 , MS , Ml , 2 , 0) 

CALL EDDY ( 26 , 51 , 1 , MSI , 1 , 1 ) 

CALL EDDY ( 52, N1 , 1 , Ml , 1 , 0 ) 

The sample problem discussed in this report provides another example of the type 
of geometries that can be analyzed with the present code. 
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USE OF PARAMETER LKK 

In this appendix, several examples are given for the use of parameter LKK in 
trying more than one value of time-step parameter FDT or daitping coefficients CCP and 
CCT during one computer run. For a given problem, it may be that the calculations 
become unstable (i.e. , a negative temperature occurs in the flow field) for one par- 
ticular value of, say, FDT, and the user may have to try several values before he 
finds the right value of FDT for which the calculations remain stable. The same 
thing may have to be done with CCP and CCT. One way to find the right values of FDT 
and CCP and CCT would be to submit several runs with a particular set of values of 
these parameters, but to avoid this, the code has a parameter LKK at line NASCRIN 224 
that allows the user to try several values of FDT or CCP and CCT in one run of the 
code. To do so, the user has to change lines NASCRIN 219 through NASCRIN 224 as 
illustrated in the following examples: 


Example 1 

FDT = 1. 

CCP = 0. 

1000 CONTINUE 

LKK = LKK + 1 

CCP = CCP + 0.1 

CCT = CCP 

IF(LKK.EQ.6)STOP 

The code begins calculations with FDT = 1 and CCP = CCT =0.1. If the cal- 
culations develop a negative temperature in the field, the code stops, prints out the 
flow-field variables at each grid point, and then returns to statement 1000. It is 
restarted with FDT = 1 and CCP = CCT = 0.2. If the calculations remain stable, 
the code will terminate in a normal way either on a convergence criterion or on LMAX 
without trying any more values of CCP and CCT. But, if a negative temperature devel- 
ops again in the flow field, the calculations will restart with FDT = 1 and 
CCP = CCT = 0.3. For this example, the code will try five values of the damping 
coefficients in increments of 0.1 before stopping. The number of attempts is deter- 
mined from the statement 

IF(LKK.EQ.6)ST0P 

and is always one less than the number in parentheses (in this case, it is 
6-1=5). 


Example 2 

FDT = 1 . 

CCP =0.5 
1000 CONTINUE 

LKK = LKK + 1 
CCT = CCP 
FDT = FDT - 0.1 
IF(LKK.EQ.4)STOP 
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This example allows the user to try three values of FDT, namely 0.9, 0.8, 
and 0.7, with a constant value of CCP = CCT = 0.5. 


Example 3 


FDT = 0.9 
CCP =0.5 
1000 CONTINUE 

LKK = LKK + 1 
CCT = CCP 
IF(LKK.EQ. 2)ST0P 

This example uses only one set of values of FDT, CCP, and CCT, namely FDT = 0.9 
and CCP = CCT = 0.5. 

If the code develops a negative temperature very early in the calculations, the 
chances are that the value of FDT needs to be reduced. But if the negative tempera- 
ture occurs after several hundred time-steps, it may be caused by insufficient damp- 
ing. This may not always be true but provides the user with a guideline in changing 
FDT, CCP, and CCT. Also, if the user wants to try several values of FDT, he should 
start from a higher value and decrease it gradually so that the allowable time-step 
is as large as possible. On the other hand, to keep the damping contribution as 
small as possible, the user should start with smaller values of CCP and CCT and grad- 
ually increase them. 
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ADJUSTMENT OP TIME-STEP AND DAMPING COEFFICIENTS 

In this appendix, examples are presented to show how to gradually increase FDT 
or decrease CCP and CCT over a number of time-steps from their initial values. 


Example 1 

Suppose the calculations are started with FDT = 0.5, and the user wants to 
increase it to 0.9 over 1000 time-steps. Ihis can be done by inserting the following 
lines after line VISCOUS 51 for viscous flow calculations or after line INVICID 60 
for inviscid flow calculations: 

IF(L.GT.1000)GO TO 400 

FDT = FDT + ((0.9 - 0. 5)/1 000) *2. 

400 CONTINUE 

The factor 2 in the esgiression for FDT occurs because the change in FDT is made 
after every two time-steps. By inserting the above three lines, FDT will change from 
0.5 to 0.9 in 1000 time-steps. FDT will then remain constant at 0.9 for subsequent 
time -steps. 


Example 2 

In the second example, suppose the user wants to use a damping coefficient 
of 0.7 for the first 1000 time-steps and then wants to reduce it to 0.5 over the next 
2000 time-steps. Ihis can be done by inserting the following lines in the code at 
the locations discussed in the previous example: 

IF(L.GT.3000)G0 TO 400 

IF(L.LT.1000)G0 TO 400 

CCP = CCP - ((0.7 - 0.5)/2000)*2 

CCT = CCP 

400 CONTINUE 

These statements change the damping coefficients from 0.7 to 0.5 between time- 
steps 1000 to 3000. Ihe damping coefficients remain constant at 0.5 thereafter. 
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SAMPLE OUTPUTS 


Geometry Output 


N 

M 

X 

Y 

XXI 

XETa 

YXI 

YETA 

AJ 

1 

61 

.00000000 

.04000000 

.00500000 

.OOOOOOOO 

.OOOOOOOO 

•00002055 

.00000010 

1 

60 

.00000000 

.03997510 

.00500000 

.oooooooc 

•OOOOOOOO 

.00002983 

•00000015 

1 

59 

•OOOOOOOO 

.03_99-3-8'5’7.. 

. .OO5Q0000 . 

.50005000 

.OOOOOOOO 

.00004324 

•00000022 

1 

58 

.00000000 

.0398B665 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00006259 

•00000031 

1 

57 

.00000000 

.03981099 

.00500000 

.OOOOOOOO 

.OOOOOOOO 

.00009C40 

.00000045 

1 

56 

.00000000 

.0397P187 

,00500000 

.OOOOOOOO 

.OOOOOOOO 

.00013018 

.00000065 

1 

55 

.00000000 

.03954509 

.00500000 

.OOOOCOCO 

•OOOOOOOO 

•00018662 

.00000093 

1 

54 

.00000000 

.03932106 

,00500000 

.OOOOOOOO 

•OOOOOOOO 

.00026582 

.00000133 

1 

53 

,00000000 

♦C 3 90033 9. 

. ...005550110.... 

._i50J35Q505 

,50000000 

•00037522 

•00000186 

1 

52 

.00000000 

.03055786 

•00500000 

.OOOOOOOO 

.OOOOOOOO 

.00052288 

.00000261 

1 

51 

.00000000 

.03794250 

.00500000 

.OOOOOOOO 

.OOOOOOOO 

.00071578 

.00000358 

1 

50 

.00000000 

.03711035 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00095622 

.00000478 

1 

49 

.00000000 

.03601659 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00123655 

•00000616 

1 

48 

.00000000 

.03463154 

.00500000 

.OOOOOOOO 

.OOOOOOOO 

.00153337 

.00000767 

1 

47 

.OOOOOOOO 

.03295600 

.CO5OQ000 

•OOOOOOOO 

.OOOOOOOO 

.00180529 

.00000903 

1 

46 

.00000000 

.03104684 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00199974 

•00001000 

1 

45 

.00000000 

.02900000 

.00500000 

•OOOOOOOO 

•OOOOOOOO 

.00207073 

.00001035 

1 

44 

.00000000 

.02695316 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00199974 

•00001000 

1 

43 

.00000000 

.02504200 

,00500000 

•OOOOOOOO 

•OOOOOOOO 

.00180529 

.00000903 

1 

42 

.00000000 

.02336P46 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00153337 

.00000767 

1 

41 

.00000000 

.02198341 

. 00500000 

.OOOOOOOO 

•OOOOOOOO 

.00123655 

•00000616 

1 

40 

.00000000 

.02088965 

.00500000 

•OOOOOOOO 

•OOOOOOOO 

.00095622 

.00000478 

1 

39 

.0,0000000 

.02005750 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00071578 

•00000358 

1 

30 

.00000000 

.01944214 

•00500000 

.OOOOOOOO 

•OOOOOOOO 

•00032288 

•00000261 

1 

37 

.00000000 

.01699661 

.00500000 

.OOOOOOOO 

.OOOOOOOO 

,00037522 

•00000166 

1 

36 

.00000000 

.C1867894 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

.00026582 

•00000133 

1 

35 

,00000000 

.01845491 

.00500000 

.OOOOOOOO 

. OOOOOOOO 

.00018662 

•00000093 

1 

34 

.OOOOOOOO 

.01829813 

.00500000 

.OOOOOOOO 

•OOOOOOOO 

•00013018 

•00000065 

1 

33 

.00000000 

.01018901 

•00500000 

•OOOOOOOO 

•OOOOOOOO 

.00009040 

.00000045 

1 

32 

.OOOOOOOO 

.01811335 

•00500000 

.OOOOOOOO 

•OOOOOOOO 

•00006259 

.00000031 

1 

31 

.00000000 

.01806103 

.C050006o 

•OOOOOOOO 

•000 000 CO 

•00d04324 

•00000022 

1 

30- 

•OOOOOOOO 

.01002490' 

.00500000 

•OOOOOOOO 

•OOOOOOOO 

•00002963 

•00000015 

u. 

29. 

.. jO’oooooo'o 

.01800000 

.00500000 

.09000000 

•OOOOOOOO 

.00002443 

.00000012 

1 



. nosnnnnn 

. nnnnnnno 


^00OO?94i 

.nnnnnni *i 

1 

27 

.00000000 

.01793941 

.QQiQOOQO..: 

.OQOOOOOU 

•OOOOOOOO 

•00004494 

•00000022 

1 

26 

.00000000 

•01788348 

.Q05000Q0. 

,00000000 

•OOOOOOOO 

•00006653 

.00000034 

1 

25 

,00000000 

.C1779835 

•00500000 

•OOOOOOOO 

•OOOOOOOO 

•00010416 

.00000052 

1 

24 

.00000000 

.01766927 

.00500000 

•OOOOOOOO. 

•OOOOOOOO 

.00015752 

,00000079 

1 

23 

.00000000 

.01747481 

.00500000. . 

.50000000 

•OOOOOOOO 

.00023643 

.00000118 

. 121 . 

_..QOOaQQOQ_ 

.0171846Q 

.QQ5QQQQQ 

.OOOOOOOO .Donooooo 

,005350 81.. 

-.,050110115. 

1 

21 

•OOOOOOOO 

.0167^7^^. 

•00500C00 

,50000000 

-.OOOOOOOO 

.00051199 

•00000256 

1 

20 

.00000000 

.01614187 

_.055QOOOO_. 

.,05050000 

•OOOOOOOO 

.00072909 

•00000365 

1 

19 

.00000000 

.015^8Q27^ 

.. . 05150000. , 

.,50555500 

•OOOOOOOO 

.00100279 

•00000501 

1 

16 

. OOOOOOOO. 

_, Q1.412^ , 

_j_ooioaoQO.. 

,05500050 

•OOOOOOOO. 

.00131509 

.00000656 

1 

17 

. .00000000 _ 

.01265206 

.00500000 

.,^55550500. 

,00000000 

•00162060 

•oooQoeio 

1 Ife 

.QQOQQnnO 

.01090722 

.00500000 

.OQQQQQQQ __.J)QQQQQQQ 

,Q01B101.7_ 

-.00000915 

1 

15 . 

...OOQOOQQO,. 

.0O90Q0aQ_ 

_.0Q100055__ 

_,05D0QQ0Q. 

.00000000 

.00193626 

•00000968 

1 

14 

.00000000 

^00709278.. 

.. ...051^0 QCM10_ 

,50000000 

•OOOOOOOO 

.00185017 

•00000925 

1 

13 

.00000000 

^00534794 

_ .,55550000... 

.,00505000 - 

,00000000 

,00162060 

.00000810 

1 

12 

.00000000 

_. 00387689 .0050 QQ£1Q. _ 

_,J1550Q500 . 

..OOOOOOOO _ 

.00131509 

•00000656 

1 

11 

..OOOOOOOO. 

._QQ27L973_ 

_..555550Q5_ 

_,50QQ55Q0 

•OOOOOOOO 

.00100279 

•00000501 

1 10 

.QPQQQQQQ 

.001 8581 3 

.00500000 

.00000000 

.00000000 

. QQQ779Q9 

.000001^6 

1 

9 

.00000000. 

.^0012 42.45 

. i5Q500000_ 

l5505QQOQ_ 

,00000000 

.0005U99 

^00000256 

1 

8 

.00000000 

.0QQ615L4iL 

.C05QQQQQ 

_,QQ50D050_ 

.OOOQQQQQ .. 

-,00035085 

.50000175 

1 

7 

.00000000 

.00052519 

,^00555550 

•OOQOQOOQ 

^QOOOQOQO 

,00023643 

•QOOOOlie 

1 

6 

.00000000 

.00033073 

.00505000 . 

. ,00000000 

•OOOOOOOO 

.00015752 

•00000079 

1 

5 

.00000000 

•0002DI65_ 

,00500000 

•OOOOOOOO 

,00000005 

,00010416 

,00000052 

1 

4 . 

.ooooaoaa_ 

.QQQ11652 

_ .QQ5QQQQQ 

..,50050005 . 

-,00000 0 00. . 

•J} 0056653 

,50000034 

1 

3 

.00000000 

.00006059 . 

_ ,00100050 

.OOOOOOOO 

•OOOOOOOO 

•00004494 

.00000022 

1 

2 

.00000000 

. .Q00Q2396. 

.. ,05155505- 

.,110005000 

. .OOOOOOOO 

,00002941 

• 000005L5 

1 

1 

^OOOOOOOO. 

..COOQDQOQ 

._,0555555Q_ 

_,5Q00Q0QC_ 

.00000000 . 

. ,00001922 

,00000010 



APPENDIX E 


Geometry Output for Lower Surface of Strut 

The geometry output for the lower surface of the strut is printed only for 
N STRUT = 1 . 


N M 

1 29 

2 29 

3 

^ 29 

5 29 

6 29 

7 29 
e 29 

9-Z9. 

10 29 
U ^9 

12 29 

13 29 
1^ 29 

15 29 

16 29 

17 29 
le 29 
19 29 

. 20 29 
.^1. 29 

22 29 

23 29 

24 29 

25 29 

_ 26 .2.9. , 
27 29 - 
._2B_29 . 
. 29 ^9 . 
- 30 29. 

. 31 29 _ 


XXI 


XETA 


YXI 


YETA 


AJ 


.C5500000 ,01000000 ,00250000 ,00000000 -,00016457 

,05750000 ,01767067 ,00250000 -.00000000 -.00032913 

. .*0600QQQJ3 . ,Q17 3417A .QQ25 QQQQ -tOUOOODQO -.00032913 

,06250000 ,01701261 ,00250000 -,00000000 -,000329X3 

,06500000 ,01668340 ,00250000 .00000000 -.00032913 

,06750000 ,01635434 ,00250000 ,00000000 -,00032913 

.07000000 ,01602521 ,00250000 ,00000000 -.00032913 

,07250000 ,01569608 .0Q250000. -.00000000 -.00032913 

OlSJiaOOi) . .. ..01536695 .0Q25QOQO .OQQQOOQQ - -OQQ32913 _ 

,07750000 ,01503702 ,0025000Q_ ,00000000 -,00032913 

,00000000 ,01.470069 ,.Qfi250QO{l tCOOCOQOO -^00032913 

,00250000 .01437956 .0025OO00-.-j00OQQ000 -.00032913 

,00500000 ,014050.43 ,002500QQ* _.ODQOOOOO -^00032913 

,00750000 ,01372129 ,00250000 .00000000 -.00032913 

.09000000 ,01339216 . .fijQ25QQilO--.^0Q£lDi)00 -.00Q32911_ 

,09250000 ,01306303 .00250000 -.00000000 -.00032913 

.09500000 .01273390 .00250000 .00000000 -.00032913 

.09750000 .01240477 .110250000 .00000000 -.00032913 

,10000000 ,01207564 . ,00250000 -.00000000 -,00002720 

,10250000 ,0123.5037 _,0i)25000.0 -, 00000000 ,00027473 

.105^10000 .0 1262510 ,00250000 .OQQQOQOQ . Q 0Q27 473 

,10750000 ,01289983 . .00250000 -.00000000 .00027473 

.11000000 ,01317456 ,00250000 -,00000000 ,00027473 

.11250000 .01344929 ,00250000 -.00000000 .00027473 

,11500000 ,01372402 .00250000 .OOOOQQOO ,00027473 

.11750000 _.013_99B25 .00250000 -.00000000 .00027473 

.12 000000 — .01427340 .00250000 -.QQOQonnn ._0Q027AJ3_ . 

^12250000- -.01.454^021 .00 250000 -.00000000 .00027473_ 

- .125HOQOO^_^Q1462293 .00250000 -.QOQQO QCQ _ ..000274 73 

.12750000 .015.09766 .00250000 .QQOQQOnO __*0QO27473 

-.13000000 . ,01537239 ,002500 00 .QOQQQOQQ jj00000557_ 


,00001922 .00000005 

.00001067 .00000005 

-.000 01 a 52 _ , Q QQQQQQ5 
.00001816 .00000005 

.00001761 .00000004 

.00001746 .00000004 

,00001711 ,00000004 

,00001676 ,00000004 

. aOUQ 16 4 1 _.Hil Q a00£L4 
,00001570 .00000004 


000015C0 

.00001430 

.00001360 

.00001269 


.00000004 

•00000004 

.00000003 

.00000003 


.00001219 . ..110000003- 
,00001149 .00000003 


00001C78 

00001006 

,00000928 


.00000003 

•00000003 

,00000002 


• 00Q010C2 .110000003 

...OHOHIO 6I__. Q00110QQ3-. 


00001131 
00001196 
•00001260 
•00001325 
.00001389 
..1100014 54 
..QQ001518 
-.00001563 
.00001612 
-.00001641 


,00000003 

,00000003 

.00000003 

•00000003 

.00000003 

.QOOOQOQ4 

.00000004 

.00000004 

,00000004 

. 1100000 ^ 4 ^ 


58 



APPENDIX E 


Flow-Field Output 


L?20700 

KAX.ERROR- 

.1383Q-E-03 

MEAN 

SORT 

ERkCR- 

•16595E-06 

L-20750 

KAX. ERROR- 

.30886E-03 

MEAN 

SORT 

ERRDR- 

.26362E-06 

L-2OB00 

MAX. ERROR- 

•15369E-03 

MEAN 

SORT 

6RR0R- 

•20901E-06 

t-20850 

MAX. ERROR- 

.29912E-03 

MEAN 

SORT 

ERR CR« 

.31551E-06 

L-20900 

MAX. ERROR- 

•25059E-03 

MEAN 

SORT 

ERRCR- 

•25991E-06 

L-20950 

MAX. ERROR- 

.29888E-03 

MEAN 

SORT 

ERRCR- 

•26I56E-06 

ND. OF ITERATIONS -20966 TIKE - 

.2635199E-03 

TIME LIMIT- 

• 2635959E-I 

03 

CCP- .-5000 

Ctl- 

- 8nnn 

EJU- 1 ,nnnn 


FOTI ■ 1 ^fmOQ 

LOCAL ERROR CPITERION- 

• 00010 

average error 

CRITERION- 

•0000001 

ORDER OF EXTRAPOLATION 

BOOY STATION. 

NX.- 27 






X y 

U 

V 

P 

T 


RQ. SH 


•0e2500 

•082500 

•082500 

•082500 

•082500 

^Qfi25-Qa 

•082500 

.082500 

.082500 

.082500 

.002500 

4.0a25QD 

•082500 

.082500 

.082500 

.082500 

.062500 

•082500 

•082500 

•082500 

.082500 

.082500 

.082500 

•082500 

.062500 

.082500 

.082500 

.082500 

.082500 

^QB?5QQ 

•Q8Z5Q0 

.062500 

.OBZ5QO 

.08Z5.00 

•06Z5Q0 

.062500 

.082500 

.082500 

.062500 

.082500 

.002500 

t082i00 

.OfiZ^iQO. 

.082500 

.082500 

t_082_500 

.082566 

tO»25QQ. 

..082500_ 

.082500 

.082500 

.062500 

.082500 

.082.500 

.082500 

• 082500 

• 082500 
•08250C 
.062500 


..oaia54 

•031839 

.031816 

.031783 

.031735 

.031667 

•Q31 56B 

.031927 

.031227 

•030996 

.030559 

.030035 

• 0Z9 39i 
.028979 
.027920 
•026216 
•029927 
•023638 
.022935 
.021381 
.020509 
.019820 
.019296 
•013906 
.018628 
.018928 
.018266 
.018188 
.018119. 
.01B07L 

•018016 
•018000 
.0X9362- 
.019339 
..Q1^9X 
.019230 
-..019X13- 
.013989 
.013773 
.013955 
.0X259 7 

• 01Z35_6 
. 1 QII 995 

•01090X 
.009102 
^00t68J| 
006269 
.009966 
■.003_872 
.003011 
. ^502370. 
.001912 

. 00159 . 9 . 

.001378 

.00X233 

.001137 

.001079 

.001032 

.001006 

.000987 


-OQQQQ 

191.93205 
336.07852 
578.07776 
8 05. 2976 9_ 
986.93322. 
nn.7?75Z. 
1290.99651 
1365.aeC7X 
15Q0. 035X7 
1587.67329: 
1653.987X0. 
1ft7Q.36R9fl 
1683.79256 
1680.99560 
1679.98531 
1679.96388 
1689.2792^ 
1676.61699 
1669.09793 
1669.07806 
1669.92529 
1663.79095 
1696.28555 
1626.9X619. 
1558.12728 
1937.55773 
125Q.X91Z9_ 
1071.71691.. 
865 .60208 

-Z67..5AQ40 
. . ..QOOQO 
7 57,99202 

.509^7X917 

850 ,75098 

1083.60635 

.1965.73013 
156X. 26636. 
1639._Q57Q9. 
1 650.56900 
1668.6 2196 
1697.68176 

X650 ..16967. 

16 92.36726 
1729^X96 
173X1.97619 
1733t70X5_5 
1736 i996 59 
1713.7916 6 
I66X . 1 8952 
1653,50732 
.1619,^101 
X_6_2 1.71661 



,00000 
-16.06998 
-36.95959 
-69.58265 
_-85. 39118 
.-109.12011 
-118,79615 
-138.72676 
-19a. 85078 
.tX 69. 13759 
-182.92320 
-1B9:.17761 
-ZQ3.Q3762 
-181.11519 
^2^19209 
-206.16029 
-189.19661 
-199.13599 
.rJX7.aai79 
-150.37669 
-183.67210 
-215.63091 
-175.97593 
-70.29511 

-X 6 . 1120 a 

-26.30672 

-X3.Q7965 

■^13.89921 

-6.03502 


1,668X1 
-..53586. 
-.QQOQO 
- 31.7 0999 
-82.29 051 
-111 .00615 
-151.65619 
-167.50022 
-XX2_.OX796_ 
t216. 95.890 
-210. 51868. 
-238.77675 
-215.22299 
-273.37145 

■^ 54,77 X33- 
^ 37.26265 
-90.02671 
-10.26960 
2.6B921 
r3iPQ226 
67.658 12 
225.9992 1 
269.26_6iC 
29X. 89582 _ 
X6X._950O2 
219.39321 
181.60771 
136. •9999 9 
100.05599 

99.80071, 

■OQOOO 


21117. 65Z77 
21117.65277 
21090.51608 
21170.87359 
21015.02696 
21266.95906 
20912.22281 
21389.75892 
20859^98179- 
2X929.67558 
21121.9Q7Q7 
21298.61539 
22328,21261 
20268.58505 
21009.27238 
22323.90896 
21900.67316 
21239.56059 
23869.62837 
26157.35799 
25003.17521 
23381.30219 
27399.37958 
37976.16622 
.9X798.91X18 
91791.89687 
93290.19975 
92510. 20Q69 
93165. 6913X 
92192._992X6 
93Q23,qQ2Dq 

92an.xax92 

92971.XaX90 

29-9IX.Xa6.99 

29612.Zn662 

29222.2575-6 

29878..68764 

23969,99326 

25199 .9iea.9 
X9005-.91005 
29981^18220 
26X91.25372- 
2X122.211X1 
2S5I8. 92230 
26 639.66 689 
1702 7.00295 
1X9 9 8 ,.90129 
10919. 19jai 
103U .91829 
9588t62L52 
X2659.7119X- 

21219^019X2. 

27X18_.9_9898 
25120.X827j9 
. 2636.9..9598X 


25970.60636 

25622.29331 

25027.66985 

25761.99902 

25761.99902 


.19X9.911526. 

.Q5XL2 

1945683,08563 

,QQQQQQ 

,nonQnn 

1939.90576 

.05112 

1995883.08563 

•OOQQQO 

,000000 

1920.07701 

.05175 

1926967.36132 

•000002 

•000000 

1350.09297 

.05932 

1369209.38656 

.000017 

•000000 

1250.97902 

.05856 

1256106.17966 

•00U069 

.000000 

1127.02993 

.06575 

1132101.56738 

.000198 

.000000 

1022,32666 

.07127 

1026927.13113 

,000385 

,nnnnnn 

905.19391 

.08239 

909216.55292 

.000995 

.000000 

779.06809 

.09307 

777551.39909 

.000505 

.000000 

621.31771 

.12015 

629113.69232 

.000626 

.000000 

509. 05603 

.19959 

511396.78013 

.000596 

•QOOOQO 

920.69120 

.17599 

922589.30916 

.000261 

.000000 

391 .97Q9Q 

,1 qaha 

X.93332xX6396- 

.0000 az 

• onoiioQ 

375.99293 

.18783 

377689.90010 

.000010 

•000000 

379.17616 

.19266 

361565. 60Q67 

.000002 

.000000 

386.89891 

.20107 

388589.72723 

•OQQOOl 

•000000 

381.93901 

.19599 

383155.98809 

.000000 

.000001 

378.12336 

.X9567 

379829.91351 

^QOOQ. 

•DO Q 001 

3.91.320Xa. 

^21259 

X930aX.0AaX6. 

-.XLOQOOQ. 

.^OOQjQO 

909.32357 

.22592 

906193.02259 

•OOQQQO 

.000001 

397.31669 

.21997 

399109.61809 

•^QOOOOO 

.000002 

390.87965 

.20892 

392633.58388 

•000002 

•000003 

907.72313 

.23368 

909557.88697 

•000017 

•000001 

998.05501 

.29532 

950071.25703 

•000199 

. 000002 

980.6289X 

.30302 

982791.29191 

-.QQOSXO 

.Qonxox 

577.99980 

.25193 

500600.79916 

•000616 

.000002 

739.09716 

•20910 

792372.86875 

.000519 

•000000 

998.35927 

•15618 

952626.86225 

•000399 

•OOOOQl 

1117.57226 

•13958 

1122601.33088 

•0Q0399 

.000002 

1279.77820 

.11151 

1285537.7008? 

•000109 

.000003 

1 4? 3, 1 0908 

,10534 1929508, Q9375. 

•000018 

•.QQQQQX 

_X98 a. 724-95- 

-xXQQSX. 

.1995973.93600 

•000001 

•000003 

.198B.7299X 

.1DQ52- 

X995973. 93608 

.000000 

•0C0003 

2X35.15858 

.Q63B7_ 

13.91X66.78876 

.000000 

.000000 

1223..139Q66_ 

.06738 

IZiaiAax 3X939. 

-.^000011 

..QOOOQO 

1X99.7159-7 

. .Q739XJXi98a9^ 69X93 

.0000.70 

-•OQQOQO 

-988.09997 

.08773 

992995.69629 

.000182 

•000000 

825.60987 

T lOp 6 

-tiZ.93?5.XX5Za 

■•000209 

,nnnnnn 

612.69228 

•13850 

635539.39513 

.000270 

•000000 

522.66908 

.16003 

-5X5021.08968 

.000179 

.000000 

939.395X6- 

.20093 

936299.71807 

.000039 

.QQOQOQ- 

912.71018- 

.21810 

9X9589.87790 

.000009 

,000000 

-326.91005- 

x20-970- 

. 398766.95803 

.000000 

.000000 

4?n,«;n?i8 

,?3631 

47?'^q4,943 67 

,000000 

..QQOOQQ 

119^20681- 

-.22910. 

j916070,7569X 

•000000 

.000000 

163^72526 


_3653_62, 02588 

•000000 

-.000000 

317x85192 

,12605 

1X9281,75610 

.000000 

,000000 

1.03.23391 

.12010 

_305 100,75092 

•a-ooooo 

,000000 

-301.61352 

.11912- ^29X0.77609 

-•DOOOOO- 

•OOODQl 

301.95676 

.11064 

30331^a363Z^ 

.000000 


336^96697 

xl3109_ 

_3i8O0X ,16151 

■QQQOQO 

-,-QOQQM 

901^21522 

■20046 

-A05X93 .36988 

■QOQOQQ 

_.GOOQOQ 

91L.83389 

,22961 

913687.19712 

•OOODOl 

•COOQOO 

.197,01618 

,22091 

_i98063. 52915. 

,00.0016. 

•DOQOQQ 

99IxX2OX0 

.20591- 

-99.9_X82,.362IQ. 

-,i^0120 

,0000.00. 

566,9 5673 

.15729 

569508.03069 

.0Q0X9X 

.QOQOQQ 

612,80505 

.11X33 

816962.67060 

•000196 

•000000 

1051.57282 

. 08990. 

10-56309.89938 

,000111 

-lOOOOOO- 

1229.50098 

.07319 

1235119.09396 

•000018 

,000000 

1325. 1U09 

.06779 

1331079.08599 

•OOOOQl 

*000000 

1325. U109 

,06779 

1331079.08599 

..QOOOQO 

,0^00000 
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Elow-Eield Oitput on Lower Surface of Strut 


The flow-field output on the lower surface of the strut is printed 

nn = 1 


NSTRUT = 1. 


j. KJl. 


xs 

YS 

US 

V5 

PS 

TS 

RQS 

• 055000 

.018000 

.00000 

.00000 

21150.10588 

1022.61893 

.07206 

.057500 

.017671 

.00000 

.00000 

26800.30404 

1204.73457 

.07753 

•060000 

.017312 

.00000 

.00000 

26062,07177 

1258.61660 

.07215 

.062500 

.017013 

.00000 

.00000 

26132.16556 

1268.71339 

.07177 

•065000 

.016683 

.00000 

.00000 

25678.79687 

1306.53963 

• 06646 

.067500 

.016354 

.00000 

•00000 

24707.57060 

1209.19170 

•06678 

.070000 

.016025 

• DQOOQ. _ . 

...*00000. 

. 25546.14460. 

1330.53078 

*06690 

.072500 

.015696 

•00000 

•00000 

25998,06986 

1330.99596 

• 06806 

.075000 

.015367 

.00000 

.00000 

25403.51061 

1315.06317 

.06731 

.077500 

.015030 

,00000. _. 

-*0.0000 

_ 25671.50292 

1347.91056 

•06636 

.080000 

.014709 _ 

.00000 

_.*OODQO _ 

— 25094. 74709 

1342.85261 

•06511 

•062500 

.014380 

-. *00000 

__ .OOOOQ 

-24473.19644 

1335.15058 

• 06387 

.085 000 . 

..014050.. 

— *00000 *00000 2565 8 *71111 . 

.13X3*160X3 

•06511 

.087500 

.013721 . 

_.OQ0OO_ 

. _*0QQQQ . 

25440.28089 

1361.34970 

•06511 

.090000 

.013392 

. ..QODQQ. 

__.DOQOO,. 

_ 25335*41168 , 

1357.54659 

•06503 

.092500 

*013063.. _ 

.QQQQQ _ 

.00000 _ 

__25D11*3A130 .. 

.1383..89I77 

.06297 

.095000 

.012734 . 

. - .qqqqq 

0.0000.. 

. .24701.53615 .. 

. 1368*0^0297. _ . 

.06292 

.097500 . 

*012405 

.QQQQQ 

.—*00000 _ 

255 2 6*2103 8-, 

.1401*97981 . 

*063-44 


135398^.1977^ 1353984.197735 


-*.lilQDniL_ *HlZ0Zij , 


.102500. 

.105000. 

.107500 

•IIQQOO 

.112500 

>115GQ O 

..H75Qa 

..1200Q0_. 

*122500 

.125000. . 


aOJ.2350 . _ _.OilOJOO ..QQQQQ ___ _ 9323 . 1 1 7fife 

.*012625 .OOQOJ3 *00000 7.95A*1250a_ 

.Q1290U-. *00000 *00000. 7233, 4?4Rft 

.013175. *00000 *lOOQOO__„.. 6951*06013 _ 

*013449 .00000 .QQQQQ O340.8322X_ 

* 0 13 7 2 4 .QQQQQ .qqqqq 


.. iaia.00272 
.1315*5 0334. 
.1334*56212. 
.33Q9*1923Q_ 
1353*51634_ 
^324*26413- 
■1353.23405 


*013393 _,*Q0QQQ *00000 104.9 4.* 36 405 _ 10X4.23447 

. 01 42X3_ .QQQQQ . 00000 14690,87377 .70312 

.014 548 *00000 *00000_. _20691*04227 1391 *.46116 . 

.014823 .._.OOQOO .OOOQQ 24326*21X14.. _1450*0QX4a.. 


5HS HS VISS 

1965 1027220.719653 .*000042 

8051 1210155.880506 .000046 

7540 1264280.575397 .000048 

Q142 1274422.601416 .000046 

5387 1312419.053874 .000049 

5706 1294993.057856 .000048 

7092 1336518. 17Q923 *000049 

<»612 1336905.446119 .000049 

•000049 
.000049 
.000049 
•000049 
*000050 
•000050 
•00005C 
•000050 
•000050 
*000051 
..*000049 
•000049 
.*000049 
^*000049- 
*000050 
*000069 

.02107 1353323*6070.7- 135 9 32 3. ^Q7Q^g 

.03402 1079671.221^6 1079671.221560 ^ *000043 

-*Q354X_14522O0.7ei34 1452208..731338 _ .000052 

.*O5X01339X722*73432.1397722*X3432I .000050 

-*059 90 14 56532* 51753_JJ» 5 653 2.5175 9.3 * 00 0052 


— *.0.4394 -1354Q63*X3L23 1354 Q66. 7312 R3 
.02469 1321423. 7055Q ^321423. 705490 
. *02077 1340668*10302-1340668*103019 
..01925^1315083*66731 1315083*667314 
- .01789 1359607.16406.1359607*164064 
. -.01668 133Q223.31745 1330223. 31 744Q 






Mass 

Flow Rate 

Output 




BODY 

station 

NO.- 

33 

X( METERS )• 

.09750 

MASS 

FLOW 

RATE( KG/M/StC )■ 

6.07562 

BODY 

STATION 

NO ■ ■ 

34 

X( METERS )?. 

.10000 

MASS 

FLOW 

RATE(KG/M/S60- 

8.05485 

BODY 

STATION 

N0.» 

35 

X( METERS)- 

.10250 

MASS 

FLOW 

RATE(KG/M/SEC)- 

8.05039 

BODY 

STATION 

ND.« 

36 

X(METERS)- 

.10500 

MASS 

FLOW 

RAT£(KG/M/SEO- 

8.05796 

BODY 

STATION 

NO* ■ 

37 

X(METERS)- 

.10750 

MASS 

FLOW 

RATE<KG/M/SEO- 

8.03889 

BODY 

STATION 

NO.- 

38 

X(METERS)- 

•11000 

MASS 

FLOW 

RATE(KG/M/SEO- 

8.05932 

BODY 

STATION 

N0.« 

39 

X(METERS )• 

.11250 

MASS 

FLOW 

RATE(KG/M/SEC )• 

6.04369 

BODY 

STATION 

NO.- 

40 

X(METERS)- 

.11500 

MASS 

FLOW 

RATE(KG/M/SEO- 

8,06730 

BODY 

STATION 

NO.- 

41 

X(METERS)- 

.11750 

MASS 

FLOW 

RATfc(KG/M/SEO- 

6.05163 

BODY 

STATION 

NO.- 

42 

X(METER5 )• 

.12000 

MASS 

FLOW 

RATE(KG/M/SEC )- 

8.07953 

BODY 

STATION 

NO.- 

43 

X(METERS)- 

.12250 

MASS 

FLOW 

RATE(KG/M/SEO- 

8.06094 

BODY 

STATION 

NO.- 

44 

X(METERS)- 

.12500 

MASS 

FLOW 

RATE(KG/M/SEO- 

0.09378 
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Flow Spillage Output 


BODY 

BODY 

BODY 

BODY 

BODY 

.aODLY. 

BODY. 

BODY,. 

BODY, 

30DY 

BODY 

-BODY 

BODY 

.BODY. 

BODY. 

BODY. 

.BODY. 

BODY 


STATION NO, 

station no.' 

STATION NO, 
STATION NO.« 
STATION. NQ.' 
SJkJJLDh. K 
STATlCm ND.*. 
STAIIDN m.! 
STATION. NCt? 
STATION. NQi 
STATION Nlh.- 
<?TATIDN Nn,» 
STATION^KIU.- 
S.TAT1QN..NQ, 
STAJ10N._NQa! 
STA TION NQ.« 
STATIDN_ KCU« _ 
STATION NO.' 


26 

XP- 

.092376 

LOCAL MASS SPILLED- 

.12079 

TOTAL 

MASS 

SPILLED 

UPTO 

XP- 

27 

XP- 

.095263 

LDCAL MASS SPILLEO- 

.12302 

TOTAL 

MASS 

SPILLED 

UPTO 

XP- 

26 

XP- 

..09.6150 

LQCAL MASS SPILLED- 

•12671 

TOTAL 

MASS 

SPILLED 

UPTQ 

XP- 

29 

XP- 

.101036 

LOCAl MASS SPILLED? 

.15362 

TOTAL 

MASS 

SPILLED 

UPTO 

XP- 

30 

X2-. 

•103923 

LOCAL MASS SPILLED- 

.17421 

TOTAL 

MASS 

SPILLED 

UPTO 

XP- 

11. 


.1 n^ftin 

IQCAL JUSS-S RULED- . 

.19476 

TOTAL 

lUSS 

SP-ILLEO 

UP.TO. 

XP- 

32, 


.i09I?.9.7 

LOCAL MASS 5PULED- 

*20270 

TOTAL 

MASS 

SPILLED 

UPTO 

XP- 

33. 


■112563 

. LOCAL MASS PUL ED- 

.22472 

TOTAL 

MASS 

SPILLED 

UPTQ 

XP- 

34. 

He. 

X115A71 

LOCAL JIAIS IP.ILLED? 

.24052 

TOTAL 

MASS 

SPILLED 

UPTQ 

XP- 

35_ 

JL1-. 

_ai8_15.t 

LOCAL MAES. S.PJL LEO- 

.±24475 

total. 

MASS 

SPILLED 

UPTO 

XP- 

36_ 

,x^ 

,rl2l24.4 

LOCAL AAS5_SEILLED? 

.24369 

total mass 

spilled 

UPTO 

XP- 

3X 

XP- 

■12413Q 

LOCAL MASS SPILLED- 

.x23U4 

total JIAES spilled UPTQ XP- 

15. 

11?_ 

■127Q17 

LOCAL MASS_3P1LLED- 

.23119 

.TOTAL. 

MASS 

EPILLEO 

UPTQ 

XP- 

19. 

XP- 

.129904. 

LOCAL MASS_.SP1LLED- 

.22547 

.TOTAL 

J1AS5 .EEILLEO 

UPIQ_.X.P- 

kO. 

H±. 

■132791 

LOCAL ilAS S. 1 PI LLED- . 

• 2204.8. 

TOTAL MASS 

.EPlLLED 

UP.TQ..XP-. 

41 

Xl±- 

a316.77. 

LOCAL MASS SPILLFH • 

.21061 

.lQIAUMA5S_EEiLL£D 

UPTQ 

XP? 

-42_. 

XP- 

.136564 

LOCAL J1A3S__1PULE0? 

. .19.84 3 

TOTAL 

.MASS 

EPILLEO 

UPTO. 

XP-. 

43 

XP- 

.141451 

LOCAL MASS SPILLED- 

•16332 

TOTAL 

MASS 

SPILLED 

UPTO 

XP- 


Pressure Ratio Output on Lower and Upper Boundaries 


H XP 


1 

..OOQOQ 

.2L 

■ .005-77 

3 

_.0U55, 

4 

. 01732. 

5 

. .C23Q9. 

6 

.±02.6.87. 

7 

.a034Q.6 

il 

.03868 

9 

.X.04272 

10 

...0,4619_ 

u 

..OA907. 

12 

.05196 

13 

.05485 

14 

.±05174. 

15 

.06.062 

16 

.06351. 

17 

±06640. 

18 

.06926, 

19 

.07217 

30 

..C.I5Q.6. 

21 

.07794 

22 

.08083 

23 

.08372 

24 

..08660, 

25 

.08949. 


CENTER LINE PR. RATIO (M = 

ItOOOQ 

1x110.70 

1.Q06L 

I.QOIZ 

KQ109 

J^OQll 

I.Q125 

1.00Q5 

-I.01A5. 

. x99I9 
.1.022.3 
• 9938 
1.0200 
c99ao 
1.0079 

1.0029 
.9973 
.99A9 

1.0030 
,x59Q0 

.9900 

.9968 

1.0^30 

2.1081 

2.6353 


1 ) SIDEWALL PR. RATIO 

.L.OQOO 

1.7feA7 

2.2301 

2.2139 

2.1625. 

2.1627 

2..i.7.68 

2.1529 

_2.1.66ii 

2.1424 

2.1016 

2.0947 

2.0974 

2.1517 

2.1632 

2.1563 

2.1223 

2.1324 

2.0776 

2.1131 

2.0951 

2.1114 

2.1145 

2.1314 

2.1406 


(M 


.00366 

.00421 

•004Sfi 

•00496 

•00546 

xOn5J99 

•00656 

•00718 

.00785 

.00855 

.00926 

_x01Q63 

.01129 

.01193, 

.Q121t_ 

.01315 

.01370 


= Ml) 
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Pressure I^tio Qitput on Lower and Upper Surfaces of Strut 

The pressure ratio output on the lower and upper surfaces of the strut is 
printed only for NSTRUT = 1. 


N 

yp 

16 

.06351 

17 

.06640 



19 

*07217 

20 

.07506 

21 

.C7794 

22 

.08083 

23 

.08372 

24 

.08660 

25 

.08949 

26 

.09238 

27 

,09526 

28 

.09615 

29 

.10104 

30 

.10392 

31 

.10681 

32 

.1C97C 

.33 

_ai258 

34 

.11547 


.11836 

36 

.12124 


STRUT UPPER SURFACE PR. 

1.4932 

1.2633 

1.1509 

1.1512 

1.1044 

1.2087 

1.2125 

1.2956 

1.9811 

3.1793 

4.2971 

4.9100 

5.0094 

4.9075 

4.6010 

4.3770 

4.1791 

3.9674 

-4.1153 

3.5351 


RATIO 


STRUT LOWER SURFACE PR. RATIO 

2.1150 

2.6608 

2.6062 

2.6132 

2.5679 

2.4708 

2.5546 

2.5998 

2.5404 

2.5672 

2.5095 

2.4473 

2.5659 

2.5440 

2.5335 

2.50U 

2.4702 

2.5526 

1.7000 

.9323 

.7954 


Program Termination Message 


-OtLJlUYliaLI TItlE. CQNVEEGENLE.CRITERliDN 
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